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Preface 


This  thesis  pertains  to  a  very  significant  problem,  the  formation  of 
secondary  lobes  in  antenna  radiation  patterns  due  to  dielectric  radomes. 

The  problem  is  not  understood  well  enough  to  even  suggest  a  sure  approach 
to  obtaining  a  solution.  This  study  is  only  one  of  several  studies  which 
could  be  made.  Yet,  it  is  meant  to  be  a  foundation  which  other  studies 
could  be  built  upon  to  eventually  engineer  a  definite  solution.  Several 
interesting  complexities  of  the  problem  are  realized  by  this  study  and  are 
discussed  to  set  that  foundation  to  work  from. 

The  symbols  used  throughout  this  paper  are  the  same  as  those  used  in 
the  literature  sighted  or  are  defined  as  they  occur.  The  primary  refer¬ 
ence  used  was  Time-Harmonic  Electromagnetic  Fields  by  R.  F.  Harrington 
(Ref.  2),  and  the  symbols  follow  almost  without  variation  from  that  text. 

It  should  be  pointed  out  that  this  thesis  contains  theory  on  series 
solution  and  moment  method  solution  formulations.  But,  most  of  the  time 
for  this  study  was  spent  learning  about  and  working  on  the  Hewlett-Packard 
21MX  M  series  minicomputer.  This  system  was  chosen  to  be  used  because  of 
its  easy  accessabil ity.  A  sincere  thanks  is  extended  to  ASD/ENAMA  for  the 
unlimited  use  of  their  system  as  well  as  for  the  use  of  other  resources 
and  also  to  the  support  from  that  branch  especially  through  William  J. 

Kent,  electronic  warfare  engineer. 

An  error  was  discovered  having  to  do  with  the  results  from  the  series 

solution  presented  herein.  The  theory  Is  correct,  but  the  calculations 

of  the  coefficient  F  is  double  what  it  should  be  for  orders  not  equal  to 

n 

zero.  This  error  should  not  change  the  results  for  the  scatterer  of  six 

ii 

_ _  J 


tenths  wavelength  outside  diameter  where  the  zero  order  term  should  domi¬ 
nate.  But,  for  a  diameter  of  sixty  wavelengths  the  error  could  become 
significant  and  should  be  checked  against  the  results  in  Appendix  £.  The 
results  from  the  moment  method  are  considered  to  be  correct. 

A  special  thanks  goes  out  to  my  advisor,  Major  August  Golden,  Jr.. 

He  provided  the  motivation  and  guidance  which  was  needed  throughout  this 
study.  And,  the  scrutinized  reading  of  this  thesis  by  Bill  Kent  and 
Captain  Thomas  Johnson  was  also  sincerely  appreciated. 

Robert  K.  Schneider 
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This  thesis  is  dedicated  to  my  wife,  Linda,  and  to  my  family  at  home. 
They  all  had  confidence  that  the  hard  work  would  pay  off  and  were  the  in¬ 
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Abstract 


The  problem  of  scattering  by  thin  cylindrical  dielectric  shells  of 
large  circular  cross  sections  is  approached  by  two  methods:  (1)  an  infi¬ 
nite  series  of  eigenfunctions,  and  (2)  the  method  of  moments.  Numerical 
results  are  presented  for  shell  radii  of  0.3x,  3. Ox,  and  30x,  the  source 
being  an  electric  line  current  near  but  external  to  the  shell.  Computer 
programs  are  presented  which  implement  these  two  solutions.  When  the  scat¬ 
tering  structure  does  become  large  limitations  on  numerical  results  are 
encountered  due  to  computer  memory  and  speed  limitations.  Other  difficul¬ 
ties  are  also  encountered  in  an  analysis  of  such  scattering  problems  and 
are  presented  and  discussed  along  with  recommendations  to  resolve  such 
difficulties. 
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I .  Introduction 

Scattering  by  perfectly  conducting  obstacles  is  exhaustively  discuss¬ 
ed  in  the  literature  (Ref.  1;  2;  3;  4;  and  5).  The  same  is  not  true  for 
scattering  by  dielectric  bodies.  Several  studies  have  been  done  dealing 
with  scattering  from  planar  dielectric  surfaces  (Ref.  1;  2;  3:  4;  5;  6; 
and  7).  Pathak  (Ref.  7)  gives  an  exhaustive  analysis  of  the  diffraction 
of  a  TMq.  Surface  wave  by  a  dielectric  slab  terminated  and  flush  mounted 
in  a  perfectly  conducting  surface.  This  thesis  is  motivated  by  the  study 
of  radomes  which  are  usually  curved  surfaces.  Studies  of  the  scattering 
by  curved  dielectric  surfaces  have  been  rather  limited.  Most  of  the  arti¬ 
cles  found  in  the  literature  concentrate  on  scattering  by  solid  dielectric 
cylinders  of  circular  or  elliptical  cross  section  (Ref.  8;  9;  10;  11;  12). 
Some  studies  have  been  done  concerning  other  curved  structures  (Ref.  13  to 
20). 

Kuester  and  Chang  (Ref.  14)  investigate  the  continuous  radiation  of 
a  wave  as  it  travels  along  a  uniformly  curved  section  of  open  waveguide. 
They  present  a  technique  for  determining  the  radiation  loss  at  a  discontinu 
ity  in  curvature.  A  very'  modest  study  of  theradome  scattering  problem 
has  been  done  modeling  with  a  cylindrical  shell.  Lewin,  Chang,  and  Kuester 
(Ref.  13:95)  examine  the  case  of  plane  wave  incidence  on  the  scattering 
surface,  which  lends  insight  to  the  problem  but  is  not  exactly  the  case 
when  the  source  is  assumed  to  be  very  close  to  the  scattering  surface. 

Thiele  (Ref.  21:306)  explains  the  theory  behind  the  method  of  moments 
as  do  Harrington  (Ref.  22)  and  Richmond  (Ref.  23).  Harrington  and  Richmond 
apply  this  method  directly  to  scattering  from  a  dielectric,  circular  cylin¬ 
drical  shell  for  small  geometries;  that  is,  cross  sections  of  less  than 
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40/u  wavelengths.  Richmond  details  the  technique  of  setting  up  the  equa¬ 
tions  to  calculate  the  scattered  field  and  discusses  employing  the  Lag- 
range  Interpolation  method  (Ref.  24)  to  reduce  the  size  of  the  matrix 
arrived  at.  Richmond's  study  only  applies  to  geometries  much  smaller  than 
the  sixty  wavelength  cross  section  (see  Appendix  A)  of  concern  in  this  re¬ 
port.  Also,  Harrington  (Ref.  2:198)  develops  the  basics  for  setting  up  a 
series  (or  eigenfunction)  solution  for  this  problem  using  cylindrical  wave 
functions . 

This  thesis  develops  an  infinite  series  of  eigenfunctions  as  a  solu¬ 
tion  for  the  scattering  of  electromagnetic  waves  incident  from  a  nearby 
electric  current  filament  upon  a  dielectric,  cylindrical  shell  of  circular 
cross  section.  A  computer  program  is  used  to  generate  numerical  results 
from  this  series  solution.  These  results  are  compared  to  Richmond's  (Ref. 
23:338)  for  the  geometry  indicated  there.  The  solution  by  the  method  of 
moments  is  also  developed  and  programmed.  Results  are  compared  to  those 
from  the  infinite  series  solution.  The  geometry  is  then  enlarged  to  ap¬ 
proach  the  sixty  wavelength  cross  section  of  interest.  The  objective  is 
to  compare  the  results  obtained  from  the  infinite  series  solution  and  those 
from  the  moment  method  solution,  and  thus  obtain  confidence  in  the  solution. 
Finally,  difficulties  encountered  from  either  solution  method  are  enumer¬ 
ated  and  discussed. 
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II.  The  Problem  and  Approach 


The  radome  and  antennas  configuration  to  be  modeled  is  shown  in  Fig¬ 
ure  1.  The  antennas  are  located  six  inches  behind  the  radome  in  the  me¬ 
tallic  fuselage  and  have  beam  maximum  at  45°  from  the  skin  of  the  air¬ 
craft.  The  resulting  antenna  radiation  pattern  is  shown  in  Figure  2. 


Figure  1.  Antenna/Radome  configuration. 


Figure  2.  Present  radiation 
pattern. 
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Figure  3.  Recommended  radiation 
pattern . 


Note  In  Figure  2  the  secondary  lobes  at  plus  and  minus  thirty  degrees, 
their  magnitude  being  only  about  ldB  down  from  the  primary  lobe.  The 
fact  that  the  patterns  overlap  off  the  zero  degree  line  is  likely  due  to 
asymmetry  in  product  manufacturing  and  installation  of  the  aircraft  an¬ 
tennas  and  radome.  The  recommended  radiation  pattern  appears  in  Figure 
3.  The  secondary  lobes  have  been  decreased  to  about  6dB  below  the  prima¬ 
ry  lobe  and  cross  over  occurs  at  zero  degrees.  A  comparison  of  the  ideal 
and  typical  patterns  radiated  appears  in  Figure  4. 


Figure  4.  Formation  of  the  secondary  lobe  from  the  left-hand 
antenna  relative  to  the  ideal  pattern. 


The  main  concern  here  is  the  cause  of  the  undesired  secondary  lobe. 


This  thesis  does  not  answer  that  question  but  investigates  an  approach 
which  might  be  useful  in  engineering  an  answer.  As  a  first  step  toward 
that  answer,  the  problem  is  modeled  as  shown  in  Figure  5. 


Radome- 


— ECM  Antenna 


« 


-Dielectric  Cylindrical 
Shell 


-Metallic  Fuselage 


Figure  5.  Modeling  the  radome  as  a  dielectric  cylindrical  shell 
and  the  ECM  antennas  as  line  sources. 


This  model  is  simply  a  dielectric,  cylindrical  shell  of  circular  cross 
section  with  an  external  line  source.  The  model  does  not  account  for  the 
metallic  fuselage,  but  it  may  still  lend  some  insight  to  how  electromagnet 
1c  radiation  scatters  from  large,  curved  structures.  In  addition,  this 
model  and  approach  may  bring  out  some  of  the  problems  which  will  be  en¬ 
countered  in  trying  to  answer  the  question  of  lobing  due  to  large  struc- 


Ill .  The  Basic  Theory 

A.  Separation  of  Variables/Eigenvalue  Method 

The  geometry  to  be  analyzed  is  shown  in  Figure  6.  It  is  assumed  that 
the  cylindrical  shell  and  the  line  source  are  infinite  in  z  and  that  there 
is  no  z-variation  in  any  quantity.  All  quantities  have  an  eJwt  time  depen¬ 
dence  assumed  and  suppressed.  Since  the  source  is  a  z-directed  electric 
current,  the  potential  due  to  the  source  is  a  vector  potential  A,  with 

A.  *  0  (7  »  0;  A  ,  A  *  0)  (1) 

z  P  4> 

Note  that  in  this  study  H  *  v  x  A  which  is  consistent  with  Harrington  (Ref. 
2).  The  current  on  the  wire  is  given  by  (also  see  Appendix  B) 

J  =  I  6(F)  Z  (2) 

It  is  assumed  that  the  permeability  in  each  region  is  that  of  free  space. 
Also,  it  is  assumed  that  the  dielectric  material  is  a  perfect  insulator, 
having  a  zero  loss  tangent. 


a 

i 

flj  ■  1 1^1 1  1 K 

/  /f/r 

k|SSq^^^B 

1  Vti.il 

mi 

_ 

Figure  6:  Model  for  the  eigenvalue  solution,  including  the 
numbered  regions. 


The  fact  that  there  is  no  variation  with  z  and  that  there  Is  only  a  nonzero 
scalar  potential  given  by  (1)  implies  that  Az  must  be  a  solution  to  the 
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scalar  Helmholts  equation  for  source-free  regions  (3),  where  the  surfaces  in 
Figure  6  are  not  in  any  region  but  present  boundary  conditions  to  those  re¬ 


gions. 


-  2  2 
7  Az  +  k  Az  “  0 


A  solution  to  (3)  in  cylindrical  coordinates  is  (Ref.  2:199) 


where 


\  ■  Bn(kp»>  ”„<♦> 


IT  *  »•  kz  '  o 


Bn(kpp)  is  some  solution  to  Bessel's  equation  of  order  n  and  hp ( 4> )  is  sin  <|>, 
cos  ,  or  some  linear  combination  thereof.  The  separation  parameter  equa¬ 
tions  are 

kP.2  +  kz2  =  ki2;  1  =  °*  !*  n»  111  (6) 

substituting  (5)  into  (6)  yields 

k  -  k.;  i  =  0,  I,  II,  III  (7) 

pi  1 

From  (1)  and  £5)  the  electric  field  intensity  can  be  written  as 
Ez  =  -h  Az 

where 

A  A 

ZQ  =  juv.  The  magnetic  field  is 

H  =  3Az  ♦  1  aAz  -  (9) 

9p  9  P 

For  the  continuous  geometry  shown  in  Figure  6,  one  would  (correctly)  expect 
the  potential  to  have  even  symmetry  about  the  origin,  hence  a  cosine  varia¬ 
tion  in  <(». 

Az  *  Bn^k1p^  cos(n$)*  n  an  integer  (10) 

From  (8)  and  (10)  the  total  field  for  a  given  region,  i,  and  azimuthal  index. 


n.  Is  given  by 

Ez  “  costn^  01) 

i  ,n 

Referring  to  Table  5.1  (Ref.  2:203),  the  potentials  for  the  respective  re¬ 
gions  can  be  written  as  follows: 


REGION 

ELECTRIC  VECTOR  POTENTIAL 

0 

■  f„"n(2!  (V)  «>(»♦) 

u,n 

(a) 

I 

Az  *  ^nW^  +  BnVkop^  cos(n^ 

I,n 

(b) 

02) 

II 

Az  =  ^nW^  +  °nVk2p^  cos  'n^ 

II,n 

(c) 

III 

A7  *  E  J  (k  p ) 

2III,n  n  n  0 

(d) 

where 

ko 

I — — — •  2irf  2v 

*  “JVo  *  “ET  = 

03) 

k2 

*  -JvJ  =  w{voeocr  =  koJ^ 

04) 

and 

H  ^  the  Hankel  function  of  order  n  of  the  second  kind 
n 

J  is  the  Bessel  function  of  the  first  kind 
n 

N  is  the  Neuman  function  (also  called  the  Bessel  function 
n 

of  the  second  kind) 

These  potentials  automatically  satisfy  appropriate  boundary  conditions  at 
the  origin  and  the  radiation  condition.  Equation  (16)  contains  six  unknown 
constants  for  a  given  order  n,  An...Fn.  To  solve  for  these  six  unknowns  the 
following  general  boundary  conditions  are  applied  (see  Figure  7): 

"  x  (H,  -  ■  3S 

(VW*"-R,  (,5) 
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Figure  7.  Geometry  for  the  boundary  conditions  over  a 
surface  S. 

Since  the  obstacle  (the  dielectric,  cylindrical  shell)  Is  not  a  magnetic  con¬ 
ductor  and  is  assumed  to  be  perfectly  Insulating,  the  magnetic  and  conduction 
current  are  neglected.  With  this  the  following  six  boundary  conditions  re¬ 
sult: 

n  x  Ei  n  8  n  x  ny  i  ■  0,  I,  II  @  p  =  a,  c,  b  respectively  (a) 

n  x  n  «  n  x  ny  1  *  I,  II  @  p  *  c,  b  respectively  (b)  (16) 

n  x  K  .  -  n  x  H.  _  ■  J- .  @  p  *  a  (c) 

o,n  l,n  s 

—  I  * 

where  Je  ■  -  6 (<t> )  z  on  the  p  =  a  surface  as  shown  in  Appendix  B.  The  goal  is 

S  d  * 

to  determine  the  far-zone  radiation  pattern  from  the  line  source  in  the  pres¬ 
ence  of  the  obstacle.  This  field  will  be  given  by  a  sum  of  the  field  given 
by  (11)  and  (12a)  over  all  possible  modes. 

Ez(po’  *o}  "  ^  I  Fn  Hn(2)  (kopo}  cos  (n*'}  (17) 

n=o 

The  observation  point  for  the  field  in  the  far  zone  is  (pQ,  $0).  From  (17) 

it  is  evident  that  the  only  coefficient  that  need  be  solved  for  is  F  .  This 

means  that  once  the  problem  is  solved  for  F  for  each  n  from  zero  to  infinity, 

n 

the  radiation  at  all  points  in  region  0  can  be  directly  obtained. 

By  applying  (16)  six  .equations  in  six  unknowns  result  (see  Appendix  B) 
which  can  be  arranged  into  the  follow  matrix  equation  : 
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"JnV> 

"r'V1  0 

0  0 

•Hn(2)'v> 

An 

0 

Vkoc) 

Nn(koc)  -Jn(k2c) 

-Nn(k2c)  0 

0 

Bn 

I 

0 

X> 

CM 

C 

T> 

O 

Vk2b>  'W* 

0 

Cn 

0 

Wkoc> 

Wkoc)  -k2Jn’(k2c> 

-k2Nn'(k2c)  0 

0 

i 

Dn 

0 

0 

0  k2Jn<k2I’> 

Wk2b>  -koVkob> 

0 

En 

0 

J  (k  a) 

N  (k  a)  0 

0  0 

-H  (2)(k  a) 

F 

I 

n  o 

n  o 

n  o 

n 

rkoaJ 

Matrix  equation  (18)  can  be  solved  to  obtain  coefficients  An-..Fn  for  a  given 


order  n. 

Richmond  (Ref.  23:338)  gives  results  for  scattering  by  the  cylindrical 
shell  in  terms  of  echo  width  per  wavelength.  To  make  a  comparison  to  his  re¬ 
sults  the  equation  for  the  normalized  echo  width  for  this  eigenvalue  solution 
is  derived.  In  a  two-dimensional  problem  having  linear  polarization  as  that 
considered  here,  the  echo  width  can  be  defined  by  (Ref  2:358) 


L  -  lim  (2irp 

6  p-KO 


I5 


) 


(19) 


The  total  field  at  some  point  in  the  far  zone  is  the  sum  of  the  incident  field 
at  that  point  due  to  sources  without  the  obstacle  present  and  the  scattered 
field  due  to  the  polarization  currents  impressed  upon  the  obstacle  by  the 
sources 


E  =  Is  +  E1  (20) 

The  scattered  field  in  (19)  is  determined  by  (20).  The  reference  field,  E^, 
in  (19)  will  be  assumed  to  be  the  incident  field  at  the  center  of  the  obstacle 
without  the  obstacle  present,  for  consistency  with  Richmond  (Ref.  23). 

The  large  argument  asymptotic  expansion  for  the  Hankel  function  of  the 
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second  kind  (Ref.  25:364) 


».(21(k.p.) 
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can  be  used  to  rewrite  the  total  field  as 


(2n+l )n/4] 


e2(p0.^0)  rr  e'J  °P°  R 1  Fne  T  cos^^0) 

A  oo  n=o 

Recall  that  in  (22)  the  far-zone  observation  point  (pQ *4>q )  is  relative 
origin  at  the  center  of  the  obstacle,  and 


where  D  is  the  diameter  of  the  obstacle. 

The  incident  field,  t * ,  can  be  written  as  (Ref.  2:236) 

H0(Z)  <Vn> 

0 

In  (24)  is  the  distance  from  the  source  to  the  observation  point  (p 
equivalentally  (X,Y),  and  is  written  as 

Pn  *  J  (X-X* )Z  +  (Y-Y')2  ' 

which  in  the  far-zone  can  be  written  as  (see  Appendix  C) 

Pn~  p  -  acos$ 

Again,  the  large  argument  asymptotic  expansion  can  be  used 


Ho‘2)  <Vn>TT 


"* a“j[k  (p-acos<|i)-ir/4] 


irko(p-acoS((>) 


o  n-H» , 

In  (32)  p>>a  which  allows  (32)  to  be  written  as 


(2) 


<Vn> 


kopn-H» 


Substituting  (28)  in  (29) 


rr 

J**? 


e’jkop  {T  ejkoacos* 


(21) 

(22) 

to  the 

(23) 

(24) 
4>),  or 

(25) 

(26) 

(27) 

(28) 
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E^(p  »<t> ) 


V1  I  2  e‘Jkop  JT  ejkoacos* 
4toe0  <  ^ 


(29) 


From  (17),  (20),  and  (29)  the  scattered  field  at  the  far-zone  point  is 


The  incident  field  at  the  center  of  the  obstacle  without  the  obstacle 
present,  E^,  is  also  given  by  (24),  but  pn  now  becomes  the  distance  from  the 
source  to  the  origin 

pn  -  W 

The  large  argument  asympototic  expansion  cannot  be  used  in  this  case  since 
the  observation  point  is  not  in  the  far-zone.  can  be  written  as 


e*  *  -y.  ty*0.>  - 1 


(31) 


where 


(2) 


(z)  =  Jn(z)  -  j  Nn(z) 


(32) 


Substituting  (JO)  and  (31)  into  (19),  dividing  by  XQ,  and  simplifying,  the 

echo  width  normalized  to  the  free  space  wavelength  can  be  written  as 
2  4  1 

Le .  _  8e  c  — n - * -  ®  . 

\  ~  Jo  (koa)  +  No  (k0a)  1  Fn  e 

tt  f  I  n=o 


+ 


ej  (koacos<j> ) 


2 


(33) 


where 

C  -  V  ?s  3  x  108  m/sec  (34) 

y  e 
0  0 

Equation  (33)  is  a  relatively  simple  representation  for  the  echo  width 
per  wavelength.  To  evaluate  the  equation  for  specific  parameters  in  princi¬ 
ple  requires  evaluation  of  an  infinite  sum  involving  constants,  Fn,  which 
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must  each  first  be  calculated  by  solution  of  (18).  The  analysis  to  this 
point  presumes  the  existance  of  an  "exact  solution"  represented  as  an  in¬ 
finite  summation.  If  this  infinite  summation  is  an  exact  representation  of 
a  sable  situation,  it  must  be  covergent  to  the  exact  solution.  If  the  sum¬ 
mation  is  absolutely  convergent,  truncation  after  a  finite  number  of  terms 
will  result  in  a  bounded  error.  The  magnitude  of  the  error  will  depend  on 
the  number  of  terms  retained.  This  method  of  evaluation  through  truncation 
is  made  use  of  in  the  program  in  Appendix  F.  The  error  produced  by  trunca¬ 
tion  will  be  addressed  later  In  this  paper. 


B.  Moment  Method  Solution 


In  the  previous  section  a  solution  for  the  scattering  by  the  dielectric 
obstacle  was  obtained  from  an  exact  solution  to  the  Helmholtz  equation.  Con¬ 
sider  now  an  inhomogeneous  region  V  containing  source  J1  and  M^  as  shown  in 
Figure  8. 


Figure  8.  Region  V  containing  sources  and  obstacles 
Within  region  V  the  following  must  hold 
-v  x  E  =  z  H  +  M i 

_  _  -  _  (35) 

V  x  H  *  yQ  E  +  J 

A  A 

Where  the  admittivity,  yQ  =  ju>e0,  and  the  impedivity,  zo  =  jwu0,  are  functions 
of  position.  If  most  of  the  region  T  is  homogeneous  except  for  some  small 
subregion(s),  Z,  region  ¥  can  be  defined  as  all  of  V  excluding  X.  Within  re- 

A  A 

gion  ¥  then  z  and  y  are  constant.  Equation  (35)  can  be  written  as 
°w  °v 

■v  x  E  =  z  H  +  M* 

°w  (36) 

7  x  H  =  y  IT+O* 

o 
w 


where 


if  =  (z  -  z  )  H  +  M1 
0  °w 

J*  -  (y0  -  i0  )  e  ♦  ? 

w 


(37) 


J6  and  M*  are  the  effective.. currents  and  can  be  treated  as  source  currents 


radiating  in  a  homogeneous  region.  For  the  problem  being  considered  z  and 


y  are  the  free  space  parameters  with  v  -  u  and  a  ^  0  so  that 
M6  =  M1 

J6  =  jio(e-eo)  E  +  J1 

The  situation  shown  in  Figure  9  is  equivalent  then  to  that  of  Figure  8. 


(38) 


Figure  9.  Equivalent  sources  radiating  in  a  homogeneous  region  V, 

The  total  field  in  region  V  is  given  by 

E  =  t '  +  E5  (39) 

The  incident  field,  E^ ,  upon  the  obstacle  (subregion  T)  is  produced  by 
outside  the  obstacle  while  the  scattered  field,  E5,  is  produced  by  polariza¬ 
tion  currents,  J*5,  induced  throughout  the  obstacle. 

J*5  =  ju)(e-e0)  E  (40) 

In  (40)  f  is  the  field  induced  throughout  the  obstacle. 

As  in  the  development  of  the  eigenvalue  solution,  the  incident  field, 
from  the  electric  current  filament  parallel  to  the  z-axis,  will  have  only  an 
z-component.  The  z-component  of  incident  field  will  produce  a  scattered 
field  having  only  an  z-component  and  meeting  all  boundary  conditions.  From 
(39)  the  total  field  will  also  have  only  an  z-component.  The  result  then 
is  an  incident  field  from  the  filament  source  (see  Figure  6)  in  all  space. 

The  obstacle  can  be  replaced  by  the  polarization  currents  induced  on  the 
obstacle  radiating  in  free  space.  The  superposition  of  the  fields  produced 


1  >L*t; 


is  the  total  field  existing  in  free  space. 

The  field  from  an  electric  current  filament  is  given  by 
2 

dE  =‘x*—  H  (2)  (kp)  dl  z 
4o)e  0 
0 

From  (40)  the  scattered  field  is 
2 


(41) 


dE 


•s  -  cr  ho(2>  "'<■> dI 


where 


and 


dl  =  JPds'  =  ju)(e-e0)  E  ds * 


2  2 
r  =  k  c 
o 


(42) 


(43) 


In  (43)  E  is  the  total  field  within  the  obstacle  which  has  only  a  z  -compo¬ 
nent,  and  ds'  is  the  increment  of  surface  area  on  the  cross  section  of  the 
obstacle.  (The  prime  on  ds'  indicates  source  coordinates.)  Integrating  (41) 
over  the  surface  of  the  obstacle  and  substituting  (43)  for  dl,  the  scattered 


field  from  a  dielectric  of  low  loss  tangent  is 
2 

EzS(x,y)  =  //s,  ju(e-e0)  Ez(x\y’)  (kp)  dx’dy' 


(44) 


ik2 


ff%»  (er-l)  Ez(x\y')  Hq(2)  (kp)  dx'dy' 


where 


Gr  ■  U. 


p  -|(^x')2  +  ly-y'  )2 

and  (x,y)  is  the  observation  point  while  (x',y’)  is  the  source  point. 


(45) 

(46) 


From  (39)  and  (44) 
i, 


Ez(x,y)  =  Ez'(x,y)  +  (-  4*_)  ff%,  (yl)  Ez(x',y')  HQu;  (k0p)  dx'dy'  (47) 
The  procedure  now  is  to  divide  the  cross  section  of  the  obstacel  into  N 


(2) 
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'square'  cells  small  enough  so  that  e  and  the  field  intensity  are  essentially 
constant  over  each  cell  (refer  to  Figure  10).  Enforcing  (44)  at  the  center 
of  each  cell  m  results  in 

Ezm  +  4-"  ,(V’>  ez„  "cell  n  Ho<2)  <C>  d*'d*'  -  Ez'  <48> 

m  n=  l  n  m 


Ez  -MvV  (51) 


and  Ez  is  the  field  intensity  at  the  center  of  each  cell  n. 
n 


Figure  10.  Dividing  the  obstacle  into  N  nearly  square  cells. 

Equation  (48)  states  that  the  total  field  at  some  cell  m  is  due  to  the 
incident  field  at  that  cell  plus  the  scattered  field  from  all  N  cells  (in¬ 
cluding  cell  m)  at  cell  m.  Yet,  (48)  is  one  equation  with  N  unknowns.  To 
solve  for  these  N  unknown  (48)  is  enforced  at  the  center  of  all  N  cells,  re¬ 
sulting  in  N  equations  in  N  unknowns.  This  process  of  enforcing  the  integral 
equation  (44)  at  the  centers  of  N  cells  is  called  Point-Matching  and  is  a 
special  case  of  the  method  of  moments  (Ref.  21:312). 

The  evaluation  of  (48)  involves  an  integration  over  a  square  surface  of 
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the  Hankel  function  of  the  second  kind  of  order  zero.  A  singularity  exists 
when  the  argument  of  the  Hankel  function  approaches  zero,  which  occurs  when 
the  source  point  approaches  the  observation  point.  This  singularity  is  in¬ 
tegrate  (see  Appendix  D)  if  the  square  surface  is  replaced  by  a  circular 

surface  of  equal  area  with  radius  a  .  The  result  is 

n 

2ir  9p  i  n\ 

U  'o  Ho  <V>  p’  d»'  d+' 

cell  n 

-  <j/2)  I>koanHi(2)(koan)  -  2j],  m  =  n  (52) 


2  VkoV  Ho 


m  =  n 


where  p  is  given  by  (49)  and 


The  polar  coordinates  p'  and  are  based  on  a  coordinate  origin  at  the  center 
of  cell  n,  the  source  point  contributing  to  the  total  scattered  field  at  cell 


m. 

For  the  geometry  shown  in  Figure  6  the  matrix  equation  given  by  (47) 
will  be  symmetric  applying  (52)  for  m  =  1,2,3,..,N.  It  is  obvious  that  the 
diagonal  terms,  m  =  n,  will  be  equal  and  the  off-diagonal  terms  depend  only 
upon  the  distance  between  the  source  and  observation  point.  The  symmetry 
Inherent  in  this  geometry  results  in  a  Toeplitz  matrix  (Ref .  21  :340) .  This 
specialization  to  a  Toeplitz  matrix  is  very  advantageous  when  inversion  of 
large  matrices  is  necessary  and  will  be  discussed  later  in  this  paper. 

Once  the  matrix  equation  has  been  solved  for  the  field  induced  within 
the  dielectric  (44)  can  be  used  to  obtain  the  scattered  field  in  the  far  field. 
Applying  the  method  of  point  matching,  (44)  can  be  rewritten  as 
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(54) 


E2S(*,y)  -  -£o  ^  U„  -  1)  EnV>,  Ck0.„)  H0<2>(Vf,) 

7 

Note  that  in  (54)  the  fact  that  the  observation  point  is  in  the  far  field 
guarantees  that  it  will  be  different  from  the  source  point  on  the  obstacle  so 
that  the  appropriate  integral  for  m  /  i  of  (52)  is  used  and 

“n  *  l<X  -  */  +  (V  -  (55) 

From  Appendix  C  the  distance  from  the  source  point  on  the  obstacle, 

(Xn , Yn ) ,  to  the  observation  point  in  the  far  field,  (X,Y),  can  be  written  as 

pn  =  po  "  Xn  cos*‘  Yn  sin*  (56) 

where  pQ  and  $  are  the  polar  coordinates  of  the  observation  point  in  the  far 
field.  Because  the  observation  point  is  in  the  far  field,  the  large  argument 
asymptotic  expansion  for  the  Hankel  function  of  order  zero  (Ref.  25:364)  given 
by  (57)  can  be  used. 

Ho(2)(kp)  -*{^7  e”^kp~ir/4^  (57) 

Since  pQ»  Xn  and  pQ»  Yn,  the  result  is 

H0(2)(k0p0)  -jiffs;  Vin*>  •  (58) 

Substituting  (58)  into  (54). 


EzS(X,Y)~-j]2p0  e“^kopo  i  (V1)Wl(koan) 

n=1  (59) 

«eJ£k0(XnC0S*  +  Y^s1n<j» )  +  tt/4] 

As  was  done  for  the  eigenvalue  solution,  the  echo  width  per  wavelength  can  be 
determined  using  (59)  and  (19),  which  is  rewritten  here  for  convenience 


lim  (2irp  .  ) 

.  O  r-l 
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The  incident  field  referred  to  in  (19)  is  a  constant  field  intensity  to  which 
the  scattered  field  is  normalized.  So  that  comparison  of  results  can  be  made 
to  Richmond  (Ref.  23:338),  this  incident  field  is  taken  to  be  that  at  the  cen¬ 
ter  of  the  dielectric  shell  with  the  obstacle  removed,  radiated  from  the  elec¬ 
tric  filament  line  source.  This  incident  electric  field  is  given  by  Harring¬ 
ton  (Ref.  2:236)  and  is 

E,1  -  -  Jil  H  {^}(k  |p-p 1 1 )  (60) 


where  | p-p '  |  is  the  absolute  distance  from  the  source  to  the  center  of  the 
obstacle  and  is  relatively  small  compared  to  pQ.  Therefore,  (19)  is  written 


Le/  -j  -ik  o  E  (e  -1)E  a  J  (k  a  leJ^ko^Xncc>s,,,+YnS1n,,’^  2 

=  lim  (2trp  J  Tp  e  JKopo  V  n  utnanJrka  n,e  o  n  n  ) 

o  o  _  n= i  _ 

JlI  H  (2)<k  |p-p'|) 


A  "  ,)e3[ko(Vos*+Ynsin*>’ 

o  n=i  n  n  n  I  o  n 


-V  Ho(2>(k0|p-p'|) 

4ujeo  (61 ) 

It  is  also  possible  to  obtain  the  total  field  in  the  far  field  by  adding  (60), 

A  A 

evaluated  at  p  =  pQp  +  4>p ,  and  (58).  It  should  be  noted  that  the  far  field  Is 
determined  using 


o  X 


where  D  is  the  diameter  of  the  obstacle  and  x  is  the  free  space  wavelength. 


IV.  Procedure  to  Arrive  at  Numerical  Results 


A.  The  Eigenvalue  Solution 

The  theory  and  derivation  behind  both  the  eigenvalue  solution  and  the 
moment  method  solution  are  not  overly  complex,  but  several  complexities  were 
realized  in  attempting  to  obtain  numerical  results  from  these  solutions. 

The  eigenvalue  solution  involves  an  infinite  series  of  coefficients  multi¬ 
plying  Hankel  functions.  The  coefficients  are  arrived  at  by  solving  six 
simultaneous  equations.  Each  equation  is  comprised  of  Bessel  functions  of 
the  first  and  second  kind,  Hankel  functions  of  the  second  kind,  and  the  de¬ 
rivatives  of  each.  How  to  calculate  the  Bessel  and  Hankel  functions,  how  to 
perform  the  matrix  inversion  and  how  to  work  with  the  infinite  series  were 
the  three  major  questions  to  be  answered  to  obtain  numerical  results  from 
the  eigenvalue  solution. 

Calculating  the  Bessel  and  Hankel  Functions 

Evaluation  of  the  eigenvalue  solution  requires  calculation  of  Bessel 
functions  of  integer  order  and  of  widely  varying  argument.  An  argument  as 
small  as  1.57,  Involved  in  obtaining  results  to  compare  to  Richmond's  results 
(Ref.  23:338),  had  to  be  handled  as  well  as  those  as  large  as  190,  to  pro¬ 
duce  answers  for  a  scatterer  of  60  wavelengths  cross  section.  Another  con¬ 
cern,  in  addition  to  range,  was  how  to  calculate  the  value  of  these  functions 
efficiently.  The  matrix  equation  (23)  would  theoretically  have  to  be  solved 
an  infinite  number  of  times  which  meant  evaluating  the  elements  of  the  matrix 
an  infinite  number  of  times.  Easily,  a  routine  to  calculate  these  functions 
could  consume  too  much  CPU  time  to  approximate  the  summation.  A  third  area 
of  concern  was  for  the  accuracy  in  evaluating  the  functions. 


A  routine  based  on  equations  given  by  Abramowitz  (Ref.  25:369)  was  first 
obtained.  This  routine  involved  calculating  Bessel  functions  of  the  first 
and  second  kind  of  order  zero  and  one  using  polynomial  approximations  appro¬ 
priate  to  the  argument.  Using  recurrence  relations,  values  for  the  functions 
at  other  orders  were  calculated.  The  values  from  this  routine  were  the  same 
to  an  average  of  6  decimal  places  as  those  of  Table  9.4  (Ref.  25:407)  for 
arguments  less  than  3.  When  the  argument  was  increased  to  100  accuracy  was 
maintained  at  an  average  of  6  decimal  places  except  at  the  point  where  the 
order  also  became  100.  From  an  order  of  100  to  about  104  the  accuracy 
dropped  to  an  average  of  4  decimal  places  for  an  argument  of  99.95  (compari¬ 
sons  made  with  those  given  by  Aiken  (Ref.  26)).  Then  the  accuracy  returned 
to  an  average  of  6  decimal  places  for  order  greater  than  104.  The  program 
was  checked  to  determine  why  such  a  lack  of  accuracy  existed  for  the  case  of 
large  argument  and  order  being  equal.  The  program  followed  exactly  from  the 
equations  given  by  Abramowitz  and  no  programming  errors  were  detected.  One 
change  which  was  made  was  to  increase  the  point  at  which  recurison  downward 
began  for  evaluation  of  the  Bessel  function  of  the  first  kind.  The  program 
had  been  designed  to  begin  recurison  from  an  order  10  greater  than  the  order 
desired.  That  starting  point  was  increased  to  20,  40,  and  80  above  the 
order  desired.  No  change  at  all  was  detected  by  increasing  even  to  the  40 
point,  but  the  change  to  the  80  point  resulted  in  the  development  of  signifi¬ 
cant  errors.  The  errors  were  seen  as  exponents  of  -23  multiplying  the  value 
of  the  function  for  argument  and  order  being  large  and  equal.  This  error 
could  have  been  the  result  of  the  normal ization  value  (Ref.  25:385)  be¬ 
coming  very  large  and  exceeding  the  limit  of  the  machine  when  divided  into 
the  trial  values  the  result  was  very  small  numbers  of  exponent  -23. 

The  problem  with  the  previous  Bessel  routine  was  not  investigated 
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further  mainly  because  it  was  not  an  economical  routine  to  use.  It  would 
return  values  of  Jn(x)  and  Nn(x)  for  only  one  order  at  a  time  and  would 
therefore  have  to  be  called  some  2100  times  for  the  scatterer  of  sixty  wave¬ 
lengths.  Another  routine  was  obtained;  it  calculated  the  zero  and  first 
order  terms  in  the  same  way  as  the  previous  routine  did,  but  it  used  what  is 
called  "continued  fraction  formulae"  (Ref.  27:153).  This  approach  results 
directly  from  the  recurrence  relations,  and  it  produced  no  more  accuracy  than 
the  preceeding  routine  except  at  equal  large  argument  and  larger  order.  At 
that  point  accuracy  was  maintained  at  an  average  of  4  to  5  decimal  places. 

The  attractive  feature  of  this  routine  was  that  it  would  pass  an  array  of 
values  starting  from  the  order  desired  down  to  order  zero,  resulting  in 
fewer  required  calls.  Yet,  because  of  the  array  size  and  the  number  of  ar¬ 
rays  required,  core  limitations  were  exceeded  on  the  machine  used.  Another 
problem  with  this  routine  was  that  it  would  not  return  correct  values  for 
the  Bessel  function  of  the  second  kind  having  arguments  less  than  17.5. 

This  problem  was  corrected,  but  because  the  routine  actually  consisted  of 
four  separate  Subroutines  and  a  lot  of  core  was  required,  the  first  routine 
was  used  to  try  to  obtain  some  type  of  answers  for  this  eigenvalue  solution. 

Matrix  Inversion 

With  the  matrix  formed  it  now  had  to  be  inverted  to  solve  for  the  co¬ 
efficients.  There  are  several  methods  available  in  the  literature  to  ac¬ 
complish  matrix  inversion;  Gaussian  Elimination,  Gaus-Seidel,  and  Gaus- 
Jordan  are  three  typical  methods  (Ref.  24,  25).  But,  for  matrices  of  order 
n>2  a  method  known  as  Cholesky's  method  (also  named  Crout  reduction)  "requires 
fewer  arithmetic  operations"  than  either  of  the  three  mentioned  methods, 
"making  it  the  fastest  of  the  basic  elimination  methods"  (Ref.  24:198).  It  Is 
a  lower-upper  triangular  decomposition  method  and  can  be  made  economical  of 
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core  by  storing  the  lower  and  upper  triangular  matrices  in  the  same  area 
allocated  to  the  matrix  to  be  inverted.  This  method  results  in  a  simple, 
fast  program  for  matrix  inversion. 

The  Infinite  Series 

As  discussed  in  the  development  of  the  eigenvalue  solution,  the  Infinite 
summation  must  converge  to  the  exact  answer.  In  order  to  determine  the  ac¬ 
curacy  of  the  truncated  summation,  the  behavior  of  the  coefficient  for  in¬ 
creasing  order  must  be  considered.  A  definite  bound  on  the  error  due  to 
truncation  was  not  determined,  yet  the  computer  generated  numbers  for  the 
coefficient  Fn  did  show  that  Fn  became  purely  imaginary  and  then  approached 
zero  for  increasing  order,  as  was  predicted.  Approximately  350  terms  were 
needed  in  the  sum  before  additional  terms  became  insignificant  relative  to 
some  e  for  a  structure  of  60  wavelength  diameter.  It  was  opted  to  proceed 
with  the  programming,  working  with  the  truncated  infinite  series  as  de¬ 
scribed,  in  order  to  obtain  some  type  of  results. 

B.  The  Moment  Method  Solution 

The  most  serious  obstacle  in  obtaining  numerical  results  from  the  mo¬ 
ment  method  solution  was  the  demand  for  core.  The  matrix  to  be  inverted  be¬ 
came  quite  large  for  a  large  geometry  since  the  number  of  cells  must  increase 
to  maintain  accuracy  In  the  final  answer.  Also,  to  keep  the  program  as  gener¬ 
al  as  possible  so  that  characteristics  of  the  obstacle,  such  as  size  and  shape, 
might  be  changed,  the  program  was  written  as  a  collection  of  separate  sub¬ 
routines.  This  programming  method  called  for  other  arrays  in  addition  to  the 
matrix  to  be  set  up  which  quickly  exceeded  the  amount  of  core  available  to 
the  program.  Fitting  the  program  in  the  available  core  and  obtaining 
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reasonable  execution  time  were  two  pressing  problems  in  programming  the 
moment  method  solution. 

Core  Requirements 

In  developing  the  theory  for  the  moment  method  solution,  it  was  men¬ 
tioned  that  the  symmetry  of  this  particular  problem  would  yield  a  Toeplitz 
matrix.  Numbers  were  generated  for  a  simple  case  and  showed  that  indeed  the 
matrix  was  Toeplitz  in  form.  This  was  a  significant  observation  since  rou¬ 
tines  exist  which  take  advantage  of  this  special  form.  The  routine  chosen 
to  be  used  (see  Appendix  G)  does  not  require  the  passing  of  the  entire  matrix 
to  be  inverted  but  only  the  first  row  or  first  column  of  the  matrix.  Over¬ 
all  the  routine  requires  only  five  arrays  of  dimension  n  where  n  is  the  di¬ 
mension  of  the  square  matrix.  This  results  in  a  significant  reduction  in 
the  demand  for  core. 

To  combat  the  core  requirements  due  to  programming  by  use  of  subroutines 
a  feature  called  Extended  Memory  Area  (EMA)  was  utilized.  All  arrays  were 
first  defined  to  be  in  a  common  block;  that  common  block  was  then  defined  to 
be  in  EMA.  Basically,  EMA  is  disk  memory  as  opposed  to  core  memory.  A  com¬ 
piled  program  is  loaded  into  core  to  be  run,  and  arrays  which  are  not  speci¬ 
fied  to  be  in  EMA  are  allocated  space  in  core.  Only  so  much  core  can  be 
taken  up  by  any  particular  program,  64K  for  the  HP  21MX  M  series;  if  more 
core  is  asked  for,  the  loader  aborts  and  the  program  cannot  be  run.  EMA  al¬ 
lows  the  program  to  be  loaded  and  run,  but  it  also  has  limitations. 

Processing  Time 

EMA,  in  conjunction  with  designating  a  large  background  when  loading, 
will  allow  a  large  program  to  be  loaded  and  run,  but  it  also  requires  more 


processing  time.  The  system  has  to  swap  data  in  EMA  from  disk  to  core  in 
order  to  process  the  data  as  directed  by  the  main  program.  This  swapping 
takes  up  much  more  computer  time  than  working  with  data  already  in  core. 

For  the  size  of  the  arrays  required  to  analyze  a  large  structure  of  60  wave¬ 
length  cross  section  run  time  exceeded  two  hours.  Obviously,  this  was  a 
costly  program  to  run  on  this  minicomputer  as  programmed. 

C.  The  Priority  Set  for  Both  Methods 

The  intent  was  still  to  produce  two  independent  computer  solutions  for 
the  scattering  from  a  dielectric  cylinder  of  circular  cross  section  and  to 
form  an  argument  for  their  validity.  In  an  attempt  to  meet  this  goal  the 
preceding  concerns  were  incorporated  into  two  programs.  Not  all  questions 
were  answered  nor  were  all  the  concerns  met.  A  discussion  of  the  results 
from  these  two  programs  follows  along  with  conclusions  and  recommendations 
on  how  to  improve  upon  them. 
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V.  Discussion  of  Results 

Computer  generated  output  appears  in  Appendix  E,  including  plots  of  the 
scattered  and  total  field  for  a  number  of  trials.  As  emphasized  throughout 
this  paper,  the  goal  has  been  to  generate  numerical  results  from  the  two 
methods  discussed,  to  compare  those  results  to  Richmond's  results  (Ref.  23: 
338),  and  to  increase  the  size  of  the  structure  to  model  an  actual  radome. 

In  the  process  an  argument  for  the  validity  or  invalidity  of  the  answers  plus 
a  discussion  of  problems  encountered  would  be  presented.  Looking  at  Figures 
E-2,  E-3,  and  E-4  along  with  Table  E-l ,  it  is  evident  that  the  results  are 
very  close  to  being  the  same  for  the  small  structure  of  0.6  wavelength  cross 
section.  The  fact  that  equivalent  results  were  obtained  from  the  two  com¬ 
pletely  different  methods  seems  to  indicate  that  the  error  due  to  the  trunca¬ 
tion  of  the  sum  in  the  eigenvalue  solution  was  not  a  significant  determinant 
of  the  final  answer.  Both  programs  calculated  the  Bessel  functions  using 
the  same  routine  with  the  accuracy  problems  discussed  previously.  The  re¬ 
sults  thus  obtain  credibility  because  these  two  methods,  which  were  indepen¬ 
dent  and  very  different,  resulted  in  comparable  answers. 

Before  ever  attempting  any  calculations  one  would  probably  assume  that 
a  small  obstacle  of  0.6  wavelength  cross  section  would  cause  very  little 
disturbance  of  the  field  incident  upon  it.  Figures  E-5  and  E-6  are  consis¬ 
tent  with  this  assumption,  showing  a  very  smooth,  slight  variation  over  360°. 
This  structure  does  not  possess  any  abrupt  changes  in  curvature;  therefore, 
a  wave  which  might  be  trapped  by  the  structure  would  radiate  consistently  as 
it  traveled  around  the  structure.  As  the  trapped  wave  radiated  it  would  be 
losing  energy  and  would  in  turn  have  less  and  less  to  radiate.  What  differ¬ 
ence  a  trapped  wave  would  present  is  not  evident  in  this  trial,  but  this 
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argument  is  consistent  with  the  smooth  variation  of  the  total  field.  It  is 
also  possible  to  view  this  problem  as  a  superposition  of  a  circular  array  of 
line  dipoles  and  a  line  source  radiating  in  free  space.  This  is  a  physical 
picture  of  the  moment  method  solution.  Two  dipoles  having  some  current  dis¬ 
tribution  on  them  and  spaced  a  distance  d  apart  will  produce  nulls  in  their 
radiation  pattern  about  X/d  radians  apart.  The  maximum  distance  d  for  this 
first  trial  is  2c  =  0.6a.  Therefore,  nulls  should  appear  in  the  scattered 
field  spaced  by  approximately  1.67  radians  or  95.5°.  Looking  at  Figure  E-3, 
the  nulls  appear  to  be  spaced  by  about  110°. 

The  outer  diameter  of  the  obstacle  was  increased  10  fold  for  the  second 
trial.  The  thickness  of  the  dielectric  was  kept  at  0.05a  and  the  distance 
from  the  outside  of  the  obstacle  to  the  location  of  the  line  source  was  kept 
at  0.2a.  Since  the  first  trial  compared  closely  using  37  cells  to  approxi¬ 
mate  the  cylindrical  shell,  370  cells  was  used  in  the  second  trial  to  main¬ 
tain  close  to  the  same  cell  density  per  wavelength.  A  comparison  of  Figures 
E-7  and  E-8  and  Figures  E-9  and  E-10  shows  equivalent  results  existing  again 
for  the  two  different  methods.  A  couple  of  slight  differences  do  exist  though. 
First,  the  magnitude  of  the  results  from  the  moment  method  solution  are  no- 
ticably  less  than  those  from  the  eigenvalue  solution.  Round  off  errors,  in¬ 
correct  evaluation  of  the  Bessel  functions,  and/or  a  nonoptimum  number  of 
cells  could  account  for  this  small  difference.  The  second  difference  is  that 
the  total  field  of  Figure  E-9  shows  a  sudden  drop  at  about  152.5°  and  then 
a  sudden  increase  about  2.5°  later.  This  also  occurs  at  about  207.5°.  Fig¬ 
ure  E-10  shows  a  smooth  variation  in  the  total  field  in  these  regions.  The 
eigenvalue  results  do  not  show  a  squiggle  in  the  pattern  at  172.5°  and  187.5°, 
as  do  the  moment  method  results.  Both  polar  plots  were  plotted  at  every  1° 
so  the  difference  is  not  due  to  resolution. 
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An  interesting  observation  can  be  made  by  overlaying  Figures  E-5  and 
E-9  and  Figures  E-6  and  E-10  as  shown  in  Figures  E-ll  and  E-12  respectively. 
After  increasing  the  cylinder  diameter  the  result  appears  to  be  that  the 
general  shape  of  the  total  field  pattern  remains  the  same.  Effectively,  the 
circular  array  spoken  of  previously  has  been  increased  in  size  and  more  di¬ 
poles  added.  The  nulls  should  now  be  spaced  about  Vo.Ox  =  0.167  radians  or 
9.55°  apart.  The  addition  of  dipoles  to  the  array  would  cause  the  variation 
in  the  gain  pattern  as  seen.  With  the  observations  pointed  out,  it  is  felt 
that  the  results  of  this  trial  are  also  correct. 

The  third  trial  was  for  another  10  fold  increase  in  the  geometry.  The 
results  from  the  eigenvalue  solution  appear  in  Figures  E-13  and  E-14.  It 
was  not  possible  to  run  the  moment  method  program  because  the  neeo  for  3700 
cells  over  extended  EMA.  An  attempt  was  made  to  break  the  obstacle  into 
3700  cells  and  to  work  with  on  every  other  cell  in  computing  the  total  field. 
The  program  with  this  modification  would  then  load  and  run,  but  after  three 
hours  of  run  time  the  program  was  aborted.  (This  moment  method  program  slow¬ 
ed  the  system  tlown  considerably.)  Because  there  is  nothing  to  compare  the 
results  of  the  eigenvalue  solution  against,  only  general  comments  on  those 
results  can  be  made. 

Applying  the  array  analysis,  the  nulls  should  appear  at  A/60x  =  0.0167 
radians  or  0.955°.  Observing  Figure  E-12  there  are  indeed  20  peaks  in  20°. 

A  change  in  the  scattered  field  as  compared  to  the  previous  two  trials  is 
the  maximum  occuring  at  zero  degrees.  This  maximum  is  plotted  at  about  0.95 
in  Figure  E-13.  Having  normalized  all  calculations  of  the  scattered  field 
relative  to  the  maximum  magnitude  of  scattered  field,  the  maximum  as  plotted 
should  be  1.0  as  for  trials  one  and  two.  This  evident  error  immediately 
sheds  doubt  on  the  validity  of  the  results.  Another  concern  is  that  the 
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scattered  field  pattern  has  such  deep  nulls  and  so  many  of  them.  It  is  hard 
to  predict  what  is  causing  such  a  spiked  pattern  to  exist.  Symmetry  still 
holds  for  the  total  field  in  Figure  E-14,  but  again  it  is  difficult  to  vali¬ 
date  the  pattern  with  nothing  else  to  compare  to. 


VI.  Conclusions  and  Recommendations 


The  fact  that  the  two  different  programs  produced  results  comparable  to 
published  results  for  the  small  structure  of  0.6  wavelength  diameter  gives 
confidence  that  these  results  are  correct,  thus  validating  both  the  approach 
and  the  computer  program.  When  the  diameter  was  increased  to  6.0  wave¬ 
lengths,  the  eigenvalue  and  moment  method  results  were  again  nearly  the  same, 
except  for  the  small  but  noticeable  difference  in  magnitudes.  Adding  some 
56  numbers,  each  from  a  combination  of  Bessel  functions  accurate  to  only  5 
decimal  places  could  contribute  significantly  to  this  difference.  To  deter¬ 
mine  the  degree  to  which  the  results  are  correct,  one  recommended  approach 
would  be  to  bound  the  error  due  to  truncation.  A  second  approach  would  be 
to  run  these  programs  on  a  system  which  has  a  longer  word  length,  so  that 
the  machine  generated  errors  will  be  reduced.  An  additional  reason  for  going 
to  a  larger  machine  would  hopefully  be  to  decrease  the  run  time  required  by 
each  program  and  to  have  more  memory  to  work  with.  More  accurate  evaluation 
of  the  Bessel  functions  could  also  lead  to  better  results. 

The  two  methods  studied  do  present  and  address  some  of  the  complexities 
encountered  in  trying  to  solve  the  problem.  A  better  understanding  of  this 
problem  could  come  from  also  plotting  the  magnitude  and  phase  of  the  inci¬ 
dent  fields,  the  field  in  the  dielectric,  and  the  scattered  field.  This 
would  give  an  indication  of  the  interactions  of  the  fields  producing  the  to¬ 
tal  field  observed.  Yet,  the  fact  remains  that  the  point  matching  method  for 
the  structure  of  60  wavelength  diameter  required  more  memory  than  is  available 
to  the  average  system  user.  Some  method  for  retaining  the  data  generated  such 
as  on  magnetic  tape  is  needed  to  circumvent  the  core  limitations.  The  other 
methods  for  reducing  the  size  of  the  matrix  such  as  Lagrange  Interpolation 
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can  be  tested. 

With  results  for  scattering  from  the  circular  obstacle  of  6.0  wavelength 
diameter  the  scatterer  could  then  be  modeled  as  an  ellipse  and  verified  by 
an  essentricity  of  zero.  The  essentricity  could  then  be  allowed  to  approach 
one,  which  would  be  a  model  more  representative  of  the  radome  of  concern. 

The  methods  developed  to  reduce  the  size  of  the  matrix  for  the  circular  ob¬ 
stacle  could  possibly  be  manipulated  for  use  with  the  ellipse.  The  ellipse 
will  present  its  own  problems  such  as  how  the  cell  size  should  vary  around 
the  structure  and  what  effect  a  very  close  source  (relative  to  the  size  of 
the  obstacle)  will  have.  But,  these  steps  should  lead  to  an  eventual  solu¬ 
tion  to  the  secondary  lobe  seen  in  the  measured  antenna  patterns. 
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Appendix  A.  The  Model  and  Its  Parameters 


To  model  the  radome  as  a  dielectric  cylinder  of  circular  cross  section 
the  radius  of  curvature  of  the  radome  must  be  determined.  Figure  A-l  shows 
the  geometry  for  the  model  and  for  the  approach  to  be  taken. 


Figure  A-l.  Model  and  geometry  to  determine 

the  radius  of  curvature  of  the  radome. 


A  result  for  the  radius  of  curvature  can  be  obtained  using  the  following  pro 
cedure. 


Figure  A-2.  The  geometry  to  be  used  in  solving  for 

the  radius  of  curvature  of  a  circular  arc. 


The  equation  for  the  arc  S  from  p.  to  p,  is 


The  equation  for  vector  A  which  is  the  perpendicular  bisector  of  the  arc 


from  P1  to  P2  is 


.  X2”X-j 
ya  -  Y2-Yl 


xA  +  b 


(A-2) 


The  midpoint  of  segment  I  (B  to  p2)  is 


X  = 


VX2 


Y1+Y2 

Y  --V- 


Substituting  (A-3)  into  (A-2),  (X,Y)  being  the  point  at  which  vector 
segment  I  intersect,  yields 


Y-j+Y 
~2 


,t,2  _x2-x  X  +X 

W  =  -?—L  (-Lpl)  +  b 


VY1 


(A-3) 
A  and 

(A-4) 


Solving  (A-4)  for  b  results  in 


b  = 


Yl+Y2 


V  _  Y  ' 
*2 


2  +  2(Y2-Y1) 


(A-5) 


Replacing  b  in  (A-2)  by  (A-5),  the  equation  for  vector  A  is 

Y  _  X2"X1  y  Yl+Y2  X22'X12 

YA  -  -  Y2-Yl  XA  2  2^) 


(A— 6 ) 


Similarly,  the  equation  for  vector  D  is 

"2~  +  2(Y,-Y,) 


v  X3"X2  ¥  Y2+Y3 

YD  ‘  ’  Y,-Y„  XD  + 


2  2 
X  -X  c 
*3  *2 


'3~'2  u  *■  ^ v  *3“ ’2 J 

Solving  (A-6)  and  (A-7)  simultaneously  for  the  point  of  intersection 

XA 


(A-7) 


(h,k). 


(Ya  =  Yd  =  k,  XA  =  XD  =  h),  yields 


YX1  h+  VY2  X22-Xl2 

y  y  M  '  r*  ' 

Y2~ '1 


T~  +  2(Y2-Yl) 


X3'X2  „  Y2+Y3 

Y3_Y2  +  ~T~  *  2(V3-V2) 


or 


h  = 


VV3  .  X32-X22 

~T~  +  2(Y3-Y2) 


VT2  X2V  , 


X2'X1 


ly3-y2  -  y2-Yi 
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and  simplifying, 

T  [~(Y3-Yi)(Y3-Y2)(Y2-Y1)-t-(X32-x22)(Y2-Y1)-(X22-x12)(Y3-Y2) 

2  (x3-x2)(y2-y1 )-(x2-x1 )(y3-y2) 

Knowing  h  from  (A-8),  k  can  be  calculated  from  (E-6)  where 

X2-Xi  VY?  X22-Xl2 
k  =  ‘  h  +  2'('Y2-Y1 )  ^A"9^ 

The  equation  for  the  radius  of  a  circle  is 

r  =  [(Xrh)2  +  (Yrk)2]1/2  (A-10) 

Figure  A-l  shows  the  model  of  the  radome  as  a  circular  cylinder.  The 
radius  of  curvature  of  the  cylinder  can  be  arrived  at  using  (A-8),  (A-9), 
and  (A-10). 


Appendix  B.  Derivation  of  the  Eigenvalue  Solution  Matrix 


Given  the  following  representations  of  the  vector  electric  potentials 
in  the  indicated  regions 

REGION  ELECTRIC  VECTOR  POTENTIAL 


0 

Az  " 

o,n 

FnHn^^  {k0p)  608 

(a) 

I 

I  ,n 

tAnVkop'+BnVkop^  cos(n,)>) 

(b) 

II 

Az  = 

zII,n 

[CnJn(k2p)+DnNn(k2p)]  cos(n«j>) 

(c) 

III 

Az  = 

zIII,n 

EnJn^kop^  C0%M 

(d) 

n  x  Hn  -  n  x  HT  =  J 
0,n  I ,n  s 


(c) 


Ei,n  -  'I./'  \  „ 


(B.l) 


the  following  boundary  conditions  (B.C.'s) 

n  x  ^  n  =  n  x  ^(i+1)  n  ’  ®  p  =  a.b.c  respectively  (a) 

n  x  ti.  n  *  n  x  ^(i+1)  n  ’  ^s0,I,II  ®  P  *  c*b  respectively  (b)  (B .2) 

@  p  =  a 

and  the  fact  that 


(B.3) 


a  matrix  can  be  formed  to  solve  for  the  six  unknowns  A  . ..F„  Assume  no  vari- 

n  n . 

ation  with  ('j!''5'  °)- 


From  B.C.  (B.2.a)  and  (B.l)  three  equations  result 

F„Hn(2!<lV»  "  Vn<V>  +  Wkoa> 

A„W>  *  WkoC>  *  Vn<k2c>  +  Wk2C> 

Wk2b>  *  DnW>  ‘  E„JAb) 


(B.4) 
(E.5) 
(B  .6) 
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For  the  case  at  hand 


H  =  -  JL—  [— l  *  I  _ *  p] 
Ju>y  Sp  p  3$ 


5 

5P 


1  —L 

P  54> 


(B.7) 


In  both  (B.s.b)  and  (B.2.C)  n  =  p  so  that 

SA_  . 


n  x  H  =  - 


3p 


(B.8) 


From  B.C.  (B.2.b),  (B.l),  and  (B.8)  two  more  equations  result 

koCVn'<koC)+BA'(koC)]  *  Wn'(k2C)+yV(k2C»  fB-9> 

k2f\V<k2b>+BA'<k2b»  ’  koEnJn''kob)  (B-10) 

The  current  density  is  given  by 

(B.ll ) 

where  I  is  the  magnitude  of  the  source  current  and  6(F)  is  the  Kronecker  delta 
function.  The  total  current  is  obtained  by  integrating  over  the  cross  sec¬ 
tional  area  through  which  the  current  density  passes 


J  =  Jz  2  =  I6(r)  z 


!Total  3  fs  J  ds  =  1 

Integrating  (B.ll)  over  a  circular  cross  section 
/s  6(F)  ds  =  1 


(B.12) 


=  ff  6{r)  pdpdij, 


which  yields 


At  p 


6(F)  =  1  6(p)  «(p) 


Jz  =  5(P-a)«U) 


(B.l 3) 

(B  .14) 

(B.15) 


where  6(p-a)  is  the  Dirac  delta  function. 

From  B.C.  (B.2.C),  (B.ll),  (B.8),  (B.7),  and 
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(B.16) 


Mi 


E  =  Z  -  juy  A 
i  n=-°°  i 


(Ref.  1:198) 

-k„  "  {FA'<2>(M)-CVn'(koa)  +  B„Nn(ko>,11"s("*)  =  J. 


o  „  1  n  n 

n=-“ 


(B.17) 


The  surface  current  J  is  obtained  by  integrating  (B.15)  over  the  surface 

sz 

of  the  source. 

Jc  -  f  p  6 (p-a )  6(41)  dp 


I 


*t6(*> 


(B.i8) 

The  surface  current  can  also  be  expanded  in  a  complete  orthonormal  basis  set 

(B.19) 

Equating  (B.18)  and  (B.19) 


J  =  -k  Z  an  cos(n4>) 
5z  n=-® 


'ko  ”  an  cos  ("♦)  =  a 
n=-« 


or. 


I  00 

-  F"a  =  1  an  cos(n<^ 
0  n=-“ 


(B .20) 


To  apply  the  rules  of  orthogonality,  multiply  both  sides  of  (B.20)  by  cos(m <(>) 
where  m,  like  n,  is  an  integer  and  integrate  over  one  entire  interval  0  to  2tt. 


f^n  -  ttt  $U)  d<t>  *  /?Tr  "  an  cos(n<(>)  cos(m<j>)  d$ 
0  V  0  n—  n 


(B.21 ) 


which  reduces  to 


T  o 

-  •£-£-  =  /^[am  cos(m<j>)  cos(ma)  +  a_m  cos(-m<|>)  cos(m4>)  d<j> 

0 

=  7ra„  +  7ra„,m#o  (B.22) 

m  -m 

For  a  =  a  .  as  is  the  case, 
m  -m 

a  -L 

am  27TkQa 
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(B.23) 


If  m  =  o,  (B.22)  becomes 


rs  ■  2"ao 

0 


or, 


ao  2  irak 


(B.24) 


Therefore, 


I 


an  2  irak. 


for  all  n 


(B.25) 


The  complete  orthonormal  basis  set  could  also  have  been  written  as 

Js  =  “ko  *  bn  cos^n<*>) 
z  n=o 


(B.26) 


since  cos(n<|>)  is  an  even  function  of  In  the  same  way  as  an  was  determined, 
bn  can  be  found  to  be 


h  -  Ie„ 
bn  n 


2 irk  a 
o 


where 


Y1  *  n  = 

En  ""1^2,  n  > 


=  o 
o 


(B.27) 


(B.28) 


Substituting  (B.19)  and  (B.25)  into  (B.17),  the  result  is 
-K  "  (^H„l{2)(koa)  -  [AnJn*(koa)  +  BnNn’  (kQa)]}  cos(n*) 


o  '  n  n 
n=-«> 


I 


■  -  ko  2  •  5FT  cos<"*> 


(B.29) 


n*-«°  o' 

Looking  at  each  term  of  (B.29)  separately  yields  the  sixth  equation  with  which 
to  solve  for  the  unknown  coefficients  and  is 


A„OJ(M)  +  B„N„*(krta)  -  F„Hn{2),(kna)  =  - 

o 


n  n  'o'  n  n  '  o  '  n  n  '  o 
where  (B.30)  is  based  on  summing  from  minus  infinity  to  infinity.  If  the 


(B.30) 
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same  procedure  is  followed  using  { B . 26 } ,  (B.27),  and  (B.28),  the  equation 
becomes 


Vn'<V>  +  W<koa> 


FnHn(2)'<V> 


Ie„ 

2tt k  ae 


I 


2 irk  a 
o 


(B.31) 


where  (B.31)  is  based  on  summing  from  zero  to  infinity.  The  six  equations 
derived  can  be  combtned  to  form  one  matrix  equation  which  can  be  solved  for 


the  six  coefficients  A  . ..F  . 

n  n 


°n(koa) 

Vkoa>  0 

0  0  -H  ^(k  a) 

n  o  ' 

A 

n 

Jn<koc> 

N„(k0c)  -0„<k2‘> 

-Nn(k2c)  0  0 

Bn 

0 

o  Jn(k2b) 

Vk2b>  "W5  0 

X 

Cn 

__ 

koV<koc)  koV(koc>  -k2Jn'<k2c> 

-k2Nn'(k2c)  0  0 

Dn 

0 

0  k2°n'(k2b) 

k2V<k2b>  -koJn'!kob)  0 

En 

Jn'(koa) 

V(koa>  0 

0  0  -H  (2)'(k  a) 

n  o 

Fn 

0 

0 

0 

0 

0 

-I 


2irk  a 


(3.3  2) 


42 


r  Appendix  C.  Representing  the  Distance  to  the  Far  Zone  Observation  Point 


The  distance  from  the  source  point  to  the  observation  point  in  the  far 
zone,  pn>  is  shown  in  Figure  C-l .  The  stipulation  that  must  be  made  is  that 
the  distance  to  the  observation  point,  pQ,  is  much  much  greater  than  the  dis¬ 
tance  to  the  source  point,  p.  With  this  condition  met,  it  can  be  assumed 
that  for  an  observation  point  out  at  infinity  the  vectors  pQ  and  pn  become 
almost  parallel  and  therefore  p"  p-  '  (see  Figure  C-l.). 


Figure  C-l.  Geometry  for  the  distance  to  the  observation  point 
in  the  far  field  from  the  source  point. 


For  r3  parallel  to  pQ  angle  e2  =  <f>  and  e3  -  90  -  <f>  which  implies  that  04 
■  4.  The  geometry  is  such  that  r4  =  r,  and 

rl  =  r4  =  Yn  sin$  (C-l) 

Also, 


r2  =  Xn  cos<(> 

Using  (C-l)  and  (C-2),  it  is  obvious  that 


(C-2) 


Appendix  D.  Integration  of  the  Hankel  Function  over  a  Circular  Area 

Using  the  moment  method  to  evaluate  the  total  field  due  to  a  dielectric 
cylindrical  shell  or  circular  cross  section  in  the  presence  of  a  parallel 
filament  line  source,  the  integral  tn  CD-I )  must  be  evaluated  (following 
Richmond  (Ref.  23:336)). 

j k2/.  //  H  (2)(kp)  p'  dp*  d**  (D-l ) 

4  cell  n  0 

where 

p  =  | p-p  '  |  (D-2) 

see  Figure  D-l . 


X 


Figure  D-l.  Division  of  the  obstacle  into  circular  cells. 

A  singularity  exists  in  evaluating  (D-l)  when  the  argument  becomes  zero 
or  when  the  source  and  observation  point  are  one  and  the  same.  To  evaluate 
(D-l)  for  this  singularity  the  integration  becomes 

jk %  f2J  HQ(2)(kp)  p'  dp1  d*'  (D-3) 
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where 


H0(2)(kp)  =  Jo(kp)  '  j  N0(kp)*  p  =  p' 
for  a  reference  point  at  the  center  of  cell  n.  Substituting  (D-4)  into 

jk%  ;o  Ho(2)(kp,)  p<  dp’  d»'  =  {flV  fl  P,J0(kp')dp‘  d<0’ 

-  j  /q71  /“  p'N0(kp')dp'd* 

Two  equations  from  Abramowitz  (Ref.  25*484)  11.3.20  and  11.3.24  are 
t“  Ju.,(t)dt  =  zuJuCz) 

rl  tU  *  2\<*> 

Splitting  (D-5)  into  two  parts  and  substituting  (D-6)  and  (D-7)  in,  the 
ing  steps  result: 

u  =  1 

/2lT  !l  p'  Jo{kp,)dp'  d*'  *  27T  ;o  p'  Vkp')dp' 

Let 


for 


t  =  kp 1  dt  =  kdp' 


p'  =  o  t  =  o 

p '  =  a  t  =  ka 

2ir  /*  p'  O0(kp')dp‘  =  /ka  t  J0(t)dt 

Ic 

=  — (ka)  J,  (ka) 
k^  1 


And,  as  above 

flV  P1  Nn(kP')dp'  d*’  =  ^  /ka  t  N  (t)dt 
0  0  0  0  0 

-  %  (ka  N,  (ka)  +  Lilly 
vr  1  * 
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\ 


(D-4) 

(D-3) 

(D-5) 

(D-6) 

(D-7) 

follow- 


(D-8) 


(D-9) 


m mmm 


Substituting  (0-8)  and  (D-9)  into  (D-5),  the  result  is 

jk2/4  i%(ka)  J^ka)]  .  J  %(ka)  N^ka)  +  £]} 
k  k 

which  reduces  to 

j/2  (irka  J-j(ka)  -  j  irka  (ka )  -  j  2} 
or 

jkZ/4  /ZlT  /J  Ho(2)(kp,)p'  dp’  d<t> '  =  j/2(irka  H^^kaJ-j  2)  (0-10) 

for 

m  =  n. 
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4>  (  Peg'ees ) 


Figure  E-2.  Scattering  pattern  of  a  circular  dielectric  shell 
in  the  presence  of  a  nearby  parallel  line  source 
calculated  with  the  in-egral-equation  technique. 


Reproduced  frcn  Richmond  (Ref.  23:338) 
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TABLE  E-I . 

Comparison  of  the  Results  from  the  Eigenvalue  and  Moment  Method 
Programs  to  Richmond's  Results. 


For  the  three  sets  of  results: 


Permeability 

Dielectric  Constant 

4.0° 

Loss  Tangent 

0.0 

Inner  Radius 

0.25X 

Outer  Radius 

0.30X 

Line  Source  Distance 

from  Center 

0.5X 

Echo  Width/Wavelength 

Angle 

Ri chmond 

Eigenvalue 

Moment  Method 
(37  cells) 

0 

0.284 

0.282 

0.289 

70 

0.697 

0.671 

0.720 

110 

0.026 

0.022 

0.022 

180 

4.90 

_ * _ 

4.91 

5.11 

NOTE:  Richmond  (Ref.  23:338)  placed  the  line  source  at  180°  while  the  source 
was  at  zero  degrees  for  the  computer  programs.  Flipping  Richmond's 
results  over  allows  the  comparison  above. 
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electric  field  using  the  eigenvalue  solution  method. 


ANGLE  C  DEGREES  > 

Scattering  pattern  using  the  eigenvalue  solution  method 


270 

SCALE:  MAX=40  db 

electric  field  using  moment  method  solution. 


Overlaying  Figure  E-5  on  Figure  E-9  for  the  total  field  using  the  eigenvalue 
solution  method. 


electric  field  using  the  eigenvalue  solution  method. 


Appendix  F.  Computer  Program  for  Eigenvalue  Solution  Method 


This  appendix  contains  the  computer  program  which  implements  the  eigen¬ 
value  solution  method.  The  main  program,  RADMI,  the  matrix  inversion  routine, 
CLSKY,  and  the  routine  to  calculate  the  value  of  Bessel  functions  of  the  first 
and  second  kind  and  their  derivatives,  BESEL,  are  written  in  FORTRAN  4.  The 
computer  system  used  was  the  Hewlett-Packard  21MX  M  series  minicomputer. 


o  o  o  o  o  o  o  o  oooooooooooooo  o  o  c*  o  o  o  o 


II  »-  0  0004  IS  f»N 


CP  0  0  tV^ *r»  UC*  J  NU  0  0  02S  t'L!  : K=  0  0  U  M 


FTri-*,L 

CWBLk ,  0  > 

C  •**  4:  '*■  ’♦':*»-■ «:*;*  1 *:«:«;*  Me:*.*:*  +  *.* ^ r^c :*f: 
PROGPhM  KRDrii 


C 

c 

c 

c 

c 

C 

C 

C 

C 

C 

c 


<  -f  :#  :* .*  t  :*  :«.*  :f  :f  :*i :« :f  :*  :* 


ROfctfcRT  K.  SChNEI  DER 


*  THIS  PROGRAM  EVALUATES  THE  EIGENVAi  U£  SOLT ION  FOR  THE  Rh! 
OUTSIDE  AH'  I  HE  1 1*  I  TEL's  LONG  Oi  ElECi  ft' I C  CYLINDRICAL  SHElL 
RADIUS  B  AND  OUTER  RADIUS  C.  The  INCIDENT  FIELD  HAS  aS 
AN  INFINITE  LINE  CURRENT  CLOSE  TO,  BUT  EXTERNAL  TO,  THE 
AND  AT  AN  ANGLE  PHI -PRIMS  EOUAL  TO  ZERO, 


'I  A 

OF 

IT 

CY 


0  FIELD 
mMER 
‘I  OURCE 
rt£?£R 


+  h  MATRIX  EIGENVALUE  SOLUTION  HaS  BEEN  FORMULATED  LEAVING  THE  COEFFI¬ 
CIENTS.  OF  THE  FIELD  TO  BE  SOL-vED  FOR,  hN  ALGORITHM  B <'  PRESCOTT  0. 
CRGUT  IS  USED  TO  SOLVE  FOR  THESE  COEFFICIENTS . 

*  THE  FAR  FIELD  IS  GIVEN  BY  : 

£2 on  =  -jwuFn  TIMES'  THE  HANKEL  FUNCTION  OF  THE  2nd  kT.vD 
OF  ARGUEMENT  K0*RaU  AND  ORDER  n. 


THE  TOTAL  FIELD  IS  A  SUM  ON  n  OF  THE  ABOVE.  A  CONVERGENCE  OF  THI 
SUM  IS  DETERMINED  BY  THIS  PROGRAM.  ALSO,  SINCE  THE  ARGUEMENT  OF  T-i 
ABOVE  HANKEL  FUNCTION  IS  »  n,  THE  APPROPRIATE  ASYMPTOTIC  EXPANSION 
I S  USED . 


*  EPS'LN  IS  A  VALUE  INPUT  BY  THE  USER  FOR  USE  IN  DETERMINING  CCNVERGENC 

*  K 0  =  WAVE  NUMBER  IN  FREE  SPACE 

*  KB  =  WAVE  NUMBER  IN  DIELECTRIC 

INTEGER  ESTIME 

REAL  K d , K2 , MAX , NRMFLD , HP HI , J 0 

COMPLEX  HN20A, HNP20A,  SUME2, A<  6,  7  >,  F<  6  FN<  I  001  ),OI  1  ,012. 

I  SUMM , SUMP , SUM , £2 ON , D5 

DIMENSION  IBUFRC  ISO,  J2<  1  CO 

COMMON  /CMBLK/  A ,  F ,  ECHWPUK  36  06  >,  DG'RENC,  SCRPEZ'I  36  00  >, 

I  DBPRP  T  *;  3b  0  0  > 

DATA  P I /3 . 1 4 ! 59265/ 

DATA  XCR/4 . 875/ , VCR/3 . 875/ , RDS/2 . 743/ 

*  SET  UP  CALL  FOR  INPUT  FROM  DISK  FILE  "DATA I  ” 

DATA  IBUFR/2+0,2HDA,2HTA,2H1  ,3+0,2216,7*0/ 

CALL  S'POPNC  I BUFR ,  I SLU  ) 

CALL  EXEC/ 22,  1 
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O  <j  •  J  o  o  o 


mm 


i'.Crtl'-.  Jill),*;  .  1.. ,  I  .  K,iif  SLN,  J-KfeV  .  ftKr,  L'uKtlv.  .  fc.h5J.NK, 

i  .  ESTIrsE,  N,  I  SKIP 

SET  DEFAULT  FOR  ESTIME 

I  FC  EST I  HE  .  EC/ .  0  >  EST I  ME=f  C«  0 
ZERu=  0 . 0 

CALCULATE  THE  WAVE  NUMBERS 

• 

K  0=2  .  0*P  I  *  FRECi/3  .  ES 
f;2=K  O-eS&RTC  EPS  lNR  > 

COMPUTE  ARCUEfiENTS  FOR  BESEL  FUNCTION  CALLS 

X I  ~K 0«AA 

N2=K 

X3*K0*C 

X4=K2+B 

X?=K2:*C 

LOOP  ESTIME  TIMES  TO  DETERMINE  THE  EIGENVALUE  En  FOR  EACH  ORDER  n 


DO  100  t— ESTIME,  ESTIME 
K=1>  t 

EPSLNH= 1 . 0 
LH1 =L-1 
LM1  =L 


CALL  BESELC LM1 , X 1 , BB , Y , BP , VP > 
HN20A-CMPLXC  BB,-Y  > 

HNP2  0A=CrtPL  XC BP, -YP > 

AC  1 ,  1  >*CMPLXCBB,2ER0> 

AC  1 , 2  >*CMPLX< V , ZERO  > 

AC  1 , 3  >=<  0 .  ,  0  .  > 

AC  1 , 4  >=C  0  .  ..  0  .  > 

AC  1 , 5  >=<  0.  ,  0.  > 

AC  1 , C >=-HN20A 
AC  1  ,  ?  >=C  0 .  ,  0 .  > 

AC B ,  1  >=CMPLXC  BP, ZERO  > 

AC  6 , 2 >=CMRLXC YP , ZERO  > 

AC  r> ,  3  >=<  0 .  ,  0 .  > 

AC  0 , 4  >=C  0  .  ,  0  .  > 

AC  6 , 5  :>=c  0 .  ,  0 .  > 

AC  6 , S  >=-HHP 2 OA 

IFCLH1  , EQ ,  CO  EPSLNN=0.5 
•EOLMT*<  CURENT*EPSLNN  >/<  PI=*AA*-k  0  > 
AC  * , 7 >=CMPLXC  EGLMT , ZERO  > 


CALL  BESELC  LM 1 , X2 , BB , Y , BP , YP  > 
BU?M=-K0*BP 
AC  3 , 5 >=CMPLX< -BB , ZERO  > 

AC  5 , 5  >=CMPLXC  DUM 1 . ZERO  > 


ce 


Cl 


Cf.l  L  BE  SELi  LM  1  ,  X.3 .  BB ,  V ,  E'P ,  VP  > 
DUM 1  »K  0-«S?P 
OUM2-K  0*7  P 

Hi  2,  i  >=CHPLX<  BB ,  2ERG  > 

Hi  2 ,  2  >=CHPLX<  Y ..  ZERO  > 

Hi  4 , 1  >=CNPi.N<  BUM 1  ,  ZERO  > 

Hi  4 , 2  >=CMPLX\  DUM2,  ZERO  > 

CHLL  BESELi  LM  i  ,  X4 ,  BB ,  V ,  BP ,  VP  > 
DUril  =K2;*BP 
DUM2*K2*VP 
Hi  3 ,  1  >~\  0 .  ,  Ci .  > 

Hi  3 , 2  >=■  0 .  ,  0 .  > 

Hi  3 , 3  >=CMPlX<  BB  ,  ZERO  > 

Hi  3, 4  ’*=CMPlN<  V ,  ZERO  > 

Hi  3 , 6  >-i  0  .  ..  0  .  > 

Hi  3, 7 >= .  Ci.  ,  0.  > 

Hi  5,  i.>=<  0  .  ..  0  .  > 

Hi  5  2  >=i  Ci 0.  > 

t-  ~ ,  3  >«C  MF'LNi  DUM 1 ,  ZER 0  > 

>  3,  !  >«CMPLX<  0UM2, ZERO ) 

«•.  .  -•  >«<  0 .  ,  0 .  > 

Hi  5 ,  7  >~<  C< .  ,  0 .  > 


cc 


CHL L  BESELi LM 1 , X5 , BB , V , BP , VP > 

OUM1 *-K2*BP 
DUM2=-K2  *  YP 
Hi  2 ..  3  >=CNF'LXi  -BE ,  ZERO  > 

Hi  2 , 4  >=CMPLX\ -V , ZERO > 

Hi  2 , 5  )=i  0  .  ,  C> .  > 

Hi  2,  B  >=<  Ci .  ,  0  .  > 

Hi  2 , 7  >=i  Ci .  ,  0  ,  ) 

Hit  4 , 3  >=‘CMc'LXi  DUM  1  ,  ZERO  > 

Hi  4 , 4  >=CMPLX<  0UM2,  ZERO  .> 

Hi  4 , 5  >=i  6.  ,  0.  ) 

Hi  4 , 6  >=<  0 .  ,  0  .  ) 

Hi  4 ,  7  >=<  0 .  ,  0  .  > 

HOW  THHT  THE  MATRIX  H  IS  FORMED,  CHLL  CL SKY  TO  SOLVE  FOR  THE 
EIGENVALUES. 

M=N«-1 

ChLL  CL  SLY  i  H ,  M  > 


c: 

c  * 

C: 

c 

c 

IS 

CIS 

c 

c 

17 

c 

C: 


PULL  OUT  Fn  FROM  HRRHY  F 

IFi K  .HE.  ICO  GO  TO  17 
DO  I  S  1  =  1,  B 

WRITE-.'  6,  IS)  Fi  I  ),L,  I 

FCiF.  MAT  i  IX,  "F=  ”,  J  E 1 5 . 8, 2X,  E 1 5 ,  S,  5X,  14, 5X,  14,  /> 

CONTINUE 

WRITEi B,  70) 

K=0 

FHiL+ESTIME+1  >=F<6> 

DO  BO  h- 1 , B 
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{  * 

t  ' 


C  WR I TS*.  6 .  5  0  >  F  <  M  > ,  M ,  L 

50  FORMAT'  1;<:,  “F»“,2X,  1  E  I  5  .  S,  2.*;,  E  i  5  .  S  .  ~.v,  1-4 .  5.-'  ,  T-4  -> 

C6  0  CON;  1  HUE 

C  WRITE*;  6, 70  > 

70  F  ORMaT  1 X ,  //  > 

C 

100  CONTINUE 

C 

WRITE';  1 ,20  CO 

2 00  FORMAT1.'  IX,  "FN  EIGENVALUES.  CALCULATED  .  "  ; 

C 

C  *  CALCULATE  THE  ECHO  WIDTH  .PER  WAVELENGTH  RELATIVE  TO  THE  INCIDENT 
C  FIELD  AT  THE  CENTER  OF  THE  OBSTACLE 
C 


SFLDMX* 0 . 0 
TFLDMX=0 . 0 
I AVNTM=  0 . 0 
I HAX= 0 

NT  I MES—  I F I  )<•;  36  0 .  /DGREHC  > 

DO  700  I  =  1  ,  NT  1  AES 
I M 1  =  I  - 1 

RH I  =  I  fi  1  «DGRENC*PI  /  ISO. 

MCNVRG= 0 
SUMEZ*< 0 . , 0 . > 

SUMM*<  0 ,  ,  0 .  > 

SUMP=< 0 . , 0 .  ) 

C 

DO  55  0 ■ J=  1  , EST I ME 
UM 1*0-1 

ARGR= UM i *P I /2 . 

1  D 1  =STN<  ARGP  > 

2  02=C0S‘.  ARGP  ') 

1 F< UM 1  . EQ .  G>  GO  TO  540 

ARGM=-ARGP 

3  D3=S1N<  ARGM  > 

4  D4=CGS< ARGM  > 

5  SUMM=FN>;  UN  1  >*CGS<  -  UN  1  '*PH  I  >*CMPLX<  D3 ,  -04  > 

54  0  SUMP=FN<  UM 1  -t  EST  I ME+ 1  >*CGS<  JM 1  •*PH  I  >+CMPLX<  0 1  ,  -D2  > 

C 


SUME2=SUME3  +2  .  *P  I  •*FREG*PEftM*'C  3UMN+SUMP  > 

C: 

6  01  0=D2'*AI.r'1AG’:  FN<  UM1  +ESTIHE+  1  >  >*COS<  JM1*PHI  >*2.  ‘*PI'*FREQ*PEPM 

IF<D10  .  GT,  EF"SLN>  GO  TO  550 
HCNVRG-NCNVRG+ 1 
J2<  NCNVPG  .)=  J 

IF*;  NCNVPG  ,LT.  ICO  GO  TO  550 
DO  425  L=2,  10  . 

IF'.'  02 a.  >  .  NE  .  J2r  L - 1  >+1  )  GO  TO  427 
425  CONTINUE 

GO  TO  5 1 0 
427  NCNVPG- 1 0-L 

DO  426  L= 1 , 1 0 

IF<  L  ..GT .  NCNVPG >  GO  TO  420 
J2< L  )= JZ< 1 0-NCNVPG+L > 

GO  TO  426 

425  J2<  L  )=  0 

426  CONTINUE 

550  CONTINUE 

C 

WRITE<6,506> 


66 


c 

c 

Cl  00 Cl 

c 


FORMAT*:  IX..  "TRIED  MOPE  THriM  ESTIME  TERMS  IN  BUM  .“,/■> 
WRITE<£,52©>  1,0 

FO&Tt  AT •;  1  X  2X  ,14, 3/ ,  “  J-  “  ,  2X ,  1 4  ,  /  > 

IF*;  0  .CT.  IhnX>  IMAX=0 
IAVNTM=IAVNTrt+J 
IF*;  I  SHIP  .ME  i  >  CO  TO  530 

ARGC3'=<  <  K 0*R  >-*.'  PI/4 . 0  >  .>/*:  2 0 0*P I  > 

CO  TO  ? 

ARGCS=*:  K  0*R  >-<  P  I  /4 . 0  > 

• 

UR  I TE<  8 , 7  S  >  Ar.GC'r; 

FORMAT*;  1 X,  “ARGOi-  “,EI3.S,/> 

D 1 =COS<  akGCS  > 

WRITE*;  t- ,  77  > 

FORMAT*.' IX, “GETTING  HERE  ?",/> 

02= -8  I  M<  ARGCS  > 

WRITE*.  6,70) 

FORMAT-;  IX,  “GETTING  HERE  0“  ,/> 

03*SGRT *;  2/<  P I  •*K  0**R  >  > 

D5*£>.3*CMPLX<  D 1  ,  £>2  > 

E2 6N=P5*SUMEZ 

SCRPEZ*:  I  >=CABS<  ETON) 

IF*.  SCRPE3*;  I  ;>  .  GT  .  TFLDMX  >  TFLDMX=SCRP£Z<  I  > 

ARG=K  0-«ABS*:  AA  > 

CALL  BESEL  <  0 ,  ARC ,  BB ,  V ,  BP ,  VP  > 

*  JO=BB*BB 
YO=V*Y 

l>  1  =3  .  £8*  *4  ,  *8 .  854E-1 2**2 .  *8 .  /<  FRECi**2  . *  CURE  NT**  *2  .  *P  1**3  .  > 
ARG=KO*AA*COS*:  PHI  > 

D2=C0S*:  ARC  > 

C>3=S1H<  ARG  > 

C>4=2 .  *Pr*FREC/*CURENT/*:  4  .  *3  .  £8**2  .  *8 . 854E- 1 2  > 
t>  1  2=D4* CMPLX<  D2 , 03  > 
t>1 1  “SUME2+D 1  2 
xmag=cabs*:di  i  >*cabs<di  i  :> 

ECHWPLK  I  >=D t  *XMAG/*;'  JO+YO> 

WR ITE*.'  8,  1  000  >  ECHWPWt  I  > ,  PH  I 

FORMAT*;  IX,  “ECHWRW=",2X,  IE15.3,5X,  “PH  1  =  “  ,  2X ,  F? .  3,  /  > 

IF*.'  ECHWPUK  I  )  .GT.  SFLDMX  )  SFLDMX=ECHWPW<  I  > 

CONTINUE 

SCRPE2< NT IMES+ I >=SCRPEZ< 1  > 

WR I TE<  6 , 1200) 

F ORMAT ■;  I  X ,  ////  > 


wRITfc'.  1  ,  > 

FORMAT *:  "PLOT  THE  SCATTERED  F 1  ElD :‘fi  Y.'N“  > 
READ-.'  t  ,  I  2  >  I*MS 
IF*.' JAMS  .EG.  !H  >  GO  TO  505 
IF< I AMS  .HE.  1HY>  GO  TO  50? 

PLOT  THE  ECHO  WIDTH  PEP  WAVELMG7H  VS  hHgLe 


CALL  ECO^T < SFLDMX > 

DO  900  1=1,.  NT  I  ME ’S 
PHI**:  1-1  .  >*DGR£NC 
WRITE*:  6,  S00>  ECHWPW<  I  >,PHX 
FORMAT*.  1  X ..  "  ECHWPW*  “  ,  2X ,  1  £  1 5 
CONTINUE 


X  .  11  PH  I  =  “  ,  2X  ,  P  ?  .  3  , 


WPI TE< 1 ,29) 

FORMAT*: ‘‘PLOT  THE  TOTAL  FIELDS*  Y/N"  > 

READ-:  1,12)  I  AH'S 
FORMAT*:  A 1  > 

IF*:  IANS  .EG.  1H  )  GO  TO  555 
IF*:  IANS  .HE.  1 H  T  >  GO  TO  950 

WRITE*:  1,13) 

FORMAT*:  “NEW  GRIDS*  Y/N"  ) 

READ*:  i  ,  12  )  IANSS 
IF-<  IANSS  .EG.  1 H  )  GO  TO  560 
I FK  IANSS  .HE.  1HY>  Gu  TO  5*35 

CALL  POLAR*:  0  >  . 

CALL  LABL<  0  > 

WRITE*:  1  .  2S  ) 

FORMAT < "ENTER  PEN#  FOR  DATA.*") 

READ*:  1  ,  *  )  I  PND 
WRITE*:  25 , 22  )  I  PHD 
FORMAT*:  “SP“ , I  1  > 

DO  600  1=1, NT IMES+1 

PHI=>:<  1-1  >*DGfi£NC+tSO .  >=*F‘I/1S0. 

DBPRPT*:  I  )=20 . ‘♦'ALOGT*:  SCRPE2*:  I  >/TFLDMX  .) 

PHI 1=< 1-1  >**DGRENC 
WR I TE<  6 , 33  )  DBPRPT*:  I > , PH . I 
FORMAT< 1 X , “ GA I N  -  POWER  =  “ , E 1 5 . 8 , “  db “ , 5X , 
F?.3,X> 

I  f*:  dbprpt*:  i :»  .  l  t  .  -40.)  dbprpt*;  i  :>= -4  0 . 

DUMMY**:  DBPRPT*;  I  )+4  0  .  V40. 

I  X=  I  NT >;  *:  NCR  -DUMMY*COS<  PH  I  >*<*RDS  >**  1  0  0  0  .  ) 

•  IY=TNT *,  K  VCR- DUMMY*  SIN*:' F’HI  )•<• RDS >■*  1  000.  > 

I  FOAMS  .NE.  1 H  .AND.  IANS  .HE.  1 HY )  GO 
WRITE*:  25, 24  >  IX,  IV 
FORMAT*:  “  PA  ",  1 5  “  ,  "  1 5  "  ;  FD  "  > 

CONTINUE 

WRITE*:  25, 26) 

FORMAT*:  "PUU  ) 


“ANGLE 


TO  60  0 


ChLL  LaBEl*;  1  > 


935 


"54  0 

c 


CS5  0 

•545 

C 


C 

c 


c 

•54? 

c 

948 

C 

945 


C 

C 

951 


C 


00  94  0  1=1,  NT  I  MI'S. 

I-t  )-*DGREWC 

UR I TE<  6 , 935 >  DBBRPT < 1  >  ,PHI 

FORMaTC  IX,  “POWER  GAIN  *  “,£15.9,"  db 11 , 5X .  “ANGLE  =“, 
F? . 3 , “  DEGREES", /) 

.  CONTINUE 


UR ITE< 6,70  > 

WRITE*;  6 , 545  >  SFLL'ftX 

FORMAT*;  1 X,  "THE  ECHWPW-MAX  =  “  ,  E  1  5  .  S,  //  > 
WRITE*-'  b ,  957  >  T  FLO  NX 

FORMAT*;  IX,  "THE  TOTAL  F I  ELD-MAX®  “,£15. 85 
WRITE*:  6,70  > 


I AVNTM= 1 AVNTM/NT I  ME S 
UR  lit*;  6 , 94  ') 

FORMAT*;  IX,  "THE  AVERAGE  NUMBER  OF  TERMS  IN  thE  INFINITE  SUM 
WRITE*:  fc  ,  94  S  >  NT I ME S 

FORMAT',  1 X,  "FOR  “,  15,  “  NUMBER  OF  ANGLES  BETWEEN  0  ii-  360  IS" 
WRITE*.'  6 , 949  )  IAVNTM 
FORMAT*;  1 X ,  1 5 ,  "  TERMS .  "  > 

WRITE*;  6,70> 

UR I T£< 6 , 95 1 >  I MhX 

FORMAT*:  IX,  "AND,  THE  MAX  NUMBER  OF  TERMS  WAS  ",  15,  "  TERMS. 
CALL.  EXEC*;'  23 ,  5HSMB  ,4,1  SLU  > 


C 


END 

END* 


(')  '  j  O  «'J  T  If,  s&  (■-  <X>  O  O 


V  T=00004  I'-'  Or-!  CROC 039  US  INS  00005  BLK't-  R=  0 0 0 0 


FTN-! .  L 

*EMa<  CMBLK ..  0  > 

C  if  if if if if  if  if  if  if  if  *•  f  *  if  if  f  if.  if  *  (if  *  if if:  *  f  if  if  if  if  if  if  if 
SUBROUT I NE  CL  SKV< N , M > 


C  *  if  1*  if  if  if  if  if  if  :*  if  if  -f  if  if  if  if  f :  f :+  if  if  if  f :«  f  *  if  if  if  if  :*  if 

c 

COMPLEX  A<6,?>,F<6>, 

1  SUM 

C 

C  * 

COMMON  /CMBLK/  h,  F ,  ECHUPUK  30  0  0  >,  DGREmO ,  SCRPE2-  30  00  >, 
1  DBPRPT<3C0CO 

C 

DATA  P I /3 . 14J 59205/ 

C 

C  i*  CALCULATE  FIRST  kOw  OF  UPPER  UNIT  TRIANGULAR  MATRIX 
C 

DO  3  o.-2 ,  h 

A>;  1  ,  J  >-A<  1  ..  J  >/«■;  1  ,  1  > 

*  CALCULATE  OTHER  ELEMENTS  OF  U  AND  L  MATRICES 

DO  S  I=2,U 
0=1 

DO  5  II=J,N 
SUM=<  0  .  O',  0  .  0  > 

DO  4  K*>  t  , JM 1 
S'UK=SUrt+ A< 1 1 , K )  f Av  K, J > 

A< 11,0  >=A< II,  J >— SUM 
IP1=I+1 
DO  7  0  0= I P 1  ,  M 
SUM=<  0. 0, 0.  CO 
IM4  =  I  - 1 
DO  iS  K=  1  ,  IN  1 
SUM=SUM+A<  I  ,  K  >:kA<  K ,  JO  > 

A<  I ,  J  J  >=•:  A<  I ,  JO  >— SUM  >/A<  I ,  I  > 

CONTINUE 


if  SOLVE  FOR  F<  I  >  BY  BACK  SUBSTITUTION 
C 

F<  N >=A<  H,N+ 1  > 

L=N-  1 

DO  1  0  NN=  1  ,  L 
SUM=< 0.0,0. 0 > 

I=H-NH 

IPI=I+1 

DO  9  J= IP 1 , N 

9  S'UM«SUM+A<  I  ,  J  >fF';  J  > 

1  0  F'.'  1  >‘«A<  I ,  M  J-SUM 

RETURN 
END 
END* 


! 

J 


*.ruri 
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I  T-  mu 00 -4  IS  ON  OR 0 0 OS'S  USING'  00  000-  OL.k'J  *=0000 
FTH4  ,  L 

*EN«<  CMBLK ,  0  > 

C  •*  '■¥  :*  :*  *  ♦=  :f  -*  :*  ••»■  :*  *  ■+  f  :f  -*  *  •+  *  ¥  :*:* 

SUBRGUT  I  HZ  ECCP  T < XMAX  > 

C  +  f  +  ft  :■*+.  +-t  -  ♦  :«-*:*-*  :*.:*■ 

C 

•  COMPLEX  «■;  6 ,  ?  ') ,  F<  t-  > 

C 

C 

COMMON  /CMBLK/  A,  f ,  CCj-lWPU-:  3000  >,  DC  RENO,  '3CRPE2X  30  00  >, 
1  DBPRPT<  3000  > 

C 

C 

DaTA  P I /3 .14)  5S-2E.5.'' 

NX OP* 2 . * 1 000 
NY  0F'=  1  .  75  *  1  000 
C 

NXSTP* 0 . S* 1 0  0 0 
NYSTR=0 . 5*  1  0  00 
C 

NXSTPS*S 
NYSTPS* 1 0 
C 

I PEN*  t 

UPITE<25,3>  I  PEN 
C 

NXE  P*HX  OR +NXSTRS*  NXSTP 
WRITE’C  25,  1  >  NX  OP,  NY  OP 
UP  I  TE-;  25 ,  1  >  NXEP ,  NY  OP 
WRITE’,' 25. 2) 

C 

MY£P*NYOP+NYSTPS*NYSTP 
WRITE-',  25,1  >  NX  OR ,  NY  OP 
WRITE-;  25,  1  )  NX0P,NYEP 
WRITE-;'  25,2  > 

C 

DO  100  J* 1 , NXSTPS 
I X*NX 0R+ J*NXSTP 
IY1*NY0P-50 
IY2=NY0P+50 
WRITE-;' 25,  1  >  IX,  I V 1 
WRITE-; 25,  1  >  IX,IY2 
WRITE-:  25, 2  > 

100  CONTINUE 
C 

DO  200  1=1 , NYSTPS 
I  Y=MY  0F+  J  -uNYSTP 
I X 1  “NX  OP— 5-0 
1X2* NX OP +50 
WRITE-;  25,  1  J  I XI  ,  IY 
WEITE< 25, 1 j  1X2, I Y 
WRITE-;  25, 2  > 

200  CONTINUE 
C 

CALL  ECLBL 
CALL  LAPEL-:  1  > 

C 

HPTS*  IFIX-:  1  SO  .  XDGREHC  > 


3>  1PEN 


IPEH-=7. 
WRITE*:  20 , 


00  400  IPHI=1,MPTS 
I  X*  I F  I  ;>X  •;  I  PH  I  - 1  . 
I  V«  I  r  I  K  ECHWKu*-: 
UR  I  TEC  25,  t  IV 

CONTINUE 


>=«  OCkENC/2  0 .  NXS.TP+ 


I  PH  I  VXMAX 


>/u  .  1  =*MYS . . 


WRI  TE-,'  25, 2  > 

FGRriAT<  "PA"  ,  15“  ,  *•  15,  "  j*PD "  ;> 
FORMAT*  “FU“  .') 

FORMAT*:  » IN j  SP “  ,  I »  > 

END 

END* 


. 
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cn  o  ui  o  o  o  o  o  o 


COMPUTE  THE  BESSEL  FUNCTION  OF  ORDER  H,  H  hN  INTEGER  NXJR 
WITH  REAL  ARGUMENT 

ALL  E0UATI0N  REFERENCES  TO  ABRAM0WIT2  AND  S7EGUN 


D 1  MENS  I  ON  C0<  7  >,  C 1  <  7  > ,  D  0<  ?>,DK?),E  0<  ?>,Et  <  7  > ,  G  0  •  7  > ,  G 1  <  7  ) 
DATA  C  Cl/ 1  .  0,  -2 . 2499997,  1  .26502  00, 31630*6,  .  4444  79E-1 , 

’•*—  .  39444E~2 ,  .  2 1 E- 3/  . 

DATA  Cl / . 5, - . 562499S5, . 21 093573, 39542S9E-1 , .4433I9E-2, 
*-.3170  IE-3.  . 1 1 09E-4/ 

DATA  00/ , 79709456, - . 77E-6, - . 55274E-2, - . 9512E-4 , . l 37237E-2 
J*—  .  72S 05E—03,  .  1 4476E- 03/ 

DATA  D 1 / . 797S9456 , . 1 56E-5, . 1 659667E-1 , .171 05E-3 , - . 24*5 1 1 E 
•t  .  i  i  3653E  -  02 ,  -  .  2  0  033E—  03/ 

DATA  E0/- . 7053991 6, - . 4 1 66397E- 1 , - . 3954E-4 , . 262573E-2, - . 54 
.  29333E-  03,  ,  1  3550E—03/ 

DATA  E  i  /-2 . 356 1  944 9 ,  .12499612,  .  565E-4  ,  -  ,  637879E-2,  .  7434SE 
* . 79924 E- 03, - . 291 6-6E- 03/ 

DATA  CO/.  3674669,  .  6 055937/ - .  7435033,  .  253001 2',  -  .  42'6  1  2'  1  E-0 1 
* . 4£79i 6E-02, - . 24S46E— 03/ 

DAT  a  G 1 /- . 6366 1 96 ,  . 22 1 2 09 1,2.1 6S27 1 , - 1 . 3 1 64S3 ,  .31 2395 , 

-  .4  00976E-01  ,  .  27073E-  02/ 

DAT  A  PI /3 .141 5926/ 

IFLG-0 

IF  < M . L T . 0>  IFLG= 1 
H= I ABS<  N  > 

IF<  ABS< X  >  .  L T  .  1  .  0E-1  0  >G0  TO  1  50 
IF.;  AE:S<  X  •> .  GT  .  3 . 0  :>GO  TO  5 0 
X3SG=X*X/9 . 0 
PROD* 1 . 0 
B 0=  0 . 0 
61  =  0,0 
CUM 0=0 . 0 
CUM  1  =  0.0 

SEE  EQUATIONS  9.41  AND  9.44 

DO  5  1=1,7 

B0=B0+C0<  I  >*PROD 

BI=B1 +C1< I  XfPROD 

CUM  0=CUM  0+G  0>:  I  ><*PRGD 

CUM  1  =CUri  1  +G 1  >;  I  >*PRuD 

PR0D=PR0D*X3S0 

CONTINUE 

B1  =B  1  ;*X 

XC=2 . 0 *SNGL<  DLOG<  DBL  E<  0 . 5*X  >  >  >/P I 
Y  0=XC*B0+  CUM  0 
Y I =XC*B 1 +CUr1 1 /X 
GO  TO  100 

EOS  9,4.3  AND  9.4.6 
0  THR0VX=3. 0/X 
RR0D-- 1.0 
F  0=  0 . 0 
F  1  =  0 . 0 
THETA 0=X 
THETA  I =X 
DO  55  1=1,7 
F  0=F  u+DO-C  I  >*PftuD 
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irjaua 


F1=Fi  +DK  I  >=<PROD 
1 HET AO=THETaO+EO\  I  J-t'PP.OD 
THcTrtl-THEJAl+EI*.  1  ):i*PKL:E> 

PROO=FROC>;*TriROVX 
55  CONTI HUE 

SGRX®  1  .  O/SNGl.1;  DSORK  DBLEC  X  >  >  > 

B0=SGRX*F  0*SNL-L<  DCOSC  DBLE*.  THETwu  >>> 

&  t  =SC'R.X*F  1  *SnG'L<  DC  OS  <  &£l£<  THET  ft  1  >>> 

.YO*SC'RX>*F  0*SNGL<  DS1N<  DBLE<  THETA 0  >  >  > 

V 1  =SORX*F  I  ■*SNOL<  DSIW  DBLE1'.  THETA )  >>> 

100  I F<  H- 1  > 1 0 1 , 1 05 ,110 

lOt  B=B0  . 

BP=-B i 
Y=Y0 
YP=-Y1 
GO  TO  2 0  0 
1 05  B=B 1 

BP=B 0-& 1 /X 
Y=Y  t 

YF-YG-Y1/X 
GO  TO  2 0  0 

C  FOR  RECURRSIVE  DIRECTION  COMMENTS  SEE  SECTION  9 .12, P3 

1 1 0  XN=-M 

IF*:XN.LT.ABSO;>>GO  TO  130 
C  FOR  X<  N  RECUR  DOWNWARD 

BLAST  ® 1  .  0 
BLASTF-  0 . 0 
J-H+ 1 0 

DO  115  1=1, J 
XI=U-I 

BNEXT =2 .  0*XI*BLAST/X.-BLASTR 
ELASTP^&LAST 
BLAST  -BNEXT 
IF< I .HE. 1  COGO  TO  115 
BLF-BLASTP 
Z-BNjjXT 
115  CONTINUE 

IF<  ABS< B  0  > , LT . ABSCB 1  ) >GO  TO  11? 

CORR=B  0/BLASTP 
GO  TO  1 1 S 

11?  CORR=-&1 /BNEXT 
1  IS  B=BLP«CORR 

BHMIN1=2*C0RR 
GO  TO  140 

C  FOR  X  >H  RECUR  UPWARD 

130  &LASTP=B0 
&LAST=B1 
DO  135  1=2, N 
XI® 1-1 

BHEXT=2 . 0+XI *BL AST /X -BLAST P 
BLASTP® BLAST 
BLAST-BNEXT 
135  CONTINUE 
&= BLAST 
BNMIH I =BLASTP 
140  BP=BNM I N 1 -XN*B/X 
VLASTP=V  0 
YL  AST  »  Y  t 
DO  145  1=2, N 
XI® 1-1 
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1 1  iC.Ml  :  +  V  L  1-i  S  i  —  Yl  .  mS  T  F 

VlAST  p=YLttS'r 
i'l  AST ~  v  NEXT 
MS  CONTINUE 
V=YLAST 

VP  ~  VL  A  ST  P  -  XN  t  Y/X 
GO  TO  SOti 
1 5  0  8=  0 . 0 

BP=  Ci .  0 
YsO.  0 
VP  -  Ci .  c< 

IP';  N-l  >  I  55 ,  ISO,  2  0  Ci 
tSS  P= 1 . 0 

GO  TO  2 0 0 
160  SP-0.5 

e* 0 . 5 ♦ X 

200  IF  <IKLG.EQ.CO  RETURN 
8=  <  - 1  )*-*iC*B 

N*-N 

RETURN 

ENC> 

END* 


oooo 


ip  T=00004  IS  ON  CP0003S  USING  00005  ELKS  R=  00 00 


c 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


I  .  r;A  I N  Prc-gran. 

I  Reads  i  nputs 

2>  Calls  CELL 
a  )  AN 
t.  >  NCELLS 

3)  Calls  CLCftD 
a  >  XN 
tO  YN 

4  >  Calls  ECHOUi 

a  >  UK  PH  I  .WAVELENGTH 

II.  ECHOUi  *;  need  Ei  and  En ) 

1  >  Cal  Is  FDLTL 
ci  5  En 


C 

c. 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

C 

c: 

c 

c 

c 

c 

c 


2>  Calls  FLDNC 

a  0  Ei  as  a  function  of  PHI 

III.  FLDTL  tneed  En,i  > 

11  Calls  FLDNC  with  PTOBS  =  0 
al  Emi 

21  Sets  up  niatrix.  and  calls  CL  SKY 
a>  En 

»>  Have  Er.  and  Ei  >>>  Have  ECHWPW 
IV.  MjIN 

5  >  Calls  FLOCT 

a>  Es< rauo, phi >  <<<  Knowing  En  which 
EMA  calculated  in  II, I 


is 


i  n  an  srra 
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0004)00  G  G  O  G  O  G  O  G  GO  GO  GOG  G  G  G  G 


W  T  =  00  Cum  IS  Oh  CR00036  USING  6 0020  BLKS  R-uOC'O 
FTii-4 ,  L 

*EMA<  BLKMn,  0  > 

C  't:l*  r  tf  -ft  t"  f  :rt< 

PROGRAM  PADriM 

c 


ROCfcRT  k.  'SCHNEIDER 


This  program  uses  the  MOMENT  METHOD  a?  presented  fcv  JACK  RICHMOND 
< “Scattering  by  a  Dielectric  Cylinder  of  Arbitrary  Cross  Section 
Shape",  IEEE  APS,  Nov  t ?e4 ,  p .  334  >  to  determine  the  Scat  tered 
Field  from  and  Echo  Width  of  a  dielectric  cylindrical  shell  of 
c  i  rcu  1  ar-  cr  oss  sect-  i  on  . 

*  The  dimensions  of  the  scatterer  are  read  from,  a  data  file  -  DAThFL 

INTEGER  NCELLS, NFTS, PTOBS, IPHI,HCEtLS 

REAL  XNN , VNN , EMAG I , Ji , J 0 , k 0 , K2 , V 0 , Y 1 , EMACN , 

1  ECHWPW 

COMPLEX  CMN , E I NC I , E I NC 0 , E I NC2 , 

1  ESCAT , 

2  TALM  ,ALMDA,  ALPHA,  COEF,  FAC,  Cl  ,C2,V,V  1  ,V2, 

3  T  AU ,  V I N ,  VOUT  ,  Ot IE ,  ZERO ,  TLF LD ,  A  0 ,  A 1 

DIMENSION  IBUFR< 16 j 

COMMON  /BLKMM/  AA, B, C, R, FREO, PERM, DCRENC, EPSLNR, AN, NCELLS, 

1  El NC 1 < 36 0 j , V I N< 1 S5 f  > , ESCAT *1 36 0  > , 

2  ECHWPLK  36  0 > , K 0 , K2 , VOU'i  f i S5 1  > , TALK  1651  > , 

3  *  E/HC2,  SCRPfc  2<  361  > ,  DBPRPT  >,'361  >,  A  fit  1651  > ,  A 1 *.'  1 SS 1  > 

EGU I VALENCE  < NCELLS , HZ  > 

DATA  P I r'3 .  14  1  55265/ 

DaTA  XCR/4 . S?5/, VCR/3 . 675/ , RDS/2 . 746/ 

DATA  IBUFE/2* 0 , 2HD A , 2HTh , 2HFL ,3*0, 22 1 B , 7* 0/ 

*  SET  SPOOL  FOR  INPUT  FROM  IBUFR  THROUGH  ISLU 

CALL  SPGPN< IBUFR, 1SLU> 

CALL  EXECf 22 , 1 > 

READ*;  I SLU,  *>  AA,B,C,  CURE  NT  ,  R ,  FREO ,  PERM ,  DGRENC,  EPSlNR  ,  NCELLS , 

1  HLFCEL 

*  DETERMINE  THE  CELL  STRUCTURE  AND  #  OF  CELLS 


C**** 

c 

C»*w 


W=  0 . 0 

CALL  CELL*.'  if  .> 

I Ft HLFCEL  .NE.  1>  GO  TO  50 
NCELLS*NCELLS/2 
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DETERMINE  THE  COORDINATES  Or  Tut  CENTER'S  Oc  £<»Ch  CElu 
CALL  CLCRD\W> 

DETERMINE  THE  ECHO  WIDTH  RER  WAVELENGTH 
XMAX-  Ci . 

CALL  ECHOW-;  XM AX ,  W  > 

NPTS=  I F I :  X  36 0  .  XDGRENC  > 

WRITE-,  1  .  23> 

FORMAT-;  "PLOT  THE  "SCATTERED  FIElD''*  Y/N”  > 

READ-;  1  ,  1 2  >  I  AMS 
FORMAT-;  A  t  > 

IF< IAnS  .ED.  1H  >  GO  TO  210 
IFCIANS  .NE.  1HV>  GO  TO  220 

PLOT  ECHO  WIDTH  PER  WhVELENGHT 

CALL  ECOPT < XMAX > 

DETERMINE  THE  FAR  ZONE  SCATTERED  FIELD  AT  EaCH  DCRENC 
CALL  FLDCT 

DETERMINE  THE  INCIDENT  FIELD  IN  THE  FaR  ZONE 

CALL  FLDNCX  1  ,  W  > 

TFLDHX= 0 . u 

DO  260  1=1,  NPTS 

TLFLD=EIHC1-:  I  >+ESCATX  I  ) 

SCRPE2-;  I  >*CAB3<  TLFLD  > 

4F<  SCRPE2-;  I>  .  GT  .  TFLDMX  >  TFLDr1X=SCRFEZ-:  I  > 
CONTINUE 

SCRPE2C  NFTS+ 1  !>=SCRP£2*:  1  ) 

WRITE-;  1  ,24  > 

FORMAT-,' “PLOT  THE  TOTAL  FIELD?*  Y/N"  > 

READ-;  1  ,  1 2  !>  IANS 
IF<  IANS  .£Cl.  1 H  >  GO  TO  230 
IF-,' IANS'  .NE.  1 HY >  GO  TO  270 

WRITE-;  1 ,25> 

FORMAT-,'  "NEW  GRID?*  Y/N"  > 

READ-;  1,1 2  >  IANS 
IFCIANS  .EQ.  1 H  >  GO  TO  240 
IF<  IANS'  ,  HE  .  *HY>  GO  TO  250 

CALL  POLAR-;  0  ) 

CALI.  LAPLC  0  :> 

WRITE< 1,23) 

FORMAT-;  “ENTER  PEN#  FOR  DATA.*  “  > 

READ-:  1  ,  =*  ">  I  PHD 
WRITE-:  25, 22  >  I  PND 


22  FORMATS  “8P" ,  I  i  > 

C 

Du  700  1=1, NPTS+ 1 

PHI**C<  1-1  >■«  DGREr-)C+ 1  80 .  >*PI/i8u. 

DBPRPT ■:  I  >=20 .  ;*:ALOGT<:  •3CkP£2*:  I  VTFlDMX  .> 

IF' DBF  RAT*;  I  >  .LT.  -4  0.  >  DBPRPT'.  I  >=-4  0. 

DUMMY-1  DBPRPT*:  I  >+4  6 .  >.'4  0  , 

IX*  IMT<  <  XCA -DUMMY-*  CO'S<  PHI  >:*RDS  >*  :  000.  > 
l  Y=  ]  NT 0  0  YCk-DUMMY--*o I Nr,  PHI  >:*RDS  >«•  j  000  .  > 

WR  ITS <  25 , 29  >  I X ,  I V 
28  FORMAT*:'  “PA",  15“,  “I5~;PD‘*  > 

700  CONTINUE 

WRITE1;  25, 3  1  > 

31  FORMAT1.' “PU"  > 

C 

27  0  WRITE*;  1 , 28  > 

26  Fui<MAT<  "HARD  COPY  OF  OALCULAT  I  ON'S  ?  5  Y/M”  > 

READ*:  1 ,  i  2  >  I  AMO 
IF1.  I  AM?  .EC'.  1H  >  GO  TO  800 
IF1,' IANS  .ME.  1HY>  GO  TO  500 
C 

S  0  0  WR  I  TE>:  6 , 85  0  >  T FL  Dn>. 

85 0  FORMAT1;  IX,  "TOTAL  FIELD-MAX  =  “,E15.3,//> 

C 

DO  5:70  1=1,  NET'S 
F'H  I  =<  1-1  >'*DGRENC 

WRITE1;'  6,  88 0  >  PHI ,  ECHWPuK  I  >,  DBPRPT I  > 

88  0  FORMAT < 1 X , " «NGLE=  “ , F5 . 2 , 5X , " ECHWPW*  " , E 1 5 . 8 , 5X , 

1  "DB  GA I N*  " , E 1 5 . 8 , X > 

870  CONTINUE 

C 

WRITE< 8, 27  > 

27 '  FORMAT< 1 X, XX/ > 

C 

DO  400  1=1, NPTS 
*  PHI=< 1-1 >»DGRENC 

XMAG=CABS'.  t  ESCATC  I  >  >  > 

WR I TE<  6 , 3  0 0 >  XM AG , PH I 

300  FORMAT*:  IX,  "MAGNITUDE  ESCAT="  ,2X,  1E15.S,5X,  "PHI*"  , 

400  CONTINUE 

C 

500 


s 


X , F7 . 3 , X  > 


CALL  EXEC*.  23 , 5H8MP 
END 
END* 


4, ISLU> 


o  o  o  o  o  o  o  o  o  o  t>o  oooooooooooo  o  o  o 


R=  [■  C<  0  0 


1  =  00004  IS  ON  CROC 03’?  USING  00004  ELKS 
F1N4,L 

*Er»A<  BLKMM ,  0  ') 

c  if  if  if  if  if  if  if  if  f  if  if  if  if  if  if if  :f  if  •*  t  if  :f  f  f  ;f  if  -f  if  if  if  if  if  -f  if  if  if  if 

SUBROUT  I NE  C  ELL-;  W  .> 

C  :*!*  !*!<•!♦  ■*  *  *  if :« ■■*■*+*+:*  :<■■**  ■*  *  ■*  **■+  *  ****  •*  *■■*#+  <*  ■* 

:*  THIS  SUBROUTINE  DETERMINES  THE  #  OE  CELLS  NEEDED  PHD  THE 
DIMENSION  OF  A  CELL  TO  SFttN  THE  OBSTACLE. 


INTEGER  NO ELLS, NPTS, F TOES, I PHI , M CELLS 

REAL  XNN, YNW,EMAGI ,  01  ,  JO,  KO,  K2,  VO,  Yf  , EnhGN, 

1  ECHWPU 

COMPLEX  CNN , E I MC 1 , E I NC 0 , E I NC2 , 

1  ESCAT , 

2  TAU  1  ..  ALMDA ,  ALPHh  ,  OOEF  ,  FAC ,  C 1  ,  C2‘ .  '■> ,  V 1  V2 , 

3  Thu, VIN, VOUT, ONE, EERO, HO, HI 

D I  MENS  I  ON  I  SUFR'i  1  B  > 

COMMON  SBLKM1S  HH,  B,  C,  ft ,  FRED  ,  PERM  ,  OGREMO .  EPS’-K.R  ,  hM  ,  NCELLS  , 

1  E I  NC  1  •;  3c- 0  >  ,  V I N<  1  S5 1  ),  ESCHT-;  360  >, 

2  EChwPw';  360  >,  KO,K2,VOUT<  i  ?5i  ;,THUt  IBSi  >, 

3  £  I NC2 ,  SCRPE2<  36  i  >,D&PRPT\  36  i  > ,  A  O'.  iO?i  >,AU  1S51 

EQUIVALENCE  < NCELLS, N2 > 


DATA  PI/3. 14150265/ 

*  THE  #  OF  CELLS  FOR  THE  CIRCULAR  DIELECTRIC  SHELL  SCATTERER  IS  A 

parameter  which  is  assumed  at  the  initiation  of  the  program  and 

IS  READ  FROM  THE  DATA  FILE 

THE  WIDTH  OF  EACH  CELL  MUST  BE  LESS  THAN  OR  EQUAL  TO 
0.2  *  WAVELENGTH  /  SORT  <  EPSLNft  > 

■«:  SINCE  THE  NUMBER  OF  CELLS  IS  A  KNOWN  VALUE  AND  THE  SICE  OF  THE 
STRUCTURE  IS  DEFINED,  THE  WIDTH  OF  A  'SQUARE  CELL",  W,  IS  EASILY 
DETERMINED  TO  BE 


W=2 . I *C/NCELLS 

*  RADIUS  OF  CIRCULAR  CELL  WITH  EQUAL  AREA  AS  SQUARE  CELL',  AN,  IS 


C: 

C2  00 
1  00 

1 

C 


ROUT-C 

RIN=C-W 

AN=SQRT<  <  ROUT 2 .  .  > /NCELLS  > 


WRITE< 6 , 1 00>  W, AN, NCELLS 

FORMAT  K IX, "W= “ , 2X, 1E15.S,5X, "AN= “ , 2X , 1E15,S,5X, "NCELLS* 
2X,I4,//> 


RETURN 

END 

END* 
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oo  o  o  o  ui  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o 


If*** 


D  T=00u04  IS  ON  CR00O3S  USING  60004  BLni  F:=0000 
FTN4 , L 

*EMA<  BLKMM .  0  > 

C  '*•*'*  I*:****  *c*<  H  *:«•*  . f  t.  f  *:*:*  *■**»  (*  (.:*:*■**  +  *+* 

SUBROUTINE  CLCftlKU) 

c 

*  THIS  SUBROUTINE  CALCULATES  THE  COORDINATES  OF  THE  CENTER  OF 
EACH  CELL . 


INTEGER  NCELLS, NPTS , F TOBS, I PHI, NCELLS 

REAL  XNN,  YMN,  EMaGI  ,  J1  ,  J 0 ,  K  0 ,  K2 ,  Y 0 ,  V 1  ,EnAGN, 
t  ECHkiPW 

COMPLEX  CNN .,  E INC  1 ,  EINCO,  EINC2, 

1  E SCAT , 

2  THU  i  ,  ALMDA ,  ALPHA ,  COEF ,  FhC  ,  C 1  ,  C2 ,  V ,  V 1  ..  VO  . 

3  '  T  AU , V I N , VOUT , ONE , ZERO ,  A  0 ,  A 1 

DIMENSION  IBUFRC 16 > 

COMMON  /BLkNM/  AA ,  £< ,  C ,  R ,  FREO ,  PERM ,,  DGF.ENC,  EFSLNR ,  AN, NCELlS, 

1  E I NO  1  <  36 0  >,  VI N<  37  02  > ,  ESCAT *,'  36 0  > ,  ,v  ;Nh<  37 02  > , 

2  YNNC  37  02  > ,  ECMUPlk  36  0  > ,  h  0 ,  k'2  .  VCn.iT <  37  02  > ,  T 3  T  C 

3  El  NCR ,  SCRPEZ*;  36  i  > ,  DBPPPT  6  361  >,h  0<  37  02’>,A<<37  02 

EQUIVALENCE  < NCELLS, N2 > 


DATA  PI/3. 1 4 1 59265/ 

|*  RADIUS*  OUT  TO  CENTER  OF  CELL  r. 

RN=C-CW/2. > 

*  INCREMENT  ANGLE  FOR  CELL  LOCATIONS 
THETAH-2 . *P I /NCELLS 
WRITE6 6, 50 )  THETAN 

0  FORMAT  6 1 X , “ THE  T AN= “ , 2X , 1 E 1 5 . 9 , // > 

=*  DETERMINE  THE  CORD  I HATES  OF  THE  CENTER  OF  EACH  CELL  n 


1  00 

1  0 
c 


DO  10  N= 1, NCELLS 

theta=*:n-i  >:*thetan 

XNN<  N  >«PN*COS<  THETA  > 
YNNC N  >*RN*SI N*;  THETA  > 


UR  I  TE<  6 ,  1  0  CO  XNN'.  H  > ,  YNN>:  N  > ,  N 
FORMAT'.'  1 X ,  " XNH=  *  ,  2X ,  t  E 1 5  .  * ,  5X ,  “  YNN=  "  ,  2X ,  1  £  1  5  .  ?  .  5N ,  “ N  - 
CONTINUE 


zx , 


14 


RETURN 

END 

END* 
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o  o  o  o  m  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  r*  o  o  o  o 


W  T-  u 0 (i 0**  IS  ON  Oft 00 OSS  USING  UUCU'S  CLMS  P=*0u00 
FTN-4 ,  L 

*EMA*;BLKMM,  0> 

C  <*-«**4‘>*>*‘**»*<***+-t.--*+--**--»*+-4*  «***«>**i*i» 

SUSROUT I  ME  E CHOW-;  XMaX  ,  W  > 

*  THIS  SUBROUTINE  CALCULATES  THE  ECHO  WIDTH  PER  UNIT  i>  avElEMGHT  , 
RELATIVE  TO  THfc  INCIDENT  FIELD  aT  f : ■  E  CEO T L •-  OF  THE  OESFhClE, 
FROM  h  DIELECTRIC  CYLINDRICAL  SHELL  CF  CIRCULAR  CROSS  SECTION 
IN  THE  PRESENCE  OF  a  PhD I AT  INC  CURRM  r  FI LAMENT. 


INTEGER  NCElLS,  NETS,  PTOBS,  I  PH  I ,  NCELLS 

REAL  XNN , VNN , E MAO I ,  U1  , J 0 , K 0 , k2 . V 0 . VI , EMAGN , 

1  ECHWPW 

COMPLEX  CNN , E I MC 1 ,EIHC0,EI NC2 . 

1  ESCAT , 

S'  T ftU  1  ,  ALMDA ,  ALPHA ,  C  OEF ,  r  wC  •  C 1  .CL,  V ,  V 1  ,  V2 , 

S  T  AU , V I N , VOUT , ONE , ZERO , DUMMY4 , A  0 , « 1 

DIMENSION  IEUFRC IS  i 

COMMON  SBLKNMS  AA, B, C, R, FRE 0, PERM, OCRENC, EF SlnR, AN, NCElLS 

1  E I NC 1  *'  36  0  > ,  V I  N<  ISOt  ) ,  ESCaT \  36  0  > , 

2  ECHWPUC  360  >, K 0,  !C2,  VOUTC  1051  >,ThU<  1SS1  ), 

3  E I NC2 , SCRPEZC  36 1  > , DBPRPTC 36 1  ) , A 0<  1  S?  1  >,A1 < 1 

EQUIVALENCE  < NCELLS, N2 > 


DAT  A  P I /3 . 1 4 1 59265/ 

* 

*  CALL  FLDTL  TO  DETERMINE  THE  TOTAL  FIELD  IN  THE  OBSTACLE 

DUE  TO  THE  INCIDENT  FIELD  UPON  THE  OBSTACLE 

CALL  FLDTL<  W > 

*  CALL  FLDNC  TO  DETERMINE  THE  INCIDENT  FIELD  AT  THE  CENTER  OF 

THE  OBSTACLE, 

CALL  FLDNC-;  2,  IJ> 

NPTS= I F I X< 36  0 . /DGRENC > 

DO  7  1=1 , MPTS 

WRITE*;  6, 5  >  EINC2 

FQRMATC 1 X , "ECHOW  :  E I NC2= " , 2X , 1 E 1 5 . S , 2X , E 1 5 . S , / > 

7  CONTINUE 

*  DETERMINE  THE  MAGNITUDE  SQUARED  OF  THE  INCIDENT  FIELD 

EMAG I *CABS<  <  E I NC2  > >*CABS<  <EI NC2  >  > 

C 


C 


DUMMY  1=K0*«P  1**2  .*FREG/<3  ,ES*EMAGI  > 


r.  •* 

c 

C 


LOOP  NOELLS  Ti MEG  F Of:  THt  “SCATTERED  FIELD1'  I >- 

DO  20  I  PH  I  - 1  ,  NP1S 

PHI=<  IPHI-I  .  DCRENOPI,  1 1  SO  . 


h  rty.  2\jNc.  . 


£>UnriY4=<  0  , ,  0  > 

C 

RN=C->:  U/2 .  > 

THETAN=2  .  *PI/NC£lL'3 

c 

DO  10  n=1  ,  NCELi.S 

THETA=.;N-1  >*THETAN 
KK’N’=PN  *COS<  T  NET  A  ) 

YNN=RN*$ I N< THETA  > 

ARC J t -K  0:*AN 

C«LL  E:E6'EL*;  I  /  R.Pb  vi  I  ,  E;B  ,  Y  ,  £:P ,  YP  ) 

J  1  =BB 

hRGH=K 0*'.  XNN^CCG-.'  PHI  >+YNN*$IN\  PHI  >  > 

DUi-iM'/  2~C  uC<  ARCH  > 

DUMMY'S- 6  I  N<  hRGH  * 

BUWMY4 =DUMr Y4 *■<  EPSLNR-1  .  >%VOUT<  N  Hkn*  di  rtP>..'« DUMMY 2,  Dummy  ?  > 
10  CONTINUE 

C 


c 

cc 

cc 

c 

c 

c 

1  00 

c 

20 

c 


ErlACN*CABS<  DUMMY4  >*CAB£<  BUMMY4  > 

ECHWPfcK  I  PH  I  >*DUMMY1*EMAGN 
IF<ECHUPW< IPHI >  .  CT.  XriAX  >  XMAX-ECHWPU*;  I PH!  ) 


PHII*PHI*130./PI 


UR I TE< 6,100)  ECHWPUK  I PH I  ) , PH 1 1 
FORMAT.;  IK,  “ECHUPW*“,2X,  1E15.S,5X,  "PHI=“ ,2X, F7 . 3,/ > 

CONTINUE 

RETURN 

END 

END# 


i 


t 
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I 


1 


i  f  —  0  L"  0 JO  LN  LK  l.*0  l-v’;'  (J--.  I  No  Oi'OOV  *  *'  ■  P'jj  i_.  I.i 


F  Tii-4 ,  i_ 

•fcfW',  &LKMN ,  &  > 

C  -*J-  »  C  M-1(  f'1  JJ  J  f  t  +:.f  .J  •* 

SUBROUTINE  FLDTL-;  U  > 


t  f  Jt  f  ft*:*  :f  f  f  J  J  +  J  :f-  f  *  +■  :fc  *< 


■*  THIS  SUBROUTINE  CALCULATES  THE  TuTHL  FIElD  In  THE 


INTEGER  NCELLS, i-tRTS ,  PTOBS . IPHI, NCELLS 


REAL  XNN,YNN,EMa2I ,  Jl  ,  JO ,  s 0 ,  K2 ;  VC  ,  V 1  , EMaSN, 
ECHUPW 


COMPLEX  CrIM ,  S I MO l  ,EI NC 0 ,  E I NC2 , 

ESChT , 

TAU  1  ,  ALMDA  .  riL  RHA ,  COEF  ,  FAC  ,  C 1  ,  CP ,  V  .  V 1  ,  •» 
T  A  Li ,  V I  f ; ,  VOUT ,  O'  i£  .  EERO ,  A  0 ,  A 1 


L1 1  PlicNS  I  ON  I  fc.’Ur  H\  1  C*  .> 


COMMON  XK.KMM/’  AA,B,C,R,  FREQ,  PERM..  DGPENC,  EPSi 
E I  NC  i  ■;  BE  0  > ,.  V I  Nn  1  SB  l  >  .  BSCA  7 <  76  0 
ECHWPLK  76>  (i  > ..  K  C ..  K2 ..  VOUT ■;  1  f:  !  > ,  T» 
E I  MCE* ,  SCPPEZ-.  3c.  l  .> ..  DBPRRTs  3o  i  > ,  .*■ 


LHP, AN, NCELLS, 


!  :  E  !  ), 
1  f't'  I  >  , 


Lin  i’t. )  .■> 


EQUIVALENCE  < NCELLS, N2 > 


DhTA  PI/3. 1415S265/ 


C  *  OBTAIN  THE  FIELD  INCIDENT  UPON  THE  OBSTACLE 


GAEL  FLBNC<  0,  Li  > 


DO  5  1=1, NCELLS 

UR  I  TE<  6 ,  ?  >  £  I  NC  0<  I  > ,  I 

.  FORMAT-;'  IX,  "FlDTL  :  EINCO*"  ,  2X,  t  £  1  5  .  S,  2X,  E 1  5  . 

1 4  ,  /  > 

CONTINUE 


.  5X, " I  =  * 


RN=C-<  U/2 .  > 

THETAN =2 , *RI/NCELLS 


t>0  30  M=1,  NCELLS 


M=  1 

THETAH=-:  M- 1  >■* THETAN 
XNM=RN-*COS*:  THETAM  > 
YNp!=RNm'SIN<  THETAN  > 


DO  20  N=1,NCtLLS 


C  ■*  DETERMINt  THE  L'EIGHl  Cmn  OH  En 
C 

theta=>:n-i  >ithetan 
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AD-A115  517  AIR  FORCE  INST  OF  TECH  HRIOHT-FATTERSON  AF8  OH  SCHOO—ETC  F/«  17/9 
THE  EFFECT  OF  RAOOMC  SCATTERING  ON  ECM  ANTENNA  PATTERNS. (U> 

DEC  Si  R  K  SCHNEIDER 

UNCLASSIFIED  AFIT/0E/EE/S1D-52  NL 


XNN-RN*COSY  TWSTA > 

YNN=PN * S I N-; THETA  > 

IF-; N  . HE .  1  >  GO  TO  1 0 

DUMMY 1 =K 0*  AH 

CaLL  BESEi.-:  I  ,  DUMMY  t  ,  BB ,  Y ,  BP ,  VP  > 

J 1  =BS 
Y 1  =Y 

WRITE-:  G,S>  Jt  ,  Y 1  ,  BB 

FORMAT-:  IX,  “FlDTL  s  Jl  =  "  ,  2X,  IE  «  5  .  8,  5. UY1*",2«, 1E13.8, 
“BB= “ , 2X, I E 1 5 . &,  / ) 

CMH=  1  .  +<  EPSLNR-  1  .  >*<  <  P I  -*K  0-*-AN/2  .  >*C  Cnf  :_X>.  Y  i  ,.j1  >>+1  .  > 

A*:  M  ,  M  >=Cr.N 

if-:m  .HE.  i>  GO  TO  £0 
TALK  M  >=CMN 
GO  TO  CO 

DUMMY  1 =P 1 *K  0  *  AM/C . 

rmn=sqrt>:  •:  xnm-xnn  >*<  xnm-xnn  >  + 

-:  VNM-YHN  >*<  YNrt-VNN  >  i 
ARGH=K  0*RMN 
ARG  J=K 0-aN 

CALL  BESEL-:  1  ,  ARGJ,  BB,  Y,  BP,  Y P  > 

Jt=BB 

CALL  BESEL-:  0 ,  ARCH,  BB ,  Y ,  BP ,  YP  > 

JO-BB 
V  0= Y 


CMN*DUMMV  1  «<  EPSLNR- 1  .  >*  J1  *CMPLX<  YO,  J0> 

C  I F< M  .ME .  1  >  GO  TO  Cu 

TALK  N >=CMN 
C 

C  A-:  M,  N  >=CMN 

20  CONTINUE 

C'  ACM/W-H  >=-£JNC:0*;'M  > 

C30  CONTINUE 

C 

C  *  NOW  THAT  THE  MATRIX  HAS  BEEN  FORKED,  SOLVE  FOR  En 
C 

MCELLS=NCELLS+ t 

C  WRITE: 0,25 CO  MCELLS, HCELLS 

25  0  FORMAT-:  IX,"  MCELLS*  “  ,  2X ,  14, 5X,  ”  NCELL  S=  ”  ,  2X  ,14,  //  > 

C  DO  400  M- 1 , HCELLS 

C 

C  DO  300  N=1, MCELLS 

C  WRITE<S,350:>  A<  M ,  H  > ,  M ,  H 

350  FORMAT*:  IX,  “FLDTL  :  A=  ”  ,  2X,  1  El  5 .  S,  2X, El  5 . 9 , 5X,  "M=  ”  ,  2X ,  14  , 
1  5X, ’N=“,2X, 14, /> 

C30u  CONTINUE 

C  WRITE-:  6, 375  i 

375  FORMAT-:  IX,/// > 

C400  CONTINUE 


DO  BOO  1=1 , NCELLS 

WRITE< 6, 700 >  TALK  I  >,  I 

FORMAT-.'  1 X ,  “TAU=  “  ,  2X ,  1 E 1  5 .  S> ,  2X ,  E 1  5 . 9 , 5X ,  1 4 ,  /  > 
CONTINUE 
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Wi  O  O  O  O  O  (V 


Cf-iUL  ffti-iT'.  1  #  XNURil #  lER!* 

WRITE*:  1  ..500  >  I  Eft 
00  FORMAT* IX, “  IEft**“ ,  2X, 14  > 

Er.  ARE  CONTAINED  IN  THE  ARRAY  VGUT  OF  DIMENSION  NCEU.S 

DO  1  00  I=1,NCELLS 
WR I  TS<  6 , 2  0 0  >  VGUT >:  I  > ,  I 

00  FORMAT*:  IX,  "FLDTL  :  VGDT=="  ,  2X ,  1 E 1  5 . 9 , 2X.  El  5 . 9 , 5X,  "  I«  "  , 

1  2X,I4,X> 

Cl  00  CONTINUE 

c 

RETURN 

END 

END* 


Vi 


8  i 

'  I 


•'  1 

i  '  ’ 

r '  i 


i 
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i  n  r-v — 


« 


ooooo  O  OOO  OOO  O  O  OOf 


:  T-i'CuuM  IS  ON  CROC*  039  USING  00009  Si.KS  R=000u 
FTH4,L 

XtMf'N'  SI  KMM,  0  > 

C  *■  V  t  f  :**:**:*:*.+  f  If  **  *:* 

SUBROUTINE  FLDNC<  PTOBS,  W  > 

C  *****  *  ♦  f  *  *  *  *  .*  *  **  * :+  :**  *  *  *  *  *  *  * ;»  * 

c 

C  *  THIS  SUBROUTINE  CALCULATES  THE  INCIDENT  FIFLD  ON  THE  OBSTACLE, 

C  AT  THE  FAR. ZONE  POINT,  AND/OR  AT  THE  CENTER  OF  THE  OBSTACLE . 

C 

C 

IN i EGER  NCELLS, NPTS, PTOBS, IPHI , MCtLLS 
C 

REAL  XNN ,  YNN , Ei'iAG  1,01,  JO,  KQ,  K2 ,  Y 0 ,  Y 1  ,  EMA GN , 

1  ECHWF'W 

C 

COMPLEX  C'MH,EINC1  , E I NC0, E INC2, 

1  ESCAT , 

2  T AU !  ,  At  MDA ,  At  PHA ,  CC'EF ,  F  AC ,  C 1  ,  C2 ,  V ,  V  l  ,  V2 , 

3  ‘  TAU,  VIN.VCUT,  ONE, ZERO,  A0,  A1 
C 

D I  MENS  I  ON  I BUFRL  1  6  > 

C 

COMMON*  /BLKMM/  AA ,  8 ,  C ,  R ,  FREQ ,  PERM ,  DGRENC ,  EPSLnR  ,  AN ,  NCELLS , 

1  •  EIHCi  *C  360  >,  VIH*;  1851  >,  ESCAT<  360  >, 

2  ECHWPLK  38 C<  > ,  K 0 ,  K2 ,  VOLfT <  1851  >,TAU<  1851  >, 

3  El NC2 ,  SCRPE-ZL  36 1  > ,  DBPRPTC  381  > ,  A  0>;  1  85 1  >,A1<  IS 
C 

EGU I  VALENCE  < HCELL S , N2  > 

i- 


O A  TA  P 1/3 . 1 41 5926 5/ 

CUgENT« J . 0 

NPT  S=  I  FIX*;  360.  /DGRENC  > 

+  WAVE  #  IN  FREE  SPACE  AND  IN  THE  OBSTACLE 

K 0-2  .  **PI=('FRE0/3 .  E8 
K2*K  Ci=*SQRT*;  EPSLNR  > 

*  EPSLN0  =  8.S54E-12 

RN=C-<  W/2 .  > 

THETAN»2.*PI/NCELLS 

DUMMY  1  =  -*;  K  0**2  .  /*;  4  .  *2  .  *P I * FREQ *8 , 854E-1 2  >  >%CUREHT 

*  IF  OBSERVATION  POINT  AT  PARTICULAR  CELL,  PTOBS  =0.  IF  OBSER¬ 
VATION  POINT  IN  FAR  ZONE  AT  SOME  ANGLE  PHI,  PTOBS  =1.  IF 
OBSERVATION  POINT  AT  CENTER  OF  OBSTACLE,  PTOBS  =2. 

IF*; PTOBS  .EG.  1  >  GO  TO  20 
C 

IP*; PTOBS  .EG.  2  :>  GO  TO  50 
C 

C  *  DETERMINE  THE  INCIDENT  FIELD  ON  THE  OBSTACLE  AT  EACH  CELL 
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c 

e 

c  * 
c 


c 

c 

c 

c 

c 

t  00 
1  0 
c 

c  * 
c 
c 
c 

20 


c 

c 

200 

30 

C 

C 

c  * 
c 

50 


C 

C 
4  0 


LOCAT  I  Of  i ,  <  Xr»  ,  Vn  > 

LOOP  MCELLS  TIMES  FOR  INCIDENT  FIELD 

DO  10  I*1,NCELLS 

THETA-*.'  1-1  >*THETAN 
XNN~RN*C03< THETA  > 

YNN=RN=*  S 1  N<  TMETh  > 

DUMMY  0~<  Vi  IN—  0 .  >*<YNN-0.  > 

A RGH=K 0=*SORT< >;  XNN-AA  >*< .VNN-AA  >«-DUMMv' 0  > 

CALL  BESEL'i  0,  ARC-H,  BB,  Y,  BP,  YP  > 

J  0=BB 
Y0»Y 

Lift  I  TE<  6  /  5  :>  d0,Y0,BS 

FORMAT','  1 X  ,  "  FLDNC  :  d  0-  " ,  2X ,  1 E 1 5 .  S ,  5X ,  -  Y  0  -  “  ,  2X  .  i  £  1  5  .  S , 

1  "BB=“,2X,  1E15.S,/> 

E I NC  0=DUMMY  1  *CMPLX>;  J 0  ,  -Y 0  > 

VIN< 1 >=EINC0 

WR1TE<  1  00  >  E INC0,  I ,  VINC  I  > 

FORMAT  < IX, “E INC  0=  * , 2X, 1 E 1 5 . 9, 2X, E 1 5 . 9 . 5X, 14, 

1  5X,  "VIN=“,2X,  1 E 1 5 . 9, 2X,  E 1  5 .  V,  S  ') 

CONTINUE 
GO  TO  40 

DETERMINE  THE  INCIDENT  FIELD  AT  THE  FrR  ZONE  POINT  DUE  TO 
THE  CURRENT  FILAMENT  AT  RAU  PRIME.  USE  LARGE  ARGUENEMT 
ASYMPTOTIC  EXPANSION  FOR  THE  HANKEL  FUNCTION. 

DO  30  J=  1  ,  NPT'S 

PHI=<  J-1  >*DGRENC*PI/130. 

ARGH*K0*< R-AA*COS< PHI >  >-PI/4 . 

D 1 =COS<  ARGH > 

42= -S I H<  ARGH  > 

D3=S0RT <  2 .  /<  K0:*PI  =»R  >  > 

E I  NC  1  <  d  >=DUMMY  1  *D3*C MPLX>;  D 1  ,  D2  > 

URITE<6,200>  EINC1< d>, d 

FORMAT< IX, "FLDNC  :  EINC1»” , 2X, 1 E 1 5 . S, 2X, E 1 5 . S, 5X, u  I  =  “ 
1  14, /> 

CONTINUE 

GO  TO  40 

DETERMINE  THE  FIELD  INCIDENT  AT  THE  CENTER  OF  THE  OBSTACLE 

ARGH=K  0*ABS>:  AA  > 

CALL  BESELC  0 ,  ARGH ,  BE: ,  V ,  BP ,  YP  > 

J0*BB 

Y0=Y 

EINC2«DUMMY1*CMPLX<  dO, -Y0  > 

RETURN 

END 

END* 
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T  T-  0  0004  IS  ON  CROC*  079  US  I  NO  0001  0  BLkS  R*0OC<0 


FTfil ,  L 

SErtoC  BL  Kriri,  0  > 

C  '<•  ♦  ■+ : ♦  :  *  ■*■  ■+'•* '*  Hi  :f  i¥  -«*  ■¥  H»  •*  ■*’  :*  •*'  H«  *tf  •■¥  '■**  ■*  ■¥ 

SUBROUT I  HE  ECOPT*.  XrwX  > 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


*f;*  :)•  * :*•*:(<*> ^ I*.****** 

*  THIS  SUBROUTINE  PLOTS  THE  NORMAL I2E0  ECHO  WIDTH  PER  WmVElEMGTH 
VS.  ANC-Lt  PHI  ON  A  LINEAR  PLOT. 


INTEGER  NCELLS .  NPTS ,  P T OtT'S ,  IPHI , riCELLS 

REAL  XNN , VNN , EMAG I , Ji , J0,K0,K2,Y0, VI ,EMAGN, 
i  ECHUFW 

COMPLEX  CfiN.ElNC* , E IhCO, EINC2, 

1  ESCAT, 

2  T  AU  t  , ALNDA , ALPHA , COEF , PNC , C 1 , C£ , V , v  t , V2 , 

3  THU , V IN , VOUT , ONE , I ERO , H  0 , A 1 

D I  MENS  I  CiN  I  BUFRC  1  6  > 

COMMON  S&LKm/  AA , B ,  C, R.EREO,  PERM ,  DG'RENC,  EPSLnR, AN,  NCELLS, 

1  E I NC 1  <  36  0  > ,  V I  N*:  1  05 1  > ,  ESC  A  T  <  36  0  • , 

2  ECHWPWC 36 0 > , K 0 , K2 , VOUT C  1  S5 1  >, T HOC l 65 1  ) , 

3  El NC2 ,  S CRRE2C  361  >, DBRRPTC  36  1  >,h 0<  i  S5  i  > ,  A 1  <  1  $5 1  > 


EOU I VALENCE  C  NCELLS , N2  > 


DATA  P I /3 . 1 4 1 59265/ 
NX OP =2 . *  1 000 
NY OP* 1 .75*1000 

* 

NXSTP=  0  .  S'*  1000 
NYSTP* 0.5*1 000 

NXSTPS=9 
NYSTPS* 1 0 

IPEN= 1 

WRITE*:  25, 3  >  I  PEN 


NXEP=NXOP+NXSTPS*NXSTP 
WR I  TEC  25,1  >  NX  C<P ,  NY  OP 
WR I TEC  25 , 1  >  NXEP , NY OP 
WRITE*.'  25, 2  > 


NYEP=NYOP+NYSTPS*NYSTP 
WR I TEC 25 , 1 >  NX OP, NY OP 
WR I TEC  25 , 1  >  NX  OP , NY EP 
WRITEC 25, 2 > 

DO  100  J=1 , NXSTPS 
I X=NX  0P+  U*NXSTP 
IY1-HY 0P-50 
IY2*MY0P+50 
WRITEC 25, 1  )  IX, I VI 
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WO-  O*  O  C*  C*  O  l\J 


Wft 1 TE<  25 , 1 >  IX, IY£ 
we-ITE*;£5,£:> 

) v 0  CumT 1  HUE 

C 

DO  £00  1=1, HYSTPS 
I  Y=NY  0P+  I-*-HYSTP 
IX1=NX0p-50 
IX£=HX0P+50 
WRITE<£5,1>  1X1, IV 
WRITE-;  £5,  1  >  IX£,  IV 
WRITE*;  £5  ,  £  > 

00  CONTINUE 

CALL  eclbl 
CALL  LABEL*;  i  > 

HPTS«* IF I X< 1  &  0 . /DGRENC > 

IPEN=3 

WRITE*; 25 , 3 >  IFEH 
DO  400  IPHI=i ,  MPTS 

IX=IFIX<>;  i'RHI-1  ,  >*DGRENC/£0,:«NXSTP+££t00.  > 
IY*IFIX<  <  ECHWRUi*;  I PHI  >/XMMX  >/0  .  1  *MYSTP+  1 750  .  > 
WRITE*; £5,  1  >IX,  IV 
00  CONTINUE 

WRI  TE<  £5 ,  £  ') 

FORMAT*;  “  PA  ",1 5  -  ,  "  1 5 ,  “  ;  RD  “  > 

FORMAT*;  "PLi"  > 

FORMAT (  “  I N ;  S'P  " ,  1 1  > 

END 

END* 


* 
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J  O  O  O  (J  O  O  O  O  O  *_>  o  o  o  o  o  o  o  o  o  o  o  o  o  o 


7  7  =  00004  IS  OH  CP: 00 OSS  USING  00004  ELKS  P=  0 0 0 0 
FTN4 , L 

*EHh<  BLKMM,  0) 

C  ■*-*'•*■*■*  -•*  ;f  :*  :*  ■+  :*  :*  •'*  :*•  :*:*  *  :f  f  ;f 

SUBROUTINE  FLDCT 

^  =*♦ :*  •*■»!**:*  s* :« :*:*:*:*  i*:*:*:*:#  ♦:*:*  :*:*:*:*  .*  .■* 


*  THIS  SUBROUTINE  USES  THE  MOMENT  METHOD  TO  DETERMINE  THE  SCATTERED 
FIELD  FROM  SOME  DEFINED  OBSTACLE 


INTEGER  HCELLS, NETS , PTOBS , IPHI , NCELLS 

REAL  XNN ,  YHN ,  EMAG  1,01,00,  K  0 ,  K2 ,  Y  0 ,  Y 1  ,  EMaGM  , 
1  ECHWPU 


COMPLEX'  CMN,EINC1  . EINCO, EIHC2, 

1  ESCHT , 

2  ,  T AU 1 , ALMDA , hL PHA , C OEF , F AC , C 1 »C£ , V , V 1 , V2 , 

3  TAU , V I M , VOUT , ONE , EERO , DUMMY® ,  a  0 ,  A 1 

DIMENSION  IBUFR<16> 

COMMON  /BLKMM/  AA, B, C, R,  FREQ,  PERM  ,  DGREnC,  EPSlNR,  AN,  NCELLS , 

1  EINCl<36CO,VIN< 1S51 >, ESCaT< 360 >, 

2  ECHWPWC  36  0  > ,  K  0 ,  K2 ,  VOUT <  1 25 1  > ,  TAU1;  1  S 5 1  > , 

3  El NC2 ,  SCkPE2\  361  > ,  DBPRPT-.  36 1  > ,  A  0%  f  25 1  > , A  K  f  S5 1  > 


EQUIVALENCE  < NCELLS, N2> 


DAT  A  P I /3 . 1 4 1 592652 

*  THE  LARCE  ARGUEMENT  ASYMPTOTIC  EXPANSION  FOR  THE  HANKEL  FUNCTION 

IS  USED.  SEE  RICHMOND  - 

*  LOOP  NPTS  TIMES  TO  OBTAIN  THE  SCATTERED  FIELD  AT  EACH  DGRENC 


HPTS= I F I X<  36  0 . /DGRENC  > 


RN=C:-I  W/2  ■  > 

THETAN=2 . *P I /NCELLS 

DO  20  I  PHI- 1, NPTS 

PH  I =•:  I  PH  I  - 1  >*DGRENC*P  1/ 1  SO . 
DUMMY 1 - <  K  0 «R >-< P I /4 . > 
C>t>mY2--C0S<  DUMMY  1  > 
DUNMY3=-S I Hv DUMMY  1  > 

DUMM Y4=SQRT<  P I  *K  0*  0 . 5/R  > 

*  LOOP  NCELLS  TIMES 

DUMMVS=< 0 . , 0 .  ) 

C 

DO  10  N* 1 , NCELLS 

THETA-IN- 1  >*THETAM 
XNH-RNkCOS-I  THET  A  > 
YNN«*PN*S1N*I  THETA  > 


91. 


S>Ui1MY5=K  i"m  •;  *-:NHfrOS<  PH  I  >+  S'HM+S I  H<  PH  1  >  > 
I»UriMY6=f-0S<  L'-uitr;  v-ei  > 

C'UfiMN'7 -f  1  N\  DUrlKYS  > 

AfcG  J 1  =KCi»H>i 

C  hLL  REBEL (  1  ,  rtftG  J 1  , BB , V , BP ,  VP  > 

J  f  =BB 

t>UHMYS»BUHMV8+<  EPSLNR-1  .  >»VOUT<  H  >**»■!■*  J1*CMPL 
1 0  COMT I HUE 

C 

ESCrtT*:  IPHI  >*DUMMY4=»C3,iPLX<  DUMMY3 , DUMrt V*  >:«-DUrl(t 
C 

20  CONTINUE 

C 

RETURN 

EnT> 

EN£>* 
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0000000  00000000  0  00000  0  00  0  Of 


FTH4,l. 

♦EMA*'  FLKMM ,  0  > 

C  +  +  +  V  f  t  :**  *:*  *-  4-  *  f  +  +  M  *.  >  -f 
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C: 

c 

C  *  From 
C 

C  "Hr.t.inna  Theory  and  Dor  igh“ 

C  Warren  L .  'hi  utzmar.  and  Gary  A.  Thiele 

C  John  Wiley  &  Sons,  New  York,  ISS'i 

G  Appendix  G.7  pp.  57S-5S i 

C 


PURPOSE 

TO  SOLVE  A  SYSTEM  INVOLVING  A 
ONLY  5N  STORAGE  LOCATIONS  FOR 


T0EPLIT2  MATRIX. 
AN  N  BY  U  MATRIX. 


TFl2  PEC'UI 


*  REMARKS' 


A  toeplitz  matrix  has  the  first  row  equal  to  the  first  column. 
All  elements-  along  the  main  diagonal  are  equal.  A.-.y  diagonal 
off  the  main  diagonal  will  have  this  same  property. 

■*  DESCRIPTION  OF  PARAMETERS 


..  N2 
'  TAU 


AG,  A1 
VIN 


MM 

XNORM 

IER 


-ORDER  OF  MATRIX 

-FIRST  ROW  OR  COLUMN  OF  THE  T0EPLIT2  MATRIX  < VECTOR 
LENGTH  N2i 

-VECTORS'  OF  LENGTH  N2  NEEDED  FOR  SCRATCH  AREA 
-FOR  THE  MATRIX  EQUATION  <2X1  >*<V>,  VIN  IS  V.  <2,1, 
*  AND  V  MAY  BE  THOUGHT  OF  AS  GENERALIZED  IMPEDANCES, 
CURRENTS',  AND  VOLTAGES,  RESPECTIVELY!.  V  IS  A  N2  BY 
MM  MATRIX. 

-NUMBER  OF  COLUMN  VECTOR'S  ON  THE  RIGHT  SIDE  OF  MATRIX 
EQUATION  <  2  >•;  I  >=•;  V  >  <  USUALLY  1  >  . 

-UPON  RETURN  THIS  IS  INFINITE  MATRIX  NORM  OF  INVERSE. 
-ERROR  CODE  WHICH  IS  0  IF  NO  TROUBLE. 


C 

C 

c 

c 


INTEGER  NCELLS ,  MPT  S' ,  PT  OBS ,  IPHI,MC£LLS 

REAL  XNN , YNN , ENAG I, J1 , JO, KG, K2 , Y G , Y 1 , EMAGN , 

1  ECHWPW 

COMPLEX  CMN,EIMC1 , E INCG, EINC2, 

1  ESCAT, 

2  TAU1  ,ALMDA,  ALPHA,  OOEF, FAC, Cl  ,  C2 ,  V ,  V 1  ,\>2 , 

3  T  AU , V I N , VOUT , ONE , ZERO , A  0 , A 1 

D I MENS I ON  I BUFR< 1 6 > 

COMMON  XBLKMM/  AA ,  B ,  C,  R ,  FREQ ,  PERM ,  DOPE  NO ,  EPS'LNR ,  AN ,  NCELL 

1  EINC1<36®>, VIN< 1S51  > , ESCaT < 3S ft > , 

2  ECHWPW< 36 0 >, K  G , K2 , VOUT < 1 05 1 ), 1 AU< 1S5J  > , 


• J  O  Q  O  O  —  —  O  O  >N  •->  O  O  O  O  O  O  O  O  O 


E1NC2,$CF:PE2<3£1  >,C>BPRPT<  361  >,A0<  '£51  >,A1<  1S5i 
EGU I  VALErJCE  <  NCEi.LS ,  NZ  > 


DAT  A  OMtA  1 E  0 ,  OE  0  V ,  ZERO  A  Ot  0 ,  0£  0  ')/ 

DAT  A  ONNE/ 1 E  0/ ,  2R  f'G/  0 .  E  0/ 

UR I TEC  6 ,  S 0  >  HZ,  MM 

FORMAT*;  IX,  “N2=,,,2X,I4,5X,“rtM*“  ,2X,  14,/) 

DO  15ft  1=1. HZ 

WkITE< 0,  1  00)  TAU*.  7  >,  VIN*.  I  ),  i 

FORMAT*;  IX,  “TPLZ  :  TAU=  "  ,  2X ,  1 E 1 5  .  S ,  2X,  E 1 5  .  £  .  5X,  "VIH=“,2X,  1E15. 
2X,E15.S,5X, “1  =  “,  14 ./) 

COHT I HUE 
N=H2- t 
1ER=0 

*  NORMALIZE  INPUT  MATRIX 

TAU1=TAU<  1  ) 

DO  200ft  11=1, N 

0  00  T  AIK  1 1 >=T AU< 1 1 + 1 )/T AU i 

*  the  FOLLOWING  CALCULATES  THE  ITERATIVE  VARIABLES  TO  OBTAIN 

ao*;n>  and  almda 

*  NOTE  ;  VECTOR  A0<  I  >  HAS  I  ELEMENTS  AND  IS  STORED  AS  AO*.  i,J>, 

0=1  ,N 

ALMDA=ONE-TAlK  1  )*«TAU<  1  ) 

AO*;  i  >=-tau*;  i  > 

1=2 

1  KK= I - 1 
ALPHA=ZERO 

DO  2  M= 1 , KK 

LL*=I-M 

2  ALPHA=Ai.PHA+AO< M  >*TALK  LL  > 

ALPHA=-*;' ALPHA+TAU*;  I  )) 

I  F>;  CABS*;'  ALPHA  >  .  EC; .  0 .  D  0  >  GO  TO  15 
COEF=ALPHA/ALMDA 
ALMDA=ALriOA-C:OEF**ALPHA 
DO  Z  0=1, KK 
L=  I  -  J 

Z  A 1  *;  J  )=A 0<  J  J+COEF*  A  CK  L  ) 

DO  ?  d= 1 , KK 
7  ACK  0  >=A  1  <  J  > 

A0<  I  )=COEF 

IF< I  ,GE.  N)  GO  TO  5 
4  1=1+1 

GO  TO  1 


UO 

50 


*  THE  FOLLOWING  COMPUTES  VALUES  OF  EACH  ELEMENT  OF  THE  INVERSE 


NH=*;  NZ+ 1  >/2 
FAC=ALMDa+T AU 1 
XMORM=ZRRO 
NP=HZ+ 1 
DO  51  1=1, HH 
XHM=2RF;0 


r>  c-  th  r>  otio  m  o  *>  r*  w  r*  o  o 


if-:  i  .he.  i  ;  co  to  52 

ft  1  *'  1  >--OhE  FrC 

XNh=Cft&S'.  HK  !  >  > 

DO  53  J=2.N2 
A1<  J>=ftO*.  J-1  .VFfiC 

53  X* <M=CABS*C «!>;  J > >-*-XNM 

GO  T  0  54 
52  XNM=ZkRO 

JH= I - t 
C1=A0*;'JH) 

NNPI=NP-I 
C2=Au'i*:NHPi  ) 

DO  55  JJ-1 ,N 
J-HP-JJ 
INPJ=WP-J 
JL- J-1 

ft  I  *;'  J  >=«  1  >;  JL  >+*:  C  l  *ft  CM  JL  >-C2-*'ft  CM  i  HP  J  >  «C 

55  XMM=CABS>;  ft  1  *;  J  >  >*HNr? 

ftl*:  1  >=A£m  1-1  .V’FAC 
XNM=XNH+C'ABS*.  ft  1  •;  1  )> 

54  IF<XNM  ,GT .  XNOF.fi  >  XN0RM=XNM 
C 

■*  MATRIX  MULTIPLY 


25 


50 


51 


DO  56  11=1, MM 
ir>*<  i  i-i  >+N2 
V=ZERO 
V1=2ERCi 
DO  5?  J= 1 , NZ 
HIDJ*ID+J 
V2=V  I  N<  N 1 D  J  :> 

v*v+V2*a  i  •:  j  > 
KNP J=MP- J 


V 1 = v 1 + V2**A 1 < KNP  J  > 
NID1--ID+I 


VpUT>;  NIDI  >=V 

UP  I  Tl<  6 ..  225  >  VOUT*;  N I  D  I  > ,  N I D I 
FORMAT*:  IX,  "TPLZ  :  VOUT*" ,2X,  1  El  5 . 3, 2X, E 1 5 . 8, 5X,  *‘NI  = 
KIDNP I  =  JD+NP-I 
VOUT >:  K I  DNP  I  >=V  1 

•  UR  1  TE<  6 , 250  >  VOUT >;  KIDHPI  >,  KIDHPI ,  VOUT*;  NIDI  >,HIDI 
FORMAT' 1 X , " TPLZ  :  VOUT* “ , 2X ,1E15.8,2X,E15.S, 5X , “K I • 
3X ,  "VOUT  =  " ,  2X ,  1 E 1 5 . 3 , 2X , £  1 5  -  3 ,  5X ,  "NI=",2X,  14 
WRITE<  6, 251  >  VOUT< KIDHPI  >, KIDHPI 
FORMAT< 1 X , "TPLZ  :  VOUT* “ , 2X , 1 E 1 5 . 3 , 2X , E 1 5 . 3 , 5X , “ K 1  * 


“,2X, 14, X. 


“,2X. 14, 
,  X  > 


:-x,  14,  x> 


1  CONTINUE 


DO  650  1=1, NC'ELLS 
WRITE*;  6. 6 00  >  VOUT*;  I  > ,  I 

00  FORMAT < IX, "TPLZ  : :  VGUT=“,2X, IE  1 5 . 8, 2X, El  5 . 8, 5X, " I * - , 2X , 14, X> 

50  CONTINUE 

RETURN 

15  WRITE*:  6,700) 

700  FORNftT< IX, "ERROR  HAG  OCCURRED.  MATRIX  IS  STRONGLY  MOMS  JUGULAR “  ) 
IEF.=  I 
C 

RETURN 


cn  o  m  o  o  o  o  o 


/ 


L  T=00i'03  IS  OH  C K (i 0 07.9  001  Hi  0001;'  ELKS  R=0000 
FTN4,L 

SUBROUTINE  BESEL*.  N,  X, B,  Y,  BP, VP  > 

C 

COMMUTE  THE  BESSEL  FUNCTION  OF  ORDER  N.  M  «N  INTEGER  N 
WITH  REAL  ARGUMENT 

ALL  EQUATION  REFERENCES  TO  aB.RhHOWITE  AND  OTEGUH 


D I  ME  NS  I  ON  C  0<  ?  > ..  C 1  <  7  > ,  D  0*.  7  > ,  0 1  ■;  7  > ,  E  0<  7  > ,  E 1  \  7  >  .  0  0<  7  ■> ,  G 1 
DATA  CO/1 . 0.-2.2499997, 1 . 2656203 , - . Si 63666 . . 444479E-1 , 
*- . 39444E-2, .21E-3/ 

DATA  Cl/ .  5,  -  .  56249905,  . 2  J  093573,  -  .  39542S9E- 1  ,  .443319E- 
*-  .31761E-3,  .  1  1  09E-4/ 

DATA  DO/  .  79700456,  - .  77E-6,  -  .  55274E-2,  -  .  95 i  2E-4 ,  .  1  37277 
■*-.  72605E-03,  .  f  4476E- 03/ 

DATA  D1/. 79700456,  .  156-E-5,  .  t  659667E-1  ,  .  171  05£-3,-.2495 
* . 1 1 3653E -  02 ,  -  .  2 0 033E -  03/ 

DATA  E0/- .  7053901 0,  -  .  4  I 66397E-J  ,  -  .  3954E-4 .2625 73E-2,- 
. 29333E—  03,  .  13S5CE-03/ 

DATA  E 1/-2. 356 l 9449,  .  12499612,  . 565E-4 , - . 637079E-2,  .743 
*  ,  79024E- 03 ..  — .  OR  1  6©E  —  03/ 

DATa  0  0/ . 36  74  669 ,  .6  055937 ,-.7435  030 ,  . 253  0 0 1 2 , - . 4 20 1 2 1 E 
* .42791 6E-02, - . 24O46E-03/ 

DATA  G1/- .6366190, .2212091,2. ! 60271, -1 ,310403, . 3 1 2395 , 
*-.400976E-01  ,  . 27S73E-02/ 

DATA  P I /3 . 1 4 i 5926/ 

IFLG-0 

IF  <N,LT. 0)  IFLG=1 
N=IABS<N> 

IF<  ABS<  X  > ,  LT .  I  .  0£~  1  0  >60  TO  I  5 0 
I F<  ABS<  X > , GT . 3 , 0  >G0  TO  50 
X3S  u  -  X;*  X/9 .  0 
PROD= 1 . 0 
B  0=  0  . 0 
B I  =  0*.  0 
CUM  0=0. 0 
CUM 1=0,0 

SEE  EQUATIONS  9.41  AND  9.44 

DO  5  1=1,7 

B0=6  0+C0<  I  >*PROD 

B 1  =B  1  +C 1  <  I  >*  PROD 

CUM  0=CUM  0+G  OK  1 )*RRGD 

CUM  1  =CUM  I  +G 1  <  I  ):*PRGD 

PR  0D=Pft  OD:*  X3  SQ 

CONTINUE 

B1=B1*X 

XC=2 . 0!*SNGL<  DLOG<  DBLE<  0 . 5*X  >  >  :>/P  I 

V  0*XC*B  0+CUM  0 

V I  =XC:*B  1  +CUM 1  /X 
GO  TO  1 00 

EOS  9 i 4 , 3  AND  9.4.6 
0  THR0VX=3 . 0/X 
PR0D= 1 . 0 
F  0=  0 . 0 
F  1  =  0 . 0 
THETA 0=X 
THETA 1=X 
DO  33  1=1,7 
F  u=F  0+D0<;  I  >=*PRGD 


;0R=  0 

<  7  > 

e-2 , 

1  IE-2, 

. 54 1 25 
4SE-3, 
-01  , 
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F  i  =F i +D i <  1 >*PkOu 
THETA  C<=THS7A 0+E  f.K  I  rvPFOD 
THETA  i  =TMETA  1  -» E  i  ■:  I  >'PkOD 
PROD=PF;OD  *  THRGVX 
CONT I HUE 

SQRX=  1  .  G/SNSL<  DSGRT-:  OBLEC  X  >  >  > 

B 0=SC;kX+ F 0:«SnGL-;  DCOS<DBlE<  THETAu  >>  > 
Bl*SC;F.tf*F  1  *3NGl< DCOSl DBLE-.  THETa 1  >>  > 

V 0*SOF>X*F 0*SNGLC  DSIN<  DBLE-;  THETA  0  >  >  > 

V  t  =SGRX*F  1  *SNGL<  DSIN<  DE-LEC  THETA l  >  >  > 
IF<  H- 1  >1  CU  ,  1  05, 1  1  C* 

B»B0 
£:p=-B  l 
Y=Y0 
VP=-YI 
GO  TO  £00 
B-B1 

BP=B  0-B 1  XX 
Y*Y« 

YF-Y0-Y1XX 
GO  TO  200 

FOP  RECURRSIVE  DIRECTION  COMMENTS  SEE 
XH=N 

I F< XN . LT . ABS<  X >  >GO  TO  130 
FOR  X<  N  RECUR  DOWNWARD 
BLAST* (.0 
BLASTP=  0 . 0 
U*N«- 10 

DO  115  1=1  ,  J 
XI  =  J-I 

BNEXT*2 . 0*XI*BLAST/X-BLASTP 

BLASTP=BLAST 

BLAST  =BNEXT 

IF< I .HE. 1 0 >Gu  TO  115 

BLP=BLASTR 

2-BNEXT 

CONTINUE 

IF<  AB3<  B0>  .LT  .  ABS-1B1  >  >GO  TO  11? 

CORR*fi O/BLASTP 
GO  TO  IIS 
CORK— B1/BHEXT 
B=BLP'«C'ORR 
BNNIN1*2»C0RR 
GO  TO  140 

FOR  X>N  RECUR  UPWARD 

BLASTP=B0 

BLAST=B1 

DO  135  1=2,  N 

XI* 1-1 

BW£XT=2 ,  0*X  I  :*BL AST XX-BL ASTP 

BLASTP=BlAST 

BLAST-BNEXT 

CONTINUE 

B=BLAST 

BMMIN I  * PL ASTP 

Ftp -- BUM  IN  I  -XN+BXX 

YLASTP--Y0 

VLAST*Y1 

DO  M5  1=2, N 

X 1  =  I  —  1 


SEC  i  I  ON  •?.  12,  P3S5 


Appendix  Q.  Computer  Program  for  Moment  Method 


This  appendix  contains  the  computer  program  which  implements  the  moment 
method  solution.  This  program  is  a  collection  of  routines  written  in  FORTRAN 
4.  Each  subroutine  is  self-explanatory;  the  road  map,  RDMAP,  gives  some  indi¬ 
cation  of  how  the  final  answers  are  arrived  at.  The  computer  system  used  was 
the  Hewlett-Packard  27MX  sertes  microcomputer. 


oooooo  o  o  o  o  o  o  o  o  O  O  O 


T=  0 0 0 Cm  IS  ON  CROC OSS  USING  00  CBS  BLKS  ft=0000 

FTH-* ,  L 

CMBl k' ,  0  > 

0  * 

PROGRAM  KHL'Tl  I 

C  *  *■  ♦  45  :*  s*  ’*■*  :*:*:*:*  :f :* :4 

C 

C  ROBERT  K.  SCHNEIDER 

C  SMtSSSMSSMMStKa* 

C 

C  *  THIS  PROGRAM  EVALUATES  THE  EIGENVALUE  SOLTI Oh  EOF:  THE  RaDIATED  FIELD 
C  OUTSIDE  AH  INFINITELY  LONG  M  ELECTRIC  CYLINDRICAL  Shell  OF  I*-'NER 
C  RADIUS  B  AND  OUTER  RADIUS  C.  THE  INCIDENT  FIElD  HAS  sS  ITS  SOURCE 
C  AN  INFINITE  LINE  CURRENT  ClOSE  TO,  BUT  EXTERNAL  TO,  "'HE  CYLINDER 
C  AND  AT  kN  ANGLE  PH  I  -SEINE  EQUAL  TO  EERO, 

C 

C  '■*  A  MATRIX  EIGENVALUE  SOLUTION  HAS  SEEN  FORMULATED  LEAVING  THE  ‘OS'rFl- 
C  Cl  ENTS  OF  THE  F  I  ELD  TO  EE  SOLVED  FOR.  AN  ALGORITHM  Sr'  PRESCOTT  D. 

C  GROUT  IS  USED  TO  SOLVE  FOR  THESE  COEFFICIENTS. 

•'*  THE  FAR  FIELD  IS  GIVEN  BY  : 

E'en  =  -jwuFn  TIMES.'  THE  HANKEL  FUNCTION  OF  T.iE  End  KIND 
uF  ARGUEMENT  k  0;*RhU  AND  ORDER  n , 

THE  TOThL  FIELD  IS  a  SUM  ON  r,  OF  THE  ABOVE.  A  CONVERGENCE  OF  THI 
SUM  IS  DETERMINED  BY  THIS  PROGRAM.  ALSO,  SINCE  THE  mEGUEMENT  OF  TM 
ABOVE  HANKEL  FUNCTION  IS  »  n,  THE  APPROPRIATE  ASYMPTOTIC  EXPANSION 
IS  USED. 

*  EPSLN  IS  A  VALUE  INPUT  BV  THE  USER  FOR  USE  IN  DETERMINING  CONVERGENC 

*  KO  =  WAVE  NUMBER  IN  FREE  SPACE 

* 

*  KB  =  WAVE  NUMBER  IN  DIELECTRIC: 


C 

C 

C 

c 

C 

c 

c  * 
c 

c 

c 


INTEGER  E ST  HIE 

REAL  K  0 ,  KB ,  MAX ,  NRMFLD ,  NPH I  ,  V C> 

COMPLEX  HNB0A ,  NMPB  DA ,  SUM  EG ,  A<  6 , 7  > ,  F<  S  .) ,  FNC  ICC)  > ,  D  )  1,012, 
1  SUMM , SUMP , SUM , EG  ON , DC 

D I  ME  NS  I  OH  I BUF  R<  I  S  > ,  JB<  1  0  > 

COMMON  /C'MBLK/  A,  F  ,  ECHWPWC  30  CO  >,  DG'RENC,  SCRPEZC  3S0G  >, 
t  DBPRPT <  36  C<  0  > 

DATA  P I /3 . 1 4 1 5026,5/ 

DATA  XCR/4 . S?5/, YC R/3 . 875/, RDS/B . 748/ 

SET  UP  CALL  FOR  INPUT  FROM  DISK  FILE  "DATA  1  11 

DATA  IBUFR/2*  0 , 1'HDA,  BHTA  ,  BH1  ,  3*  0 , 22 1  & ,  7+  (</ 

CALL  SPOPNC  I  BUFF: ,  ISLU> 

C:aLL  EXEC*:  XL,  1  ) 
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kEaC'C  ISLU, *  >  h*>  B .  C  ,  Cur'EN  ' , R, EPSL  N ,  FRE'*' . 
l  ESTIr.E,.  fi,  IS.'-.  IP 

SET  DEFAULT  FOR  ESTIHE 

IF>;  ESTIHE  .  E&  .  0  )  ESTIME=f  0  0 

2SPO=0 . 0 

CALCULATE  THE  WAVE  NUMBERS 

K  0-2 . 0*P  I  »r  RECS/3 .  ES 
K2=K  &  ♦’SORT >:  EPSlNR  > 

COMPUTE  ARGUEMENTS  FOR  BcSEL  FUNCTION  CALLS 

X 1  ~K  0*AA 
>•  2-k  u*& 

X3=K0*C 

X4»K2<*B 

Xz~k2*C 

LOOP  ESTINE  TIMES  TO  DETERMINE  THE  EIGENvALl 


DO  100  L=-ESTINE,ESTIM£ 

K-K+1 

EPS'LNW*  1.0 

LM1=L-1 

Lr11=L 

CALL  BESELC  LM 1  , X 1  ,  BB ,  Y ,  BP ,  Yp  > 
HNS  0A*CMPL  X-:  Bo ,  -Y  ) 

HNP2  0A=CriPL  >'.<  BP ,  -YP  > 

«>,'  f ,  1  >=CMF*LJ'X  BB ,  ZERO  ) 

&<  1  ,2  >*CHPLX'i  Y ,  ZERO  ) 

A-;  1 , 3  >*<  0 . ,  0  .  > 

A<  1,4  >='.  0 .  ,  0  .  > 

•  A>;  1  , 5  >=<  0 .  ,  0  .  > 

AC 1 , 6  )=-HN20A 
AC  1 , 7 >=<  0  .  ,  0 .  > 

AC  to,  1  >=CMPLXC  BP ,  ZERO  > 

AC  6 , 2  :>-CMPLX<  YP ,  ZERO  > 

A1:'  6 , 3  >=<  0 .  ,  0  .  > 

AC  6 , 4  >-<  0 .  ,  0  .  :> 
a*:  s ,  5  >=■;  o . ,  o . ) 

A1;  S ,  <S  >*-HHP2  OA 

IF< LM 1  . EO .  0>  EPSLNN=  0 . 5 

EGl LWT=C  CUREHT *-EPS'LNN  VC  0  > 

A';  6 , 7  >*CNPLX<  EGLMT  ,  ZERO  > 

CALL  BSSEU  LM  1  ,  X2  .  SB ,  Y ,  BP ,  YP  > 
DUM1  =-«(<♦  BP 

AC  3 , 5 >=CM PL X< - BB, ZERO  > 

AC  f. ,  5  >=CHPl.X<OUM  1  ,  ZERO  > 


tt-'SLNK  . 


iCH  ORDER  n  . 


CALL  BESELn  LM  1  ,  X3 .  BB ,  V ,  BP ,  VP  ,■> 

DUr-1 1  =K  fr*tK 
DUM2»r.  0*  VP 

An  2 ,  i  >=CriPLXn  BB,  EERO  > 

An  £ ,  £  >-CMPlX<  V ..  2ER0  > 

CC 

A<  4 ,  1  >=CMPl.Xn  DUM 1  ,  ZERO  > 

A>;  4 , 2  >=CriPLXn  uUr;2,  ZERO  > 

C 

CALL  BESELn LM 1 , X4 , BB , V , BP , VP  > 

DUril  -K2-+PP 
£>UM2=K2‘*YP 
An  3 ,  1  >=n  0 .  ..  C< .  > 
a*;  3 , 2  >=•:  o . ,  o .  > 

An  3 , 3  >=CMPL;<X  BB ..  ZERO  > 

A*;  3,4  >=CMF'LXn  V,  ZERO  > 

A<  3 ,  6  )=•:  0 .  ,  0 .  >  ■  t 

A<  3 ,  ?  >=•;'  u . ,  0 .  > 

CC 

a<  s ,  i  >=•;  0 . ,  0 .  > 

An  5,  2  >*».’  0 .  .  0  .  > 

An'  5 , 3  >*C MPLXn  DUM 1  ,  ZERO  > 

M 5,4  >=CMPL)« 0UM2, ZERO > 

An'  3 , 6  >=<  0  .  ,  Ci .  > 

An  5 , 7  >=<  C< .  ,  0 .  > 

C 

CALL  BESEL<  Lf!  1  ,  M3 ,  BB ,  V ,  BP ,  VP  > 

DUN  1 ~-K£*8P 
DDKS- -K2+ YP 
An  2 ,  3  >~CMPLX< -BE: , ZERO  > 

AC  2 , 4  >=CNPLXn -V , ZERO  > 
a<  2 , 5  >=•:  o . ,  o .  > 

An  2,  6  >=n  0 .  ,  0 .  > 

A<  2 , 7  >=n'  Ci .  ,  0  .  > 

CC 

A«*4, 3  >-CNPLX<  DUN  1  ,  ZERO  > 

A<  4 , 4  >=CMPLX'n' DUN2, ZERO  > 

A<  4 , 5  >=<  0  .  ,  0  .  > 

An'  4 , 6  >=<  0 ,  ,  U  .  > 

AC  4 , 7  >='.'  0 .  ,  0 .  > 

C 

C  *  MOW  THAT  THE  MATRIX  A  IS  FORMED,  CALL  CL SKY  TO  SOLVE  FOR  THE 
C  EIGENVALUES. 

C 

M-M+l 

CALL  CLSKYC  N , M  > 

C 

C  *  PULL  OUT  Fn  FROM  ARRAY  F 
C 

IF<  K  .ME.  ICO  GO  TO  17 

C  DO  IS  1=1,6 

C  WRITE'.' 6,  IS.)  Fn'  I  >,L,  I 

IS  FORMAT n IX, "F=  ",  1 E 1 5 . S,2X,E 1 5 . S, 5X> 14, 5X, 14, X > 

CIS  CONTINUE 

C  WRITE'.'  B,  70.) 

K«=  0 
C 

1 7  FN<  L+EST I ME+ 1 >=Fn  6  ) 

C 

C  DO  BO  M= 1 , B 


10? 


C  WRITS*:  6.  SO  >  F<  Pi  .* .  M  .  L 

5  0  FORMAT  -  1  X ,  *•?*“  .  £N  .  1  £  1  5  .  ®  .  £.' : ,  F.  1  5 . 3  O' ,  I  -?  .  5X ,  I4> 

CO  0  COHT I  ii’.'E 

C  WRITE.  3,  70  > 

70  F ORMaT  1  ;i ,  //  > 

C 

100  CONTINUE 

C 

WR 1  'i  Eu  1  ,  £00  > 

£00  FORMAT1;  IX,  "FN  E IGENVaLUES  CALCULATED  .  "  > 

C 

C  *  CALCULATE  THE  ECHO  WIDTH  PER  WAVELENGTH  RELATIVE  TO  THE  INCIDENT. 
C  FIELD  AT  THE  CENTER  OF  THE  OBSTACLE 
C 

SFLDMX® 0 . 0 

tfldm:-:=  o  .  o 

I  AVNTri=  0 . 0 


1 

2 


3 

4 

r 

54  0 
C 

C 

3 


4  £5 


4£7 


4  £3 
4£3 
550 

r. 


I  riAX=  0 

NT  I  MS 3=  I F I 33  0 .  /OGRE*  !C  > 
DO  7  00  l  =  ),NTlr)E3 
IM)=I-i 

PH  I  =  I M  t  i*DGRENC:*P  I  /  i  S  0  . 
NCNVRG* 0 
SUMEZ=< 0 . , 0 .  ) 

SUMM=<  0 0  .  ') 

SUMP=<  C . , 0 .  > 


DO  550  J=1,ESTIM£ 

JM 1  =  J-  I 

ARGP=  JM 1  *PI,'2 . 

C>1=S.TN*;'ARGP> 

D £~  C  O S  .  ARGF*  }  . 

IF<  JM)  .EO.  0  >  GO  TO  540 
ARGTi=-ARGP 
D3=STN<  A  PGM  > 

*  D4=COS< ARGM  ) 

SUMM-FPK  JM)  >*COS<-JPi1!*PHI  >*CMPLX< D3,  -D4  > 

SUHP=FN<  JM  1 1  EOT  I  ME  + 1  >*COS>:  JM )  *Pk  I  >*CNPLX<  D 1  ,  -D£  ) 


SUMEZ=SUMEl  +£ .  x*P  I  *FREu*PERM*>;  SUMM+SUMP  > 

Dt  0=D£*AIMAG*:FN>:  JM1+ESTIME+)  > >*COS<  JM1*PHJ  >•*£. *PI*FREO*PERM 
IF*.  DIO  .GT.  EP3LH>  GO  TO  550 
HCMVRG=NCNVRG'4  1 
j£‘.  NCNVRG'  :*=  J 

IFTNCNVR G  .LT.  ICO  GO  TO  550 
DO  4 £5  L=£ , 1 0 

IF*.'  J£<  L  ')  .NE.  J£<  L-1  >+l  >  GO  TO  4£7 
CONTINUE 
GO  TO  510 
NCNVRG- I 0-L 
DO  4 £3  L=1 , 1 0 

IF<  L  .GT.  NCNVRG >  GO  TO  4 £3 
J2<  L  >■ J2< 1 0-NCNVRG+L > 

GO  TO  4 £3 
J£<  L  >=  0 
CONTINUE 
COHT I  HUE 


WRITE< 3, 503  > 
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5  06 

c 

CS  1  0 

520 

C 

5 1  0 


53  U 


C 

C 

79 

C 


7 

c 

c 

77 
C 
S 
C 
c 

78 
C 


C 

C 

C. 


9 

t  t 
14 


C 

C 

C 

Cl  000 
C 

c 

700 

C 

C 

C 

1200 


FORMAT*;  IX,  "TRIED  MORE  THmN  £ST1;:£  TERMS  IN  SUh.“,/> 
v.'R  I  j  L-;  6  ,  52  0  }  I ,  J 

FORMAT-:  1 X ,  “  I- “  ,  2:-. ,14, 5X ,  “  J*  "  ,  2:- : ,  I -i  .  /  > 

IF-:  J  .OT.  I  MAX  >  I  MAX*  J 
I AVN  I  r!=  I AVH  fM+ J 
IF-:  I  SKIP  .HE.  t)  00  TO  530 

AROCS=<  *:  K0*R  >— ;  PI/4 . 0)>S<206*PI  > 

GO  TO  7 

ARGCS=*;  K  0-R  >-■;  P I  /4 . 0  > 


WRITE-:  0 , 79  >  AkGC  3 
FORM  AT-;  1  X,  “AftCCS*  "  ,  E 1 5 . 8,  /  > 

D 1 =COS< ARGCS > 

WRITE-:  0,7?  > 

FORMAT-:  1  X ,  “OETT I  NO  HERE  7“  ,  /  > 
02= -SIN-:  AROCS  > 

WR I TE-:  6 , 78  ) 

FORMAT-;  1 X ,  “CiETTI NO  HERE  $"  ,/  > 

D3=S0RT >:  2/-:  P I  *K  0:*R  >  > 
05=r>:v*CMPL»;  D 1 , 02  > 
E20N=EW5*SUME2 


SCRRE2-:  I  >=»CABS<  E2  ON  > 

IF<  SCRPEZ-;  I  >  .  GT  .  TFLDMX  >  TFLDMX«SCRPE2<  I  > 

ARC=K0  «1ASS-:  AA  > 

CALL  PE8EL<  0,  ARC,  BB,  V,  8P,  VP  > 

*  J0=£:8*BB 

Y0*Y-*Y 

D 1  =3 .  E8**4  .  *£;.  S54E- 12**2  .*8  .  A  FREG**2  .  *CUR£HT**2  .  ■* P I * 
ARG=K 0*AA*COS*:  PHI  > 

02=C'0S-:  arc  :> 

03=81 N<  ARC  > 

04=2 .  *P  I  *  FR  E Q * C LIRE N T . •' <  4  .  *3 .  E8**2 .  ♦•S  ,  854H  -  )  2  > 

0 1 2=D4  *'CMPLX-:  02 , 03  > 

01 1 =SUME2+D 1 2 
XMAG«CABS<01  1  >*CABS<  £>  1  1  > 

ECHWPW-:  I  >=D1*XMAGA  JO+YO  > 


WR  ITE-:  0 ,10  0  0  >  ECHWPU-;  I  > ,  PH  I 

FORMAT  < 1 X , “ECHWPW* “ , SX , 1 E 1 5 . 8 , 5X , “ PH I * “ , 2X , F7 . 3 , / > 
IF-CECHWPUK  I  >  .GT.  8FLDMX >  SFLI>MX*ECHWPW<  I  > 

CONTINUE 

■  SCRPt  3<  NT  I MES+ 1  >=SCRPE2*:  1  > 


WRITE*  6,  12.00) 
FORMAT* 1 X , //// ) 
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J  >.  ,  r  ,  .  ^  /  ) 


WRITt*  l  .  23  ' 

FORMAT*  “PLOT  THE  SCATTERED  FIElDI'H  i.'r-t”  .> 
RE  HD1.'  1  ,  I  2  :»  I  AMS 
IF* IANS  .ED.  I H  >  CO  TO  505 
IF* IaMS  .ME.  1HY>  GO  TO  507 

PLOT  THE  ECHO  WIDTH  PER  WAVELN-..TH  V G  ANGLE 

CALL  ECOPT* SFLDMX > 

DO  900  1=1, NT  I ME? 

PHI=* I- 1 . >*DGRENC 
WRITE*;  6,300  ECHUPW*  I  >,PHi 
FORMAT1.'  1 X ,  “  EC'HWPW*  "  ,  2X , 1 E 1 5 . 3 , 5X ,  "  PH  I  =  “  . 
CONTINUE 

WRITE-.'  1  ,29  > 

FORMAT1;  “PLOT  THE  TOTAL  P  I ELOv#  Y/N“  > 

READ< 1 , 1  £ )  IANS 
FORMAT1:  A 1  > 

IF< IANS  .EO.  1 H  >  00  TO  555 
IF< IANS  .HE.  1HY>  GO  TO  550 

WRITE1.'  1  ,  13) 

FORMAT-;  “NEW  GRID*'*  Y/H“  > 

READ-:  i,  12  >  IANSS 
IF* IANSS  .EO.  1 H  >  GO  TO  560 
IF*  IANSS  .HE.  1HY>  GO  TO  565 

CALL  POLAR* 0> 

CALL  LABL* 0 > 

WRITE* 1,2EO 

FORMAT* “ENTER  PEN#  FOR  DAT A . & ” > 

READ*  l,1*)  I  PND 
WRITE*  25, 22  >  I PHD 
FORMAT*  "SPMI  ) 

DO  600  1= 1 , NT I MES+ 1 

PHI=*< 1-1 >*DGREHC+1Se. >*PI/1S0. 

DBPRPT*  I  >=20  .=«  ALOGT*  SCRPE2*  I  >/TFLDMX  > 


PHI  I1**  1-1  >*DGRENC 
WRI T E* 6 , 33  >  DBPRPT* I > , PHI  I 

FORMAT*  IX,  "GAIN  -  POWER  =  "  , E 1 5 . 0 ,  "  dfa"  .  5X,  “ANGLE  =  11 , 
F7 .  3,  X  > 

I F* DBPRPT* I >  . LT ,  -4  0,  >  DBPRPT* I >=-4  0 . 

DUMMY** DBPRPT* I >+4  0.  >/4  0 , 

IX*INT< *XCR-DUMMY=*COS< PHI >=*RDS  >+ 1  000.  > 

I  Y=  I  NT*  *  YCR-DUMM¥=*S  I  N<  PHI  Y+RDS  >*  f  000  .  > 

IF* IANS  .HE.  1H  .AND.  IANS  .HE.  1 HY >  GO  TO  600 
WRITE* 25, 24 >  IX, IV 
FORMAT* " PA “ , I 5 “ , " I 5 * ; PD “ > 

CONTINUE 

WRITE* 25, 26 > 

FORMAT* "PO“ > 


■*  >^11'.*’  ■  >rrP  ■£,+— 


c 


CALL  LABEL-;  1  > 


935 


94  0 
C 


C 

94? 

C 

943 

C 

949 


C 

C 

95 1 
C 


C 


C>0  94  0  I  -  1  ,  NT  i  MES 
F'H1=*  I  -  1  >*D£RENC 
UP  ITE*  6, 935  >  DBPPPT  <  I  > ..  PH  I 

FORMAT*  1  X,  “POWER  GAIN  *  “,£15.3,"  dbn , 5X, “ANGLE  , 
F7.3,"  DEGREE’S" ,  /  > 

CONTINUE 


WRITE*;  £.,70  > 

L'K ITE*:  6,945 >  SFLDMX 

F  OP  MA  T  *  1 X ,  “THE  ECU  wp  W-  MAX  *  “,E15,3,  ,'V  > 

WRITE-:  6, 957  >  TFLDHX 

FORMAT-;  t  X  .•  "THE  TOT  AL  F I  EL  D  H X  —  "  ,  E  1  5 .  3 
WR1TE< e, 7  CO 

1  AVHTM=  I AVNTM/HT I  MES 
WRITE-:  6,947  > 

FORMAT-:  IX,  “THE  AVERAGE  HUMBER  OF  TERMS  IN  fnE  INFINITE  SUM"  > 
WRITE* 6 , 943 >  NT  I  ME  S’ 

FORMAT-;  1 X,  "FOR  "..15,"  NUMBER  OF  ANGLES'  BETWEEN  0  L  360  IS'",  /  » 
WR I TE-:  6 , 94  9  >  I AVNTM 
FORMAT-:  IX,  15,  "  TERMS  .  “  > 

WR I TE-:  6,70  > 

WRITE* 6, 951 >  I MAX 

FORMAT* IX, "AND,  THE  MAX  NUMBER  OF  TERMS  WAS  ",  15, "  TERMS." > 

CALL  EXEC*  23, 5HSMP  ,4,1 SLU > 

END 

END* 


107 


OO  O  O  <J  O  n  'J  O  O  ’T  U'i  iNSKJOO 


V  T=  0  0  ft  C'4  IS.  ON  CROC  07?  US  I  MG  0060?  ELKS  R-OOOi- 
FTW4  .  L 

*ema<  cmblk,  o  :> 

C  «  «:  «  *:*:* 

SU6RCUT I ME  CL SKY< N , M > 

C  +  **:**( •  t 
C 

COMPLEX  A<S,7  >,F<C>, 
t  SUM 


COMMON  /CMBi_k/  A ,  F ,  ECHUPUK  3S  u  0  > ,  DGREnC  .  SCRPEZ*:  j  S  0  0 
1  DEPP PT<3S00> 

DA  Ta  PI  S3.  141 SRCO?/ 

*  CALCULATE  F I  POT  K  0  0  OF  U  P  P  E  R  UNIT  TRIANGULAR  MATRIX 
DO  3  0=2:, M 

A1.  1 ,  J  >=A<  1  ,  J  >XA*;  1,0 

*  CALCULATE  OTHER  ELEMENTS  OF  U  AND  L  MATRICES 

DO  S  1=2,N 
0=1 

DO  5  1 1  =  0,  M 
SUM=<  Cl.  0,0. 0) 

OM 1=0-1 
DO  4  K= 1 , OM 1 
SUM=SUM+A<  1 1 ,  K  >*A>:  K ,  0  > 

a<  1 1 ,  o  :>=«■:  ii,o  >— sum 

IP1=I+1 
DO  7  JO- IP  1,11 
SUM=< 0.0,0.  CO 
I  tit- 1-1 
DO  *  K= 1 , I M 1 
SUM=SUM+A<  I ,  K  >*A<  K ,  00  ) 

A1:  I  ,  oO  >*<  A<  1,00  >-SUM  >/A<  I ,  I  > 

CONTINUE 

*  SOLVE  FOR  F<  I  >  BY  BACK  SUBSTITUTION 

F<  N >=A<  N,N+1 ) 

L=H- ) 

DO  1  C<  NN=  1  ,  L 
SUM— <  0. 0,0.0  > 

I-N-MN 
IP1=I+1 
DO  S  0= I P 1 , H 

9  SUM=SUM+A< I , 0 >*F< 0 > 

10  f<  i  :>=a<  i,M>-sun 

RETURN 

END 

END* 
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I  i -no 004  IS  ON  CRG003-S  USING  000 OS  Cl/S  R=0000 


FTris  .  L 

it.r.n>  Ctit-LK-,  0  ) 

f  ..+  •«  :♦  ■*'  f  t  :*:♦  :*.*  *.*:*•  :«c:|: 

SUBROUT  I  Mt  EC  OPT :  ;:r1AX  > 

C  :t  ♦  1-«:*:»::*:*:***  :f  I. :*:*.*. :* :r.  *  t  :*•*•  «  *  +-  +'  +.  * 

C 

COMPLEX  A*:  b ,  7  > ,  F<  <£.  > 

C 

c 

COMMON  /CMBLK/  A  .P,ECHWPW*I  30  00  >,  DGRENC,  i.CF:P£2< 
I  OBPRPT<  3c-  0  0  > 

C 

C 

PaTA  P 1/3 .  1  4  i  5'3*2©5/ 

NX OP =2  ,  *  1  000 
NY  Or- i .75*1 000 
C 

NXSTP= 0 . S" I  0  0  0 
NYSTP* 0 . 5^ t  000 
C 

N.vSTPS=* 

NYSTPS= 1 0 
C 

I PEN® I 

WR I T£< 25 ,  3  >  I PEN 
C 

NX£P=NXOP+NXSTPS»NX:STP 
WRITE*'  £5,  1  >  NX  OP,  NY  OP 
WR I TE< 25,1  >  NXEP , NY  OP 
WRITE*'  25,2  ) 

C 

N  YEP=NY  OP+N  YST  PS:*NYSTP 
WR I  TE>:  25 ,  1  >  K.Y  OP ,  NY  OP 
WRITE*;  25,  t  >  NX0P,MYEP 
WRI TE< 25, 2  > 

C 

DO  100  0=  1  ,  NX'STPS 
I  X=NX  0P+  J*NXSTP 
- 1 Y 1 = MY OP-5 0 
IY2=NY0P+50 
WRITE*: 25,  1  >  IX,  I V 1 
WRITE*: 25,  1  ')  IX,  IY2 
WRITE*:  25, 2  > 

100  CONTINUE 
C 

DO  200  1=1, NYSTPS 
I  Y= NY  0P+ 1  *N YSTP 
I X 1 ® NX OP-5  0 
1X2= NX OP +50 
WRITE*;  25,  1  .■)  1X1  .1 Y 
WRITE*:  25,  1  >  1X2,  IY 
WRITE*:  25, 2  > 

200  CONTINUE 
C: 

CALL  eclbl 
CALL  LAPEL*:  1  > 

c 

NPT  S=  IF  IX*;'  ISO.  /DGRENC  > 


T 


3  c- 00 


1 09 


c 

!PE:-r=3 

wk  I TE*:  25 , 3  I PEN 
C 

DO  400  IPhI  =  1  ,NPTS 

i  X=  I F I  .'O;  <  I  PK I  -  i  .  >  *  DC-F  3 NC/2  0 .  * NX3TP 
I  I F I XC  \  EvHfi'ViiK  Ir'ril  .1/ 0  i  1  * Mv" 

WRITE*: 25/  1  >IX,  IV 
400  CONTINUE 

C- 

WRJ TE<25,2> 

1  FORMAT*  “  PA “  ,  i  5 “  ,  “  1 5 ,  “  ;  PD  "  > 

2  FORMAT < "PU” > 

3  FORMAT*;  "IN;  SP ",  1 1  ) 

END 
END* 


no 


■4 1'  0  o  o 
cTP-f  1 


I- 


«  a  cn  r>  o  o  o  o  o 


L  T=00uu3  IS  OH  CROC' OSS  OS  I  HO  0001?.  BlKS.  ft-  0000 
FTN4  .  L 

SUSROUT  I  HE  EES  EL*;  N ,  X ,  S  .  V ,  DP ,  VP  > 

COMPUTE  THE  BESSEL  FUNCTION  OF  ORDER  N,  N  AN  INTEGER  H>OR«C- 
WITH  REAL  ARGUMENT 

ALL  EGUrtT  I  ON  REFERENCES  TO  ABftAMOUl  TIT  A ND  3TE0UN 

D I  MENS  I  ON  C  On  7  > ,  C 1  <  7  > ,  D  0<  7  >  .  D 1  •;  7  > ,  E  0<  7  >  .  E  K  7  > ,  0  iX  ?  >,CK  ?  > 
DAT  A  CO/1  .  0,  -0' .  24**S*7,  t  . 26362 OS, -. 3 1 63*66,  .  44447SE-1  , 

*-  .  3*->44E-2,  .  2  i  E-3/ 

DATA  Cl/. 5, -.5624**85,  .  21  0*3573,  - .  30542SSE- i  ,  .  44  331 SE-2, 

*  -  .  3 1 76- 1  E-3,  .  1  1  OSE— 4/ 

DATA  DO/ .  7*788456. 77E-6, 55274E-2, -.*5l  2E-4  ,  .  137237E-2, 

—  . 72605E— 03 r  »  1  4476E-03/ 

DATA  D I/. 7*7004 50,  . 156E-0,  . 165*667£-1 ,  ,  171 05E-3, - . 24*5 1  IE-2 
«■ .  1  I  3653E-02,  -  .  20033E— 03/ 

DATA  E  0/-  .  705390 1  6, -.  4  1  663S7E- 1  ,-,3*546-4,  .  2625  TOE -2 ,  -  .  54  ;  0 
■*- ,  2*33 3E- ©3,  .  13S58E-03/ 

DATA  £1/-2.356r*44*,  .  124**612,  .565E-4,-.63  7878E-2,  .7434SE-3 

*  .  7*024 E-  03,  -  .29  t  66E-03/ 

DATa  GO/ . 3674 SO*, . 6035*37, 7435030, . 23300 1 2, - . 426 1 2 1 E- 0 1 , 

■■i  ,  427*1  6t-02,  -  .  24S46E-03/ 

DATA  01 /- ,  6366  1  *0,  .  221 20*1 , 2 . 1  68271  ,  -•  1  .316483.  .31-3*5, 

*- . 400S76E-01 , . 27S73E— 02/ 

DATA  P I /3 . 1 4 1 5*26/ 

IFLG-0 

IF  <H.LT. 0>  I  FLO- I 
N=  I  ABS<  N  > 

IF*;  ADS*;  X  > .  LT  .  1  .  OE- 1  0  >60  TO  1  5 0 
I  F<  ABS<  X  )  .  GT  .  3 . 0  >G0  T  0  5  0 
X3SCi=X*X/S .  0 
PROD* I . 0 
B 0=0.0 
B 1=0.0 
CUM 0=  0 . 0 
CUM 1=0.0 

SEE  EQUATIONS  *.41  AND  *.44 
DO  5  1=1,7 
B 0=B 0+C 0<  I  >*PROD 
B 1  =B  1  +C1  *;  I  XkPROD 
CU11 0=CUM  0+0  0*;  I  >*PRGD 
CUM  1  =CUM  1  +0 1  *;  I  >*PRGD 
PR0D“PR0D:*X3SG 
CONTINUE 
B1=B1*X 

XC=2  .  0*8NGL<  DLOG<  DBLEn'  0 . 5*X  )  >  >/PI 
Y  0=XC*B  0+CUM  0 
Yf=XC*£:l-*-CUiTl/X 
GO  TO  1 00 

EOS  *.4.3  AND  *.4.6 
0  THRGVX=3.0/X 
PROD— 1 . 0 
F  0=  0 . 0 
FI =0.0 
THFTA0=X 
THETA  I =X 
DO  35  1=1,7 
F0=F0+D0*;  I  XfPftOD 
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T 


7 


t  00 
1  U  f 


i  o; 


c 

1  i  o 
c 


lie. 


117 
i  is 


c 

130 


135 


1*10 


r  1  =  I  1  +l>  1 1  I  PC' 

THET  A0=THE1  aO-*-ECK  J  ii'PkOD 
THETA1=THE"Ihi  +EIC  1  VvFEOD 
PR  rR 00  ♦'  T HRO  VX 

CONTINUE 

SGRX*  )  .  0/SNGL1;  DSORTC  DBLE<  X >  > 

B 0=SQRX*F 0-fSNGL*:  DCOSC  DBLEC  THcThO  >  >  > 

&  1  -  SCfK.'vi'r  1  ;*SNGLC  DCCB<D£lE<  THETA 1  >>  > 

YO=SORX*F  G*SNGL<  DSIN<DBLEt  THETAO  >>) 

V 1  *S£>RX*F  1  •«SN&'L< DS I N< DBLE<  THETA 1  >>  > 

IFcN-i  >1  01 , 1  05,  1  1  0 

£:=B0 

BR=-6i 

Y*Y0 

VP*— Y 1 

CO  TO  2 0 0 

B=Ri 

RP=B0-B1/X 

V-V1 

YP*Y0-Y1/X 
GO  TO  £00 

FOR  RECURRSIVE  DIRECTION  COMMENTS  SEE  SECTION  S.12.R3S5 
XN-H 

1F< XH . L  T . ABS< X  >  >GO  TO  130 
FOR  i«  H  RECOR  DOWNWARD 
BLAST- 1 . 0 
BLASTP--0.  0 
J=N+ 1  0 

DO  115  1=1, J 
XI  =  v'-I 

BNEXT=2 ,  0*XI*BLAST/X-BLASTP 

E'.LABTP=BLAST 

BLAST  =BNEXT 

JF< I .NE. 1 0)00  TO  115 

6LF-BLASTF 

2-6HEXT 

CONTINUE 

I F>;  ABS<  B  0  >  .  LT  .  ABS<  B 1  >  >GO  TO  117 

CORR = B 0/B L  A ST P 

GO  TO  IIS 

CORR=-BI/BNEXT 

B=BLP!*CORR 

&NMIM 1 =Z*CORR 

GO  TO  140 

FOR  X >N  RECUR  UPWARD 

BLASTP=B0 

BLAST=B) 

DO  135  1*2, N 
XI=  r-i 

BHEXT*2 . 0=*X  I  •■♦BLAST /X— BLAST P 

BLASTP=BLAST 

BLAST*BHEXT 

CONT I HUE 

B= BLAST 

BNMIN I =BLASTP 

BP*BNMIH1-XH*B/X 

YLAS'I  F-YO 

YLAST-V1 

DO  145  1*2, N 

X I  —  I  —  1 
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OOOOOOOOOOOOOOOOOOOOOOOOOOOGOGOOG  0-0  o  o 


o  o  r>  r>  o  o  o  o  -o  o  o  o  o  o  o 


1M  00dU4  IS  ON  CROC 0 USING'  00020  ELKS  Fr=00iH* 


FTN4  .  L 

*Eriw<  ELK MM ,  0  :> 

c  « •+  -i  ■*'  i’  *  *  •'*  *  ST :+ b  ■  V  :f  :*  :*•  :*  :f  :* 

PkOGk'hM  kwOpiM 

c  ■*  ;f:  4*  4*  :f.  r  •*  •*  :*  :f  :*  *  f  :+  :*  *;*.:*:*  .* 

C 

C  ROBERT  K .  SCHNEIDER 

C 

C  *  This  program  uses-  the  MOMENT  hETmOD  a*  presented  bv  oh C-h  RICHMOND 
C  <  “Scattering  by  a  Dielectric  Cylinder  of  >irfcit nary  Cross  Section 

C  Shape’1..  IEEE  HRS..  Nov  i?S-»,  p .  334  >  to  determine  the  Sc  at  verso 

C  Field  from  and  Echo  Uidth  ot  a  dielectr  ic  cylindrical  shell  of 
C  circular  cross  section. 

C 

■t  The  dimensions  of  the  scatterer  are  read  from  a  data  file  -  DhThFi 

INTEGER  NCELLS , MPTS , PTOBS, IRHI , MCEcLS 

REAL  XNN,YNN,EMAGI,  Jt  ,  JG,K(»,K2,YG,YI  .EmHGN, 

1  ECHitlRU 

COMPLEX  CMH.EINC1 , EIMC0, EINC2, 

1  ESCAT , 

2  THU  1 , NLMDH , ALPHA , COEF , FAC , C 1 . C2 ,  V,V1 , V2 , 

3  THU , V I N , VOUT , ONE , ZERO , TLFLD , H  0 , H 1 

DIMENSION  I B UF  FX  1 *  > 

COMMON  /BLKMM/  HA, B, C, R, FRED, PERM, DGRENC, EPSLNR, AN, NCELLS, 

1  EINC1<360>,  VIN<  1851  >,  ESCAT<  3©0  >, 

2  ECHUPUK  36  0  > ,  K  0 ,  K2 ,  VOUT  i  S5 1  > ,  T  HU'.'  1  S' 5 1  > , 

3  *  E INC'2,  SCRPEZt  3&  J  ),Dr*F'KF  I'i  dfi>i  1,  wu1.  !c**5t  i,HJ*.  tool 

EQUIVALENCE  < MCE LLS,NZ > 


DHTH  PI/3 .  1 4  1  5?26‘5/ 

DhT H  XCft/4 . 375/ , YCft/3 . 875/ , RDS/2 . 743/ 

DHTH  IBUFR/2* 0 , 2HDA,  2HTH ,  2HFL  ,3*0, 22 1  B ,  7* 0/ 


*  SET  SPOOL  FOR  INPUT  FROM  IBUFR  THROUGH  ISLU 


CALL  SPOPi-K  IBUFR,  ISlU> 

CALL  EXEC':  22 ,  1  > 

READ-.  ISLU,  *  J  HH,B,  C,  CURENT  ,  R, FRED ,  PERM ,  DGRENC  ,  EFSlNK  ,  NCELLS , 

t  HLFCEL 

*  DETERMINE  THE  CELL  STRUCTURE  AND  «  OF  CELLS 


C*  *** 
C 

C=*T-t  * 


U=  (i .  0 

CALL  CELL'.'  U  > 

I F<  HLFCEL  .NE.  t>  GO  TO  5C« 
MCELLS*HCELLS/2 


c  *« 

c 

030 

C 

c  * 
c 

50 

C 

C 


12 

c 

C  * 

c 

21  Ci 
C 

C:  ♦ 

c 

220 

C 

C:  * 
C 

C 

C 


260 

C 

C 

24 


C 

23  Ci 


C 

24  0 
C 

250 

26' 


DETERMINE  THE  COORDINAiES  OF  THE  CENTER'S  OF  £hCh  SEl 
CM LL  CLCF.DC  W  > 

DETERMINE  THE  ECHO  WIDTH  FEE'  WAVELENGTH 
XMAX* 0 . 

CALL  ECHOW<XMAX,W  » 

NPTS- 1 FIX< 36  0 . /DGRENC > 

WRI TEC  1 , 23  > 

FORMATC  "PLOT  THE  SCATTERED  FIELD'' £  Y/N"  > 

READ*;  1,12>  I  AM'S: 

FORMAT*;  A 1  > 

IFC  IANS:  .EC.'.  1 H  >  GO  TO  210 
IFC  IANS  .ME.  )H  Y  >  GO  TO  220 

PLOT  ECHO  WIDTH  PER  WaVELEMGHT 

CALL  ECOPT C  XilAX  i 

DETERMINE  THE  FAR  ZONE  'SCATTERED  FIELD  AT  £*CH  DGREM 
CALL  FLDCT 

DETERMINE  THE  INCIDENT  FIELD  IN  THE  FaR  ZONE 

CALL  FLDHC'  1 , W > 

TFLDMX® 0 . 0 

DO  26  0  1  =  1,  NET'S 

TLFLD=E1HC1C  I  >+ESCATC  I  > 

SCRPE2C  I  >=C  ABS:*;'  TLFLD  > 

*IF<  SCRPEZ*.'  I  >  .  GT  .  TFLDMX  >  TFLDMX=SCRPEZC  I  > 
CONTINUE 

SCP.PE2CNPTS+1  >= SCRPEZ*;  1  > 

WRITEC 1 ,24  > 

FORM AT< -PLOT  THE  TOTAL  FIELD V/N“> 

RE ADC t  ,  I  2  >  IANS 
IFC IANS  .EO.  1 H  )  GO  TO  230 
IFC IANS  .ME.  1HY)  GO  TO  270 

WRI TEC  1 , 23  > 

FORMAT C “MEW  GRID 76  YXM " > 

READC 1 , 1 2 >  IANS 
IFC  IANS  .EO.  1 H  :>  GO  TO  24  0 
IFC IANS  .HE,  1 HY >  GO  TO  250 

CALL  POL ARC  0  > 

CALL  LABLC  0 

WRI  TEC  1  ,  2S  ') 

FORMAT*.' "ENTER  REM'S  FOR  DATA. 6  "> 

READ*;  \  ,*)  IF  NO  - 
WRITEC 25, 22 >  I PHD 
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format-;  "GP" ,  1 1  > 


7  (i  (i 


31 

C 

2?  0 
26 


C 

SOO 
S3  0 
C 


36  0 

370 

C 

27 

C 


300 
4  0  0 
C 

300 


DO  700  I  =  i  , HPTS+ 1 
•  PH  I  =*  <  1-1  >s*OCRENC'+ 1  v  0  .  >*p  I ISO. 

DSPPPT <  I  >=20. +ALOGT'.  SC-RPEE*  I  >.'TF Lvri:<  > 
IP< DE-FEET*  I  >  .LT.  -40,  >  uEFEPTr  I  0 
DUMMY**  DBPEPT *  I  >+4  0.  >.-'4  0. 

IX* I NT* < >iCR-DUMMV»OOS< PHI  >*EDS  >«  i  000.  > 
I Y*  I  NT*  *  VCF:-DUnMY!«S  IN*  PH  I  >=*EDS  >*  t  000.  > 
WRITE* 23, 23 >  IX, IV 
FuRMAT  *  “  PA  “ ,  13“,  "13";  PD  "  > 

CONTINUE 

WP  I  TE*  23 ,  3  i  > 

FORMAT'.  "PU“  > 


WRITE* f , 2c > 

FORMAT* "HARD  COPY  OF  CALCULATIONS?# 
RE AD* 1, 12 >  I AMS 
IF*  IANS'  .EC.  1H  >  CO  TO  SOO 
IF* IANS  .HE.  1  Hr >  CO  TO  5  00 


V.y 


WRITE*  S, S50  >  TFLDMX 
FORMAT*  IX,  “TOTAL.  FIELD-MAX  =  “,E15.S,^ 

DO  S70  1  =  1,  NPTS 
PHI** 1-1 >*DGRENC 

WR ITE*  6 ,  SS 0  >  P.H  1  ,  E CHWPW*  I  > ,  DBPRPT *  I  > 
FORMAT* IX, “ANCLE*  ",F5.2,5X, “ECHWPW* 
"DB  GAIN*  “ , E 1 5 . S , / > 

CONTINUE 


WRITE* 6, 27 > 
FORMAT *  1 X , //X  > 


DO  400  1*1, HPT S 
*  PHI**  1-1  )**[>G'RENC 

XMAG=CABS*  *  ESC  AT  *  I  >  >  > 
wr<  I TE*  6 , 3  0  0  >  XMAG ,  PH  I 
FORMAT* 1 X, "MAGNITUDE  ESCAT = " , 2X, 1 E 1 5 . 
CONTINUE 


CALL  EXEC*23,5HSMP  ,4,ISLU> 
END 
ENDS 


no  o  o  o  oooooooooooo  o  o  o  o  o  o  o  o o o  o o o 


T=  0  0  00-1  IS  ON  CkOC'OS?  USING  0  0  OCA  c Li'.S  RsOOCC 
FTNA.l 

*EMAv  BlKMM,  0  > 

C  **■<*  +  i:*:**:**:*:*:**  f  *  *  4>  *:«  ■»  *■.,.:*:*:*:*»  f**-, 

SUBRCUT  I NE  CELL*.'  W  ;> 

t  *  f  1  f  :♦  :f  ■■*:*:*  ♦  M  :*  :*  >  t  .*  <  ■*  *f-«  :*  *  :fc 

*  THIS  SUBROUTINE  DETERMINES  THE  #  OF  CELLS  NEEDED  PND  THE 
DIMENSION  OF  A  CELL  TO  SFhN  THE  OBSTACLE. 


INTEGER  NCELLS, NETS, r  TOC'S, I PH I , NCELLS 

REAL  XNN,  YNH,ErtAGI ,  jl  ,  JO,  KO  .  K2,  YO,  VI  , EmAGN , 

1  ECHWPW 

COMPLEX  CNN . E I NC 1 , E I NC 0 , E I NC2 , 

1  ESCAT , 

2  T AU  1  ALMDa  .  AL  PHa  ,  C  OEF ,  FAC  ,  C 1  ,  C2 ,  V ,  V !  ,  V2 , 

3  ■  TAU,VIN,VOUT, ONE, ZERO, A0,A1 

D I  MENS  I  ON  I BUFR*;'  1 t-  > 

COMMON  /’BlXMM/  HA , B,C, R, FREQ, PERM , OCRENC, EPSLNE, AN, NCElLS , 

1  EINC'1  <  3S0  >,  VIN<  1051  ),  £SCAT<  3S0  >, 

2  ECHWPv*.  36  0  > ,  K  0  ,  K2 ,  VOUTC  1 S 5 1  > ,  T A'.k  1  S"?  i  > , 

3  E I NC2 ,  $CRFEZ<  56 1  >,  DBPEPT.  So  l  >,A0*.  tS5l  ">,  Ai*.  lS3t  > 

EQUIVALENCE  < NCELLS, N2 > 


DATA  P 1/3, 1415S265/ 

*  THE  -£  OF  CELLS  FOR  THE  CIRCULAR  DIELECTRIC  SHElL  SChTTERER  IS  a 

PARAMETER  WHICH  IS  ASSUMED  AT  THE  INITIATION  OF  THE  PROGRAM  mND 
IS  READ  FROM  THE  DATA  FILE 

*THE  WIDTH  OF  EACH  CELL  MUST  BE  LESS  THAN  OR  EQUAL  TO 
0.2  *  WAVELENGTH  /  SORT  <  EPSLNR  > 

*  SINCE  THE  NUMBER  OF  CELLS  IS  A  KNOWN  VALUE  AND  THE  SI 2 E  OF  THE 

STRUCTURE  IS  DEFINED,  THE  WIDTH  OF  A  SQUARE  CELL',  W,  IS  EASILY 
DETERMINED  TO  BE 

W-2 . *PI*C/NCELLS 

*  RADIUS'  OF  CIRCULAR  CELL  WITH  EQUAL  AREA  AS  'S'CUhRE  CELL’,  AN,  IS 

RGUT*C 
R I N=C— W 

AN*S6RT<  <  ROUT**** 2 .  -RIN+*2 .  VN CELLS  > 

200  WRITE*:  6,  1  00)  W,  AN, NCELLS 

1 00  FORMAT  < IX, "W=“,2X, 1E15.9,5X, “AH= " , 2X , 1 E 1 5 . 9 , OX , “NCELLS*" , 

1  2X ,14, // > 

C 

RETURN 

END 

END* 
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O  O  o  o  o  moo  o  o  o  *o  o  o 


US: IMG  00 Oii-4  DLLS  P.'=  0 0 *j 0 


D  7  =  00004  i'_.  UN  CROOuSS* 


FTn-4  ,  L 

*ENA<  BLKMrt,  0  > 

SU&RGU7 1  NS  •  CLCRP*;  US  i 


*  THIS  SUBROUTINE  CALCULATES-  THE  COORDINATES  OF  THE  CENTER  OF 

each  cell. 


INTEGER  WELLS,  NETS,  PTOBS,  I  PH  I,  MCELLS 

REAL  XNN ,  VNN ,  EMaGI ,  J1  ,  J 0 ,  K 0 ,  K2 ,  Y 0 ,  V  I  .SfViGN, 

ECHWPW 

COMPLEX  CNN . E I NC 1 , E I HC 0 , E I NOS , 

ESCAT , 

T AU  i  ,  ALMDA ,  «L P Ha  ,  C  OEF  ,  FaC  ,  C  1  ,  C  S  . V ,  V 1  ,  VS  , 

T  AU , V I N , VOUT , ONE , ZERO , A  0 , A i 

DIMENSION  IBUFR< 16 > 

COMMON  XBL  KNM/  A  A ,  S ,  C ,  R ,  FRE  0 ,  PERM ,  DGRENC  ,  ER  SlNk  .AN,  NOEL  L  S- 
EINC-1  •:  30 0  >,  7IN<  37  OS  >,  ESCaT<  3c- 0  >,  >  NN-.  37  OS  > , 
VNN  .  3702  >,  EC'HWPW*;  36  0  .  n  0,  K2,  V007<  37 02  >,  TAU< 
EINC2,  SCRPEZ*;  36 i  >,  DBPPPT*:  36 1  > ,  AO*.  37  02  ,  a  1  <;  3 

EQUIVALENCE  < NCELLS, M2 > 


DATA  P 1/3. 14152265/ 

*  RADIUS.  OUT  TO  CENTER  OF  CELL  n 

RN=C-< LI/2, > 

*f  INCREMENT  ANGLE  FOR  CELL  LOCATIONS 
THET  AN=2 .  «p  I /NCELLS 
WRITE*: 6,  5 CO  T.HETAN 

0  FORMAT*:  1 X,  “TH£TAN=  " ,  2X ,  1  E 1  5  , '? ,  //  > 

*  DETERMINE  THE  CORD I NATES  OF  THE  CENTER  OF  EACH  CELL  n 

DO  10  N= 1, NCELLS 

theta=*:n-(  :>*thetah 

XNN<  N  >»RN*COS<  THETa  ) 

YNH<  H  >=RN=»S  I  H<  THET  A  > 

WR ITE*:  6,  1  0  0  >  XNN*.  N  > ,  YNN<  N  > ,  N 

0  0  FORMAT*:  1 X ,  “  XNN=  "  ,  2X ,  I E 1 5 .  ? ,  5X ,  ”  YNN=  "  ,  2/ ,  I E 1 5 .  S ,  5X ,  “  N=  “  ,  2X 

0  CONTINUE 


RETURN 


Ei 


OOOOUOOOO  O  O  1 J  O  O  O  O  O  -ii  O  O  O  O  O  O  O  O  O  O  If/  o  o  o  o 


!  • 
i 
> 

,JpW«. 


W  T  =  0  C  0  04  IS  ON  CPC-OCC?  USING  SI- DOS  £:LLS  Fr- Ouv  :• 

FTN4 , 1_ 

•*£MA.;  EL’-'nri ,  0  > 

C  •»•« •***•*•»*!:». *****.*•«:«.:* -•!*:<  +  **■**.»..«.:+.*;*:* 

SUBPvLlT  1 NE  ECHO'.!*:  XMaX ,  U  > 

-«  THIS  SUBROUTINE  CALCULATES  THE  ECHO  iJIDTH  PEE  ONI T  NAVEL EMGM7 , 
RELATIVE  TO  THE  INCIDENT  F  I  ELD  AT  THE  CENTER  OP  ThE  OBSTACLE.. 
FROM  A  PI  ELECTRIC  CYLINDRICAL  SHELL  OF  CIRCULAR  CROSS  SECTION’ 

in  the  presence  of  a  radiating  currnt  filament  . 


INTEGER  NCELLS, NPTS, PTOBS,  IFHI ,  MCELLS 

REAL  XHN ,  VMM  .  ENACT  ,  v*t  ,  JO,  K  0 ,  K2 ,  Yu  ,  Y 1  ,  SnACN, 

1  ECH'-lPUi 

COMPLEX  CHN • El NCI , El MC 0 , E I NC -Z . 

1  ESCAT , 

2  TkUJ  ..  ALMDA  ,  ALPHA,  OGEF  ,  FAC..  C  i  .C£,V,Vi  ,V2, 

3  ThU ,  VIM ,  VOUT ..  0H£  ,  EERO,  DUMMY4 ,  r. 0,  At 

D I MENS I ON  I RUFRC t  6 > 

COMMON  XBLK'MM/  AA,  B,  C,  R,  FREQ,  PERM, DCRENC,  Ef SlNR  .  AN,  NCELLS, 

1  El HC 1  <  36 0  > ,  V I N<  1851  > ,  ESChT 36 0  > , 

2  ECHWPWC 36 0 > , K 0 , L 2 , V©UT< 1 85 1  > , Talk  1 85 1  > , 

3  E I NC2 ,  SCRPE2<  36 1  > ,  DBPRPT*:  36 1  )  ,  «  CK  1  85 1  > ,  A 1  <  1  85 1 

ECU I  VALENCE  < NCELLS , H2  > 


PAJA  P 1/3, 14 158265/ 

•■*  CALL  FLDTL  TO  DETERMINE  THE  TOTAL  FIELD  IN  THE  OBSTACLE 
DUE  TO  THE  INCIDENT  FIELD  UPON  THE  OBSTACLE 

CALL  FLDTL<  U > 

*  CALL  FLDNC  TO  DETERMINE  THE  INCIDENT  FIELD  AT  THE  CENTER  OF 
THE  OBSTACLE. 

CALL  FLDNCC 2 , W  > 

NPTS* I F I X< 36  0 . XDGRENC > 


DO  ?  1=1, NPTS 

NK I TEC 6,5)  El NC2 

FORMAT-:  1  X ,  “ECHOW  :  E I  NC2=  "  ,  2X ,  1  E  1  5 . 8 , 2X ,  E 1 5 .  S ,  X  > 
7  CONTINUE 

*  DETERMINE  THE  MAGNITUDE  SQUARED  OF  THE  INCIDENT  FIELD 


C 

c 


EMAG I  =CABS<  ■:  E 1 NC2  >  )fCABS-;  <  E I NC2  >  > 
DUMMY  I  =K  0*P  I  **2 .  *FREG\K  3  .  E8=*EMAG  I  > 
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L0‘  Ip  NlELLS  F  UK’  THE  11  0  CmTT  ER  i_D  r  iEl.O*'  If')  r He.  FhP 


C 

c 

c 


t  Cl 

c 

c 

cc 

cc 

c 

c 

c 

f  00 

c 

20 

c 


DO  20  IPhI=i,NPVS 

PHI=i;  iPHl-i  .  >*DGREr;C*FI.’,1S:Ci. 

OUnMY4=< 0 . ,0.  > 

rn=c— :  w/2.  :> 

THET AN=2  .  *P  I  /HCELLS 

Du  i  0  N*  1 ,  NCELLS 

THETA=<  N- 1  THE ThI* 

XNN=Rf-*'CC'S<  THETA  .» 

YNN«RN*SIN<  THETh  > 

ARC  J  t  *K  0*AM 

CALL  E’ESEL  >■  1 ARC  J I  EE V  /  BP,  Yr ) 

J 1  =BB 

ARGH*K  ©*<  KNN ■fCOS1,'  F'H  i  )+VHH*  SIN<  PH  I  >  > 

DUmMY2=C0$< ARCH ) 

DUMM  Y3=0  i  w  ARGH  > 

DuHMY4=DUririY4  EF’SLNR—  i  .  >*V0UT<  N  >««M*  J  l  •♦CMPL;'1'  DUhM 
CONTINUE 


EMAGN=CAF:£‘.  0UNMY4  >*CAE"3<  DUMHY4  > 

ECmiPU<  I  PH  I  >=DUrtf1Y  I  *EMAGN 
IF-.  ECHWPLK  IPHI  >  .GT.  XM«X>  XMAX-ECHWPWs  I  PH  I  > 


PHI  I -PH  I  ♦  ISC'. /PI 


UR I T£<  6 , 1 0  CO  ECHUPUK  I PH I > , PH 1 1 
FORMAT  ■;  IX,  “ECHwF’U-  " ,  2X,  »EJ5.9,5X,  "PHI=" , ZX,F? .3 , / > 

CONTINUE 

RETURN 

END 

END* 
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L  T-C  0004  IS  ON  CkCOOZS  USi NC  6  00  OS'  ELKS  k-L  'JijO 
FTH-I ,  L 

*EMA<  El. KM M,  0  > 

c  -ft  ft  :ft  >  4*  *■  -f  :ft  :ft  .ft  :ft  :ft  :ft  :f  t  .ft  :ft  :ft  :ft  .ft  ;ft  :ft:ft  -I.  .ft . ,.  ft 

SUBPOUT I  HE  FLDTLC  v  > 

C  Mftft  ft*  .ft  ft.  ft:ft  :*:*.:  ft:ftift:  ft  :ft:ft:ft:ft:ft  :ft:ft:ft:ft:ft:ft 

c 

C  *  THIS  SUBROUTINE  CALCULATES  THE  TOTAL  FIELD  In  THE  OBSTACLE 


INTEGER  NCELLS, NPTS, RTOBS, IPHI.MCELLS 

REAL  XHM , YNN , EHAG i  ,  Ol  ,  sK>, KO, K2,Y0 ,  V 1  ,ErV:C-H, 

ECHUPW 

COMPLEX  CMN , E I NO I , E I NC 0, E 1 NC2 , 

ESC  AT , 

TAU  i  ,  ALriDA ,  ALPHA ,  COEF ,  FAC  ,  C  t  ,  C2 ,  V ,  • '  i  ,  v2  , 

T AU , V I N , VOUT , ONE , ZERO , A  0 ,  « 1 

D I  HENS  I  ON  I  BUFFS  t  A  > 

COMMON  /BLKMM/  AA, B,C,  P ,  FFEO,  PERM, DGPENO,  EPSlmR,  AN,  NCELLS, 

El  NCI  <  36  0  > ,  V I  NC  1  St*  i  > .  £  S  CA  T  <  3  6  0  > , 

ECHUPUK  36ij  >,K0,K2,  VCUT*;  ('••■'1  >,TAU<  IS?)  >, 

E I NC2 ,  SCRPEZC  36  i  >,  DPP  APT <  36  1  > ,  A  Os  1  85  i  >,  A  i  ■:  i  S5  i  > 

EQUIVALENCE  <  NCELLS, H2 > 


DATA  PI/3,  14150265/ 

C 

C  *  OBTAIN  THE  FIELD  INCIDENT  UPON  THE  OBSTACLE 

C 

CALL  FLDNC<0,W> 

C 

C  DO  5  1=1, NCELLS 

C  UR ITE< 6,7  )  El NC <K I > , I 

?  FORMAT <  I X ,  "FLDTL  :  E I  NC:  0=  "  ,  2X ,  1 E 1 5 . 3 ,  2X ,  E 1  S  .  8 , 5X ,  “  I  =  “  ,  2X , 

1  14,/) 

C5  CONTINUE 

C 

RN*C-< U/2 . ) 

THE TAN=2 . *P I /NCELLS 
C: 

C  DO  30  M=l, NCELLS 

CC 

(1=  1 

THETAM-;  M- 1 >* THE TAN 
XNM=RN*COSx  THET  AM  > 

YNM=RN«-SIN< THETAM  > 


DO  20  N= 1 , NCELLS 

•«  DETERMINE  THE  WEIGHT  Can  ON  Er, 


TH£TA=C  N- 1 >»THfc7AN 


r*j  r>  0  0  0^)0  DO  o  o  o  o  o 


I 


I 


< 

I 


X)*N«{(l '  »COi>>:  THE  i  h  :> 

S I  i><  T  HE  T  a  ’ 

I F*.  N  .  N£  .  !  GO  T  Ci  1  0 

DUMMY 1=K0*hH 

CALL  BESHlC  1  ,  DUrtnY  1  ,BB,  ¥ ,  BP ,  VP  > 

J  i  =BB 
¥1  =Y 

WRITE-:  6, 0>  J» ,  ¥t  ,65 

FORMAT-:  1  X ,  “ FLDTL  :  01=",  2:< .  1  E 1 5  .  S  ,  5X ,  "  V  t  =  •  ,  2>{ ,  1  £  1 5  .  $  ,  5/ , 
“  BR=  “  ,  2X ,  i  E :»?  .  8 ,  /  > 

CrtN»l  .  «C EPSLNR- i .  >:*-<  *,  P I  :»K0:»»-.n/2  .  '>*••;  CrsPLXn  t  1  ,o1  >>♦  i  .  > 

AC  M ,  N  >=C  h.M 

I FC  M  .  ME  .  t  )  GO  T 0  £  0 
T  AIK  N  >=Chn 
GO  TO  £0 

DUMMY  1  =P  I -+K  0*»i  1/2 . 

RMN=SGRTC  C  XNM-XNN  >*<  XNM-XNN  >  + 

C  YNM-YNM  .»-*<  VMM- VMM  >  j 
AkGH=K  C'-RriN 
ARG  J=KO:-rtN 

CALL  BE'SELC  I  ,  AF.G  J,  55,  V,  BP,  VP  > 

J 1  =B8 

CALL  BESEL C  0 , ARCH , BB , V , BP , VP > 

Jd=BB 
Y  0=  Y 

CMN=DUMMV  1  *  C  EPS LHP- 1  .  >* 0 1  '»C!'1PLX<  V 0 ,  JO  > 

IFCN  ,  NE.  1 >  GO  TO  20 
TALK  N  >=CMm 

ACM,M  j=CMM 
Ci  CONTINUE 

AC  M ,  N+  f  >= -E I  AC  OC  M  ) 

3  0  CONT I  HUE 

*  NOW  THAT  THE  MATRIX  HAS  BEEN  FORMED,  SOLVE  FOP  £n 

MCELLS-NCELLSt  1 

WP I TEC  6 , £5  0 >  MCELL S , NCELLS 

50  FORMAT-:  1 X,  "MCELLS*  "  ,2'A,  14, 5X,  “NCELLS—  “  ,  2X,  14,// > 

DO  400  M= 1 , NCELLS 


5  0 

1 

C300 

C 

C4  00 

C 

C: 

c 

7  0  0 
C600 

c 


DO  300  M=1, MCELLS 

WRITEC  6,350  >  AC  ri ,  N  >,  M,  N 

FORMAT*;  IX,  "FLDTL  :  A=“,2X,  1 E  J  5. 9, 2X,  El  5 .  S',  51;,  “M=  “  ,  2X,  J  4 , 
5X, “H=“ , 2X , I 4 , / > 

CONTINUE 
WRITEC f, 3730 
FORMAT-:  TX,  ///  > 

CONTINUE 


DO  600  1=1, NCELLS 

WRITEC 6, 700 >  TALK  I),  I 

FORMAT-.'  1 X ,  ”  TaU=  “ , 2/ , 1 E 1 5 . 9 ,  2X ,  E  1  5 . 9 ,  SX ,  1 4 ,  /  > 
CONT I HUE 


U 
1  0 


C 

C 
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j 


iRi.2  1 , i-;r‘ OEM ,  ILft> 


URII&K  I  . 5C0>  I£R 
500  FORMAT*:  IK,  "  I EE="  14  > 

C: 

*  En  ARE  COMTtt Ji-ED  IN  THE  ARRAS'  vOUT 


c 

DO  100  1=1,  NO! 

c 

WRITE*;  6,200  > 

200 

FORMAT*;  IK,  “F Li 

1 

2K,  14,. 

C  i  0  0 

CONTINUE 

C 

RETURN 

END 

END* 

DTl 


VOUT =- 


* 
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0  T-OOOi:a  IS  ON  CROOO'JB  USING  000  OS  ELK'S  ReOtuC 
FTN4 ,  L 

I  EM  AC  Bi  MriM ..  A  > 

C  Jf-’f  -t  :f  :r-i  *  :*  : V ■¥ -? : V -f  -t  '•+  :t  +  ■*  =*-ff  :+.*  *  +:* 


SUBROUT I  •  *r;  RLL'NC*;  Pf  Ofc‘*r» ,  W  > 

C  ■*:**  *  -I  *-* -i  :* »  :*  *  :*  *:*:*:*:*  :f 

c 

C  *  THIS  SUBROUTINE  CALCULATES  THE  INCIDENT  FIELD  ON  THE  OBSTACLE, 
C  aY  TkE  PAR.20NE  POINT,  AND/OR  AT  THE  CENTER  OF  THE  OBSTACLE . 

C 

C 

INTEGER  NCELLS,  NETS,  PTOBS,  I  PH  I ,  MCELLS 
C 

REAL  XNN,  VNN ,  ErtAG I ,  U1  ,  JO,KO,K2,  VO,  VI ,  EiTAGN , 

1  ECHislPU 
C 

COMPLEX  CNN ,  E I NC 1  ,£I nC  C< , E  I nC2 , 

)  ESC- AT, 

2  T  AU  1  .  A!_  MOD ,  AL  PMh  ,  COEF ,  FaC  ,  C 1  ,  C2 ,  V ,  V 1  .  VO , 

3  T  AU , V I N , VOUT , ONE , ZERO, A 0 , A  t 
C 

DIME NS I  ON  I EOF R< 1 S > 


C 


COMMON 


/BLKPiM/  AA , 8, C, R,F REQ, PERM, OGRENC, EPSLNR, AN, NCELLS , 
E I N  C  i  <  3  6  0  > ,  V  IN*;  l  0  5  i  > ,  E  S  C  AT*::  >  6  0  > , 


ECHWPLK 36 0 > , K 0 , R2 , VOUT < 1 S3 1  >,T AU< f  S3 1  > , 

E I NC2 , SCRPEiX  36 1  > ,  DBPRRT-:  J6i  A 0-:  f  S? f  > ,  A  I  <  IS 


EQUIVALENCE  <  NCELLS, N2 > 


DATA  P I  S3 .141  5S265/ 


C 

CUPENT  = 1  . 0 
C  * 

NP7  S- 1 F I )«  36  u . /DGRENC  > 


C 

C  *  WAVE  «  IN  FREE  SPACE  AND  IN  THE  OBSTACLE 
C 

K 0=2  .  *PI*FREGr/3 .  ES 
K2=k  Ci  i'SQRT <  EP SLNR  ) 


C 

C  ’*  EPSLNCi  =  S.S54E-12 
C 

RN*C~ <  W/2  .  :> 

THETAN=2 . *P I /NCELLS 


C 

C: 

C  * 
C 
c 
c 

c 

c 


DUMMY  t  *  -<  K 0*  +2  .  /C  4  .*2  .•*PI*fFRE€i:*8.S54E-i2  >  >*CURENT 

IF  OBSERVATION  POINT  AT  PARTICULAR  CELL,  PTOBS  =  0.  IF  OBSER¬ 
VATION  POINT  IN  FaR  ZONE  AT  SOME  ANCLE  PHI,  PTOBS  =1.  IF 
OBSERVATION  POINT  AT  CENTER  OF  OBSTACLE,  PTOBS  =  2. 

I  PC PTOBS  ,EQ,  1>  GO  TO  20 

1 FC PTOBS  -EQ,  2  >  GO  TO  SO 

DETERMINE  THE  INCIDENT  FIELD  ON  THE  OBSTACLE  AT  EACH  CELL 
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•J  IT*  O  O  *J  *-  ■—  O  O  O  *j  O  CM  L>  O  CM  n  O  O  O  O  If* 


<  y.n ,  Vn  :> 


C  L  UUa  i  i  ON  , 

C 

C  ■*  LOOP  NCELLS  TIMES  POP  INCIDENT  FIElD 
C 

DO  10  I  ~  1  ,  NCELLS 

THETA®*.'  I-i  >+THFTAN 
XNN=RN*C0S<  THETA  > 

YMM*RN*SIN<  THE  Th  > 

DUMMY  0~<  i  NN-G  .  YNN-0.  > 

R.RGH-KCi  i-SQPT-;  JiNN-AA  >*<  XH N-AA  >4  DUMMY  0  .» 

CaLL  BESEL*;  i.  ,  akGh  ,  Oci ,  Y  ,  BF*  ,  YP  ) 

UO-BB  « 

Y  (i=Y 
C 

Oft  I  TE>;  S ,  5  >  J  0 ,  Y  0 ,  £B 

FORMAT <  IK,  "FLDNC  :  J 0* “  ,  2?. ,  1 E 1 5 .  S ,  5X ,  .  2>: . 

1  "  BB“  "  ,  2N ,  t  £  1 5 . 3 ,  '■  > 

E I NC  0=DUf1r1Y  1  :»CMPL X*.  J 0 ,  - Y 0  > 

V I  N*;  I  >--£  I  NC  0 

WRITE1;  S,  1  00  >  EIMCO,  I  .  V I I  .* 

0 0  FORMAT*:  IX,  "El NC  0=  -  ,  2X ,  1  E  1  5 . 9 , 2X ,  E  1  5  .  3 ,  CM ,  "  I  =  -  , 

f  5X,  “VIM* “  ,  2X,  IEIS.S*, 2X,E  1  5  .  .* 

0  CONTINUE 

GO  TO  40 


4.  DETERMINE  THE  INCIDENT  FIELD  AT  THE  FaR  ZONE  POINT 
THE  CURRENT  FILAMENT  AT  RaU  PRIME,  USE  LARGE  ARGl 
ASYMPTOTIC  EXPANSION  FOR  THE  HANKEL  FUNCTION. 

0  DO  30  J*1,NPTS 

PHI®*:  J-1  >*DGRENC*PI/1S0. 

ARGH*K0--*<  R-AA*C0S',  PHI  j  >-PI/4  . 

D 1  =C0S<  ARGH  ■) 
ft£r-  -S I  NX  APGH  > 

D3=S‘QRT<  2  .  /«;  K  0**P  I  R  ;  > 

E I  NC  1  <  J  >* DUMMY’  I  h*D3**CMPLX<  D 1  ,  D2  > 

WRITE*;' B,  200  >  E1NCK 0  >, J 

00  FORMAT*:  IX,  “FLDNC  :  E IHC1  =  ”  ,  2X,  t  E 1 5 .  S,  2X,  E 1 5  .  S, 

1  14, /> 

0  CONT I HUE 

GO  TO  40 

*  DETERMINE  THE  FIELD  INCIDENT  AT  THE  CENTER  OF  THE  ( 

0  ARGH - K  0*ABS< AA  > 

CALL  BESEL<  0 ,  ARGH ,  BB ,  Y ,  BF’ ,  YP  > 

J0=BB 

Y0=Y 

C 

EINC2=DUMMY1*CMPLN<  JC,-Y0  > 


I  4  , 


DUE  to 
!ErlENT 


5X .  "  I  =  ** 


JBSTaCLE 


c 

4  0 


RETURN 


T  T  -  0  0  0  04  IS  ON*  CffOOC'VS  USING  OuiH  Ct  BLKS  R«tr«?»U 


FTM  4 ,  L 

i'fcf  JrK  BL  Knri ,  0  > 


C 

C 

c 

c 

c 

c 

C 

C 

C 


c 

c 


c 

c 

C: 

c 


:f. :?  :*  :f .*  : *•  .4.  '■¥ :* :* :*  :*:*  4*  *  :*  :*  *  :*:*  :f  :*  :*  :*  :*  :* : fr  •+ :*  ■*:¥ 

5UBRGUT  I  HE  E  00*7 \  Xri«^  > 

*  rHI-S  SUBROUTINE  PLOTS  THE  NORMALIZED  ECHO  u-XOTH  PEP  WAVELENGTH 
VS.  ANCLE  PHI  ON  A  LINEAR  PLOT . 


INTEGER  NCELLS’,  NP ;  S,  PTC:£:S,  IPHI,nCELLS 

REAL  XNN,  VNN/ENAGI  ,  J1  ,  JO , KO,  K2 ,  YO,  V 1  ,  EftAGM , 

J  ECHitiPw 

COMPLEX  CMM , E I NC 1 , E I MC  0 , E I NC2 , 

1  ESCAT . 

2  T AU  I  ..  ALMDa  ,  ALPHA ,  COEF ,  F AC ,  C  )  ,  C2 ,  V ,  V  1  ,  Vi , 

3  T  AU , V I N , VOUT , ONE , ZERO , A  0 , A 1 

D I MENS I ON  I BUFR< 1 6 > 

COMMON  SBLKMS  AA ,&,  C ,  R ,  FREO,  PERM..  OGRENC,  EPSLNfc ,  AM,  NCElLS 

1  E I  NC  1  SCO  >,  V I  Nr.  1  $51  >,  ESCMT<  3  SO  ■ . 

2  ECrtWF’WC  33  0  ') ,  K  0 ,  K2 ,  VOUTC  !S5(  )  ,  TAUC  NS51  >, 

3  EINC2,SCRP£2< 361  >,DBPRPTC361  >,  A0<  i  05 1  >,A1<  1 

EOU I VALENCE  < NCELL S , NZ  > 


DATA  P I X3 . 1 4 1 5S2S5/ 
NX  Lip =2 .  *  1  000 
NY  OP*  1  .?£«■'*  1  000 

* 

NXST P= 0 . S*  J 0 0 0 
NYSTP* 0 . 5* 1 0 0 0 


C 

NXSTPS=S 
HYSTPS*=1  0 
C 

IPEH= I 

WRITEC 25, 3 >  I PEN 
C: 

NX£P»NXOP+NXSTPSm.NXSTP 
WR I T  EC  25 , 1  >  NX OP , NY  OP 
UR 1 TEC  25 , 1  >  NXEP , NY OP 
UR I TEC 25, 2 > 

C 

NYEP*NYOP+NYSTPS*NYSTP 
UR I TEC  25 , 1  >  NX OP , NY OP 
WR I  TEC  25 , 1 >  NX OP , NYEP 
UR 1  TEC  25,2  > 

C 

DO  100  0=1 , MAST PS 
IX=NXOP+J-*NXSTP 
IY(=NYf<P-50 
IY2-HY0P+50 
UR  I  TEC  25,1  >  IX, I Y 1 
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!*l  M  —  r<  *.  r.  r.  o  r>  |*j 


W&l  ?E*t  25, 1  >  IX,  TV2 
2  Tt<  5zT» /  .> 

I  0  fi  CO  til  1  ti'JE 

C 

C>U  200  1  =  1  ,  NYSTPS 
1 Y=HY 0P+ I *HYSTP 
1X1 =NX 0P-50 
IX2=r!X0P+50 
WRITE*; 25,1  >  I X  1  ,  IV 
WRITE*.' 25,  1  >  1X2,  IV 
WF  1 7£*.'  25 , 2  > 

00  CONTINUE 

CALL  ECLBL 

CALL  LABEL*;  1  > 

NP  TS» I F I X< 1 3  0 . /DGRENC > 

IP£li-3 

Wk  I  it<  25 , 3  >  I  PF.N 

DO  *5  v  0  I  PH  I  =  1  ,  ARTS 

I I F  I X-;  •;  I RH I  - 1  .  >  *DGRENC.  "2  0 . 
i  -  I F I  <  EC  HWRuK  I  PH  I  VXMAX  .V 
wpi  rE*;25,  i  >ix,  iv 
00  continue 


WRITE*;  25,2  > 

FORMAT*.  “Pa"  ,  15"  ,  “  15,  “  ,*PD“  ;> 
FORMAT*;  "PU“  > 

FORMAT'.  “  IN;  SP"  ,  1 1  ) 

END 

END* 


:*NXSTR+2  0  0  0  .  > 
‘0.1  ;*NYSTF*  1750. 
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o  o  o  o  o  o  o  o  ooo  uooooo  o  o  o  o  o 


T 


7  7=  CO 004  IS  ON  CR0003S  USING  00004  ELKS  F:=  0 0 0 0 
FTN4 ,  L 

iENi-K  BLKMM ,  0  > 

C  ♦:m-w*-H‘<i****:*w*»»**»*»*»:t 
SUBROUTINE  FLDCT 
C  ■*'■*  +  * if  *♦  *■  ***■*■*-« 

C 

C  *  THIS  SUBROUTINE  USES  THE  MOMENT  METHOD  TO  DETERMINE  THE  SCATTERED 
FIELD  FROM  SOME  DEFINED  OBSTACLE 


INTEGER  MCELLS , NPTS , FT DBS , IFHI , MCELLS 

REAL  XNN , YHN , EMAG I ,  J1  ,  00, k  0, K2, YO, Y1  ,EMAGN, 
1  ECHWPW 


C  OrlPLEX  CMH ,  E I NC 1  ,  E  J  hC  0 ,  E I NC2 , 

1  ESCAT , 

2  TAU 1  ,  ALMDA ,  a.LRHh  ..  C OEF ,  FAC ,  C  i  ,  C  2 ,  V ,  V 1  ..  V2  . 

3  TwU  ,  V I  M,  VOUT  .  ONE ,  ZERO ,  DUMnYS ,  A  0 ,  A 1 

DIMENSION  IBUFRC 1  £  > 

COMMON  /BLKMM/  AA, 3( C, ft, FREQ, PERM , DGREHC, EPSlNR, AN, MCELLS, 

1  E I  NC  T  <  36  0  > ,  V I  N>;  )  S3  i  > ,  ESCaT r.  360  > , 

2  £CHuiRu(<  3bCi>,KC<,K2,  VOUT*!  1*5 1  >,TAU<  ’S5i  >, 

3  E I NC2 , SCft‘Pfc4<  36 1  0, DE*rRR  i'i  4o1  ) , fi U\  l  &  5 1  ,',A1*.  1*c’5i 

EQUIVALENCE  <  NC  ELLS',  N2  ) 


DATA  P 1/3. 14  159265/’ 


*  THE  LARGE  ARGUEMENT  ASYMPTOTIC  EXPANSION  FOR  THE  HANKEL  FUNCTION 
IS  USED.  SEE  RICHMOND  - 


*  LOOP  NPTS  TIMES  TO  OBTAIN  THE  SCATTERED  FIELD  AT  EACH  DGREHC 

NP7S=  I F 1  X<  36  Ci .  /DGREHC  > 

RN=C~< W/2 .  } 

THETAH=2  .  *P  I  /NCELLS 

DO  20  IPHI*1,NPTS 

PH1=<  IPHI  — 1  >* DGRENC*P I / 1 S 0  . 

DUMMY  t  =<  K  O'*  ft  P I  /4  .  > 

DUMMY2=—  COSC DUMMY 1  ) 

DUMMY3=-S I M< DUMMY 1  ) 

DUMMY4*SQRT<  PI=*K  0*  0 . 5/R  > 

*  LOOP  MCELLS  TIMES 

DUMMYS=< 0 . , 0 . ) 

DO  10  IT*  1,  NCELLS 

THETA=C  H- 1 >* THE TAM 
XNHaRN *COS<~  HE7 A  > 

YMH»RM*SIN< THETA ) 
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£>UHr)rt»*i< u+:  XHK ‘ f Phi  >*WN+51Mi  PHI  »  > 
Dc.JhMVt.-si'.os';  £-ur'»~,v5  > 

DUMMN  7 -■£  1  Hi  DUMMYS  > 

ARGJi  =kO-<,rtM 

ChLL  BE;  EL<  1  ,  t-iRCi  J  l  , BB , Y , BP ,  YE  > 

J 1  =BB 

DUMMYS=DUMMYS+<  EB'j'LHR- 1  .  >*VOUT<  N  >*AM*  J 1  *  CHPUX 
CONTINUE 

ESCN  Ti  I  PHI  .>-DUHrtY4*C  HPLXi  DUMMY3,  DUMMY*  MDUmMYS 
CONT I HUE  • 

RETURN 


•_)  O  |J  ‘J  O  ■_)  O  O  O  O  O  O  O  O  *  J  G  O  O  G  G  G  O  G  o  o  o  o  o 


T=00004  IS  OH  CR0003S 


OS  I  HO  OOOOo  ELkS  P=  GOOD 

FTH<5,L 

lEriw'.'  rlkmm»  o  > 

SUBROUTINE  TPi.E>:'  r!M,  XNORM,  I£R> 

C  ■»■«■■•■*'»•:»  *:**:*:*:***:*:*:*ts*:*:^*:*:Wr*'*:|!*4 +*:**:**:*•*  t.:* ;*:*  (. :f  »*  (, 
C 

c 

C  F  r  on. 

C 

C  "Antenna  Theory  end  Derigh" 

C  ijlsi'i'tn  L.  Stutsman  sod  Cary  a .  Thiele 

C  John  Wiley  j.  So ;.s ,  Hew  York  ,  19S  J 

Append  i  x  0,7  pp  .  57S-5S  t 


■*  PURPOSE 

TO  SOLVE  A  SYSTEM  IHVuLVIHO  A  TCEPLITZ  MATRIX.  TPlZ  PEQUIFt 
ONLY  5 H  STORAGE  LOCATIONS  FOR  AN  H  BY  N  MATRIX. 

*  REMARKS' 

A  tocpl itz  matrix  has  the  first  row  equal  to  the  first  column. 

All  elements  along  the  main  diagonal  are  egwai .  Any  diagonal 
off  the  main  diagonal  will  have  This  same  property. 

*  DESCRIPTION  OF  PARAMETERS 
NZ  -ORDER  OF  MATRIX 

TAU  -FIRST  ROW  OR  COLUMN  OF  THE  TCEPLITZ  MATRIX  <  VE07 OR 
LENGTH  N2> 

A0,A1  -VECTORS  OF  LENGTH  NZ  HEEDED  FOR  SCRATCH  AREA 

VIN  i-FOR  THE  MATRIX  EOUATION  <ZXI)=<V),  VIH  IS  V.  <2,1, 

AND  V  MAY  BE  THOUGHT  OF  AS  GENERALIZED  IMPEDANCES, 
CURRENTS,  AND  VOLTAGES,  RESPECTIVELY  > .  V  IS  A  HZ  BY' 

MM  MATRIX. 

MM  -NUMBER  OF  COLUMN  VECTORS  ON  THE  RIGHT  SIDE  OF  MATRIX 

EOUATION  <Z><  I  >»<V>  *:  USUALLY  1  >. 

XNORM  -UPON  RETURN  THIS  IS  INFINITE  MATRIX  NORM  OF  INVERSE. 

IER  -ERROR  CODE  WHICH  IS  0  IF  NO  TROUBLE. 


INTEGER  NCELLS, NPTS, PTOBS, IPHI.MCELLS 
C 

REAL  XNN , YNN , EMAG I , J1 , JO, KO, K2, YO, Y1 , EMAGH , 

1  ECHWPW 

C 

COMPLEX  CMN , E I NC 1 , E I NC 0 , E 1 NC2 , 
i  ESCAT, 

2.  TAU  l  ,  ALMDA,  ALPHA,  COIF,  FAC,  Cl  ,C£,V,V1  ,  V2, 

3  TAU , V I N , VOUT , ONE , ZERO , A  0, A 1 

C 

DIMENSION  IBUFR< 1 6  > 

C 

COMMON  /BLKMM/  AA ,  B ,  C ,  R ,  FREO ,  PERM ,  DG-ENC ,  EPSLHR ,  an  ,  NCELLS  , 

1  E I  NC  1  <  36  0  > ,  V I  Nr!  1051  >  .  ESCaT <  36  0  >  . 

2  ECHWPuK  36 0  >,  K  u ,  Ki ,  vCUT <  1  05  1  >  .  T 1  S5  f  > , 


O  O  O  O  O  O  O  O  O  O  10  O  r.  r»  —  — .  r ,  ^  o  O  O  f> 


E I NC2 , SCRPEZC  36 i  > , OBPEP T *.  36 1  >.  A 0-.  1851 
ECU  I  VALENCE  <  isCEl.LS,  HZ'  > 


DATA  OHfc/C  1  E  0 ,  OE  0  '/ ,  ZERG/C OEO  ,  GEO)/ 

DATA  GNNE/1E0/,  2RRG/0.EG/ 

WRITE*:  6, '5 Co  NZ,MM 

FORMAT*:  1 X ,  “NZ*=",2 X,  14. SH.  “MM=“,2X,  14,/) 

DO  150  1=1, HZ 

WR I TEC 6, 1 uu  >  T AUC I > , V I N< I >,  I 

FORMAT*;  1  “TPL2  :  TAU=",2X,  1  El  5 . 8, £X, E 1 5 . 8, 5X,  "  V I  N=  "  ,  2:< ,  1 
1  2X,E15.S,5X,  “I**1,  14, /> 

COHT I HUE 
N=HZ-t 
I  Ek=  0 

■*  NORMALIZE  INPUT  MATRIX 

TAU1=TAUC  1  :> 

DO  2000  11=1, N 

0 0 0  T  AUC 1 1 >=T AUC 1 1  + 1  >/TAU 1 

*  THE  FOLLOWING  CALCULATES  THE  ITERATIVE  VARIABLES  TO  OBTAIN 

AOCH)  AND  ALMDA 

*  NOTE  :  VECTOR  AOC  I  :>  HAS  I  ELEMENTS  AND  IS  STORED  hS  A0*:i,J), 

J=  1  ,  N 

ALMDA=ONE-TAU< 1  )**TAU< 1 > 

.  AOC  1  >=-TAU<  1  > 

v  1*2 

1  KK=I— 1 

ALPHA=ZERO 
DO  2  N= 1 , KK 
LL=I-M 

'  2  ALPHA* ALPHA* AOC  M  >***TAUC  LL  > 

ALPHA*-*;' ALPHA+TAUC I >> 

IFC  CABSC  ALPHA  >  ,EC.'.  O.DO>  00  TO  15 
C'OEF=ALPHA/ALMDA 
ALMDA=ALMDA-COEF:*  ALPHA 
DO  3  V= 1 ;kk 
L=I-J 

3  A 1  <  J  )=A  OC  J  )+CGEF*«*A  OC  L  > 

DO  7  J=  1  , KK 

7  AOC  J  >=A  1 C  0  > 

AOC  I  >=COEF 
IFCI  ,GE.  N>  GO  TO  5 

4  1=1+1 

GO  TO  I 

+  THE  FOLLOWING  COMPUTES  VALUES  OF  EACH  ELEMENT  OF  THE  INVERSE 


0 

00 

50 


NH=C  NZ+ 1 >/2 
FAC=ALMDA*TAU1 
XHORM=ZRRO 
MP=N2+1 
DO  5»  1*1, NH 
XHM=ZRRO 
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<_>  Cl  O  «'  «|  1_>  lit  O  O  '  £>  'i>  •_> 


IF*  I  . MS'  .  1  'J  CO  TO  52 

A1  *  1  >=OmF.  T-'i  .0 
Xnr=CaBS*  hh'I  >> 

DO  53  J=2,N2 
ai*  o>=ao*  J-1  >/FAC 

53  XhM=CASS'<  HI*  J  >  >+XNrt 

00  TO  54 
51  XNM-ZfcRO 

UH~ 1  - 1 
Cl*AG<  JH  ) 

NNPI*NP-1 

C2*A0*HNPI>  • 

DO  55  JJ*1,N 
J -HF-JJ 
IMP J=NP-J 

\JL=vl-1 

hi*;  j  >=h  i  >;  jl  >+*  c  i  *  a  o*  jl  >-c2* a  <k  i  up  j  >  >/f  ac 

55  XNM-CABS*  A 1  >:  J  >  >+HNM 

A  I  *  1  >=A0*  I  -  j  >/FAC 

xnm-^nm+cabs* a i *  1  >> 

54  I F<  XNM  . GT .  XNORN  >  XNORM*XMM 

c 

C  «  MATRIX  MULTIPLY 
C 

DO  56  1 1= 1 ^ MM 
ID**  1 1-1  >*N2 
V*Z£RC< 

VI  =2ERCi 
DO  5?  J= 1 , NZ 
HJOJ-I D+  J 
V2= V I N<  H I D J > 

V*V+V2*A  KJ) 


57 


C 

225 


56 


50 


51 


KNPJ*MP-J 
V 1  =  V 1  +V2*H  1  <  MNP  J  > 

N I D I -  I D+ 1 
VOUT* MIDI  >=V 

WRITE* 6, 225)  VOUT*  N I D I > , N I D I 

FORMAT*  IX,  "TPLZ  :  VOUT-" , 2X, 1 E 1 5 . 8, 2X, El  5 . 3, 5X,  “NI«",2X,  I4,/:> 
K I DHP I  *=  I D  +NP- 1 
VOUT * K I DUP I >=Vt 


WRITE* 6, 250  >  VOUT*  kIDHPI  >,  KIDHPI ,  VOUT*  NIDI  >..MIDI 
FORMAT*  1 X ,  -  TPLZ  :  VOUT-  “  ,  2X ..  1 E 1 5 . 8 , 2X ,  E  1  5  .  & ,  5X ,  “  K I  *  "  ,  2X ,  1 4 , 
3X, "V0UT*",2X, 1 E 1 5 . 8, 2X, £  15.8, 5X, "NI-“,2X, 14, /> 

WRITE*  6,251  >  VOUT*  K I CNR I >,KIDNPI 
FORMAT* IX," TPLZ  :  VOUT  = “ , 2X , 1 E 1 5 . 8 , 2X , E 1 5 . 8 , 5X , “ K I  =  “ , 2X , 1 4 , / > 


1  COMT I HUE 


00 

50 


DO  650  1  =  1  ,  NCELLS 
WR I TE*  6 , 6  0  0  >  VOUT* I > , I 

FORMAT* IX, "TPLZ  : :  VOUT*" , 2X, 1  El 5 . 8, 2X, E 1 5 . 8, 5X,  " I* “ 
CON7 i HUE 


2X , I 4 , / > 


RETURN 

C 

15  WRITE* 6, 700 > 

700  FORMAT*  IX,  “ERROR  HAS  OCCURRED.  MATRIX  IS  STRONGLY  NONSINGULAR'1  > 
I  ER=  1 
C 

RETURN 


yi  o  oi  O’  o  o  o  o 


£>  I  MS  HS I  ON  C  0<  ?  > ,  C 1  <  7  > ,  D  O';  7  > ,  DU?  > ,  E  0\  7  > ,  S  1  *.  7  >  .  C  0<  7  '> ,  CK7  > 
DATA  CO/1  .  0,-2.2499997,  1  . 2656200 .-. 3 1 63066 ,  .444479E-1  , 

-  .  39444E-2,  .  2 1  £-3/ 

DATA  Cl/ .  5,  -  .  56249905,  .21  093573,  - .  3954209E- 1  ,  .443319E-2, 

*- .7 17c* IE-3,  .  1  1  09E-4/ 

DATA  D  0/ .  79  »■  00456,  —  .  77t  — 6  ,  —  .  55274E— 2,  —  .  95 1  2E— 4  ,  .  i  372'37*'E— 2 , 

•*  -  .  720  05E-03,  .  144 70 E  -  03  / 

DATA  D1 / . 79700456,  . 156E-5,  . 1659667E-1 ,  . 171 05E-3 , - . 2405 1  IE-2 
'*,11 3653E— 02,  -  .  2 n 033S—  03/ 

DATA  E 0/-  .  7053901 6,  -  . 4  1 0639 7E- 1  ,  -  .  3954E-4 ,  .  26257 3E-2, -  .  54  1  2 
. 29333E-03, . 1355BE-03/ 

DATA  E1/-2. 356 19449, . 12499612, . 565E-4 , - . 637079E-2 , .7434EE-3 
.  79S24E-03 ,  —  .  29 1  66E— 03/ 

DATA  CO/ .  3674669,  .  60551937,  -  .  7435030,  .25300 1 2, -. 4261 2  IE-01  , 

.  42791 6E-  02,  -  .  24846E- 03/ 

DATA  Cl/-. 6366190, .2212 091,2. 160271,-1 .316403, .312395, 
*-.40G976E-01 , .27O73E-02/ 

DATA  01/3. 1415926/ 

IFi_G=0 

IF  <ii.LT.CO  IFLG=i 
H=IA66< M  > 

IF<  ABS< ii  >  .  LT  ,  1  .  OE- 1  0  >G0  TO  1  5 0 

I  F<  ABS<  X  .  GT  .  3 . 0  >G0  TO  50 

X30C'=X+X/9 .  0 

PR 00= 1.0 

B0=0. 0 

B1  =  tv  0 

CUM 0=0. 0 

CUr11  =  0 . 0 

SEE  EQUATIONS  9.41  AND  9.44 

DO  5  1=1 ,7 

B0=B0+C:0<  I  >*PROD 

B 1  =B  1  +C 1  <  I  >*PROD 

CUM 0=CUH G+C G-;  I  >*PROD 

CUM  1 =CUM 1 +C 1 < I >*PROD 

PR0D=PR0D«X3SG 

CONTINUE 

B1=Et1*X 

XC=2 .  0*SNGL<  DLOG<  DBL£<  0 . 5*X  >  >  >/PI 
Y  0=XC*B  O+CUM  0 
Y1=XC*61 +CUM1/X 
GO  TO  100 

EOS  9.4.3  AND  9.4.6 
0  THR0VX=3 . 0/X 
PROD® 1 . 0 
F  0=  0 . 0 
F 1  =  0 . 0 
THETA C=X 
THETA 1=X 
DO  55  1=1,7 
F  0=F  0+D0<  I  >*PEOD 


F  !  =F  1  *Z>i<  i  >*-RROD 
THETA  0*TKET  »*£•»•£(  \  1  >*P-:D 
THETA 1  =THETrtl  I  >»F>:CO 
PkOD=PROD*  T  hF:OVX 
55  CONTI HUE 

SffR.x- i .  o/sngl*;  dsurtc  dsle*;  x  >  >  > 

3(i-SCjRX!*F 0;*-SNCL>.  BCOSC  i'iL£<  T NET AG  '»>> 

B I  =SC'RX*F  1  *SHGL<  DC  OS*;  R£lE*;  THET At  >  >  > 

YO=SC.*F;X*F 0*SNC-L<  DSIN*;  DBLE*;  THETAu  >  >  > 

Y1=$GRX=*F  t  *$USL<  DSINC  DELE*;  THETA  1  >•>  > 

1  00  IF*;  N-l  >101,1  05,  1  1  0 
101  B=B0 

BP=-B1 
Y=Y  0 
YP=-Y 1 
CO  TO  200 
1 05  B-B 1 

BP=B0-B1 XX 
Y-Y  i 

YP=Y0-Y 1 /X 
CO  TO  200 

C  FOP  RECURRSI VE  DIRECTION  COMnENTS  'SEE  SECTION  B.12,P3S5 

1  1  0  XH--H 

I F< XN .  L T  .  ABSY  X  >  >00  TO  130 
C  FOR  X<  N  RECUR  DOWNWARD 

6L«ST= 1.0 
BLASTP--G. 0 
J=N+ 1  0 

DO  115  1=1, J 
XI  =  v*-I 

BN EXT=2 .  G*X I *BL AS  T/X-BL ASTP 
&LASTP=BLAST 
BLUET  =BHEXT 
IF< I .HE. 1 0>CO  TO  1  15 
BLF -BLASTP 
2=BHEXT 
115  CONTI HUE 

IFC ABSC B0>.LT . ABSCB1  >>CO  TO  117 
CORE - B  0/BLAST P 
CO  TO  IIS 

117  COF:R=-B  1  /BNEXT 

1  1 S  6*BLP*C0Rft 

BNM I M 1 =2*C0RR 
GO  TO  140 

C  FOR  X >H  RECUR  UF’UARD 
130  BLASTP-BO 
BLASTER  1 
DO  135  1=2, N 
XI=I-1 

BNEXT  =2 . 0#X  I  *BLAST/X-&L ASTP 
BLASTP=BLAST 
BLAST-BNEXT 
135  CONTINUE 
6= BLAST 
BMNIMI = BLASTP 
140  BP«=BNNIN1-XM*B/X 
YLASTP=Y0 
VLAST“Y I 
DO  145  1*2,  N 
XI*I-1 
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